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NOISE STORAGE IN DELAYED 
FEEDBACK DEVICES 


Yu. S. Lezin 


The noise power at the output of delayed-feedback storage devices is determined by the 
"slowly varying magnitudes" method. 


INTRODUCTION 


The most important element of optimal filters for sequences of pulsed signals [1] is the 
summer with a feedback circuit having a delay equal to the repetition period of these signals. 
Equality of the repetition period and the delay time in the feedback circuit can be attained 
with sufficient accuracy for practical purposes by use of the methods employed in systems 
for moving-target indication [2]. However, this kind of summing device is self-exciting be- 
cause its positive feedback coefficient is unity. 

A practical approximation to such a device is a memory with positive feedback and a 
delay time T. For stable operation the memory's feedback coefficient should be less than 
unity. Figure la gives a block diagram of the memory, consisting of a summer, an attenuator 
with transfer coefficient m (equal to the feedbackcoefficient at sufficiently low frequencies), 
a delay line with delay time T and a filter F with transfer function 
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where 6 is tne time constant of the filter and AF is its bandwidth at the 3 db point. This trans- 
fer function coincides (to within a constant factor) with the frequency characteristic of an un- 
compensated resistance coupled amplifier in the high-frequency region. 

All the elements of the memory besides the filter have uniform amplitude-frequency 
characteristics extending at least up to the frequency at which the modulus of the filter's 
transfer function becomes negligibly small. Let us assume that all the elements are linear. 


Figure 1. Block diagram of a delayed feedback memory. 
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A memory for pulsed signals in such delayed-feedback devices has been thoroughly 


analyzed in Ref. 3. : 
In the present article we consider a memory for conventional white noise with a power 


spectrum 


F,(o) =a. 


The results obtained are also accurate for the passage of conventional white noise through 
a filter F and subsequent storage in a delayed feedback device (Fig. 1b). 


1. NOISE STORAGE IN A SINGLE MEMORY WITH DELAYED FEEDBACK 


The transfer function of such a memory is given by 


K* (@) K (@) ea farce tg wb 


K;() =- gas ae eN 
(0) 1—mK*(w)e—7@F 4 — mK (w) e774)” 


and the square of its modulus 


KE? (@) 


2 
Kj (@) = 1 —2mK (w) cos (@) + m?K? (@) ’ 


where 


K 0) Se (1) 


is the modulus of the transfer function of the filter F; 
p(@) = oT + arc tg w6. (2) 


Therefore, the stored noise in the memory with delayed feedback has, at the output, the power 


=| Fo) Ki (w) do ay 


0 0 


K? (@) dw 
1 — 2mK(q) cos tp (@) + m?K? (w) * (3) 


Evaluation of this integral is ordinarily quite difficult. Therefore, we make an approxima- 
tion, using the circumstance, peculiar to the devices we are considering, that the time constant 
of the filter is very small in comparison with the delay time: 


6<T. (4) 


Thanks to this we can apply the method of slowly varying magnitudes [4] to the evaluation of 
Eq. (3). It is apparent that the second term in the denominator of the integrand contains a 
rapidly oscillating component cos ¥ (w) = cos (wT + arc tg w8), but the remaining terms of the 
integrand, in comparison, change slowly with frequency. On a change of frequency equivalent 
to one period of the oscillating component, the remaining terms in the denominator and nu- 
merator are practically unchanged. We break up the integration in Eq. (3) into parts, each of 
which corresponds to one period of the oscillating component: 


2 K?* (@) do 
3 =a >) 1—2mK (o) cos tp (@) - m®K2(0) ’ 


l=0 @] 


where w, and w,,, are the frequencies corresponding to the beginning and end of the 
l -th period. 
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(or) = 20d, PH (W144) = 2m (1 + 1). 


As a consequence of inequality (4) the value of this period is approximately equal to the re- 
petition frequency. 


Q) = O41 - 0 ~Q =, 


We make a change of variable under the integral sign and let x = w — u,. 


K? (@, + x) dx 


As a consequence of inequality (4), for 0<x< Q) we have 


K (@ + 2) ~ K (@)), 


cos ) (@, + x2) ~ cosaT 


and therefore, 


2 

K? (w,) dx 
\ 1 —2mK (@,) cos-«T + mK? («)) = 
0 


Dae C d 
~s 2 OA de Ses ere a 
=p > *K (a) | 1— 2mK (@) cosy + PRK? (@) 
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Using Eq. (3.451.2) of Ref. 5 we obtain 
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This result can also be obtained by the stationary-phase method [6]. 

Since the magnitude of the 7 -th period of the function cos %(u) in view of the inequality (4) 
is extremely small in comparison to the width of the passband, and consequently, also, in 
comparison to the width of the noise spectrum at the output, then the last sum can be approxi- 
mated by the integral: 


oh =| = RO) aes (5) 
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The integrand a ea can be considered as a power spectrum of the stored noise, 


averaged over a period of the oscillating factor cos (wf + arc tg wd). Substituting Equation 
(1) into Equation (5) and evaluating, we obtain 
2 waAF (6) 


Os == 


Thus, the power of the stored noise is proportional to the intensity of the noise at the 
input and passband of the filter F. In order to determine the power of the unstored noise, 
that is the noise arriving at the output directly, bypassing the feedback circuit, it is sufficient 
to set m = 0 in the last expression: 


6x) = WaAF, 
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Therefore, the coefficient of noise storage in the delayed feedback memory, in terms of 
power, is equal to 


Ue ioe ; 


If a rectangular approximation is made to the amplitude-frequency characteristic of the de- 
layed feedback loop of the memory, i.e., if we consider the filter to be an ideal low-pass 
filter, then the noise storage coefficient is seen to be equal to [7]. 

i 1 


oO 
& 1 —m?° 


The dependence of the noise storage coefficient on the magnitude of the feedback coefficient, 
both in the case considered and in an ideal filter, are presented in Curves 1 and 2, re- 
spectively, of Fig. 2. Noise in the real memory investigated is accumulated to a much 
smaller degree than in an ideal memory with a feedback loop having a rectangular amplitude- 
frequency characteristic. This is explained by the narrowing of the power spectrum of the 
noise when it is circulated in the feedback loop 
of a real memory. In an idealized memory the 
noise power spectrum is uniformly attenuated 
only by a factor of m each time it is circulated, 
and its width remains unchanged. 


Figure 2. Noise storage coefficient (1) for a 
single memory (2) with an ideal filter, as func- 
tions of the feedback coefficient. 


2. NOISE STORAGE IN TWO AND THREE SERIES-CONNECTED MEMORIES 
WITH DELAYED FEEDBACK 


The noise at the output of two series-connected identical memories, each of which is 
encompassed by a feedback with a delay time T, and arranged as in the block diagram pre- 
sented in Fig. 1 has the power 


eet : e K4(o) do 
33 = \ Fi (@) Ki(o) do = a\ (1 — 2mK (@) costp (@) + mK? (EE 
; 


0 ( 


For an approximate evaluation of this integral we average the comb power spectrum of the 


output noise over a period of the rapidly oscillating factor cos ~(w). As a result of this we 
obtain 


co 
{1 ++ m?K? (@)] K4 (@) 
Sa TEE CTC) CAC 


Substituting Eq. (1) into the resulting expression, transposing, and applying Eq, 
(3.191.3) of Ref. 5 we determine the power of the stored noise 


am@AF 2+ m? 
~_ 


ee a m?)'/2 } ) 
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and then also the noise storage coefficient of two memories with delayed feedback, 
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ees 
Ug ee al ae (8) 
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This is less than the noise storage coefficient of two memories with rectangular 
amplitude-frequency characteristics in the delayed feedback loop: 


A 
i {me 


g" =~. 
S (1 — m?)8 


The dependence of the noise storage coefficients of two memories on the feedback coefficient 
for the case considered and for rectangular amplitude frequency characteristics are pre- 
sented in Fig. 3 (Curves 1 and 2, respectively). The ratio of the noise storage coefficients 
for one and two memories with delayed feedback is equal to 


y= 2 — 2m? 
2 


i f 2(1— mm)?" 


The dependence of this ratio on the feedback coefficient is given by Curve 3, Fig. 3. This 
ratio greatly exceeds the noise storage coefficient for a single memory. The explanation of 
this is that on passing through the first memory the noise becomes noticeably correlated, and, 
because of this, storage in the second memory behaves differently and thus, reaches a much 
greater level. 

In a similar manner we calculate the noise power at the output of three identical series- 
connected memories with delayed feedback. 


3sam®AF 8 + 24m? + 3m4 
64 (1 — m2) 


(9) 


Si = 


and the noise storage coefficient for the three memories 


i 8 + 24m? -- 3m4 (10) 
oA (i) | 8. 2n2) 2 : 


The dependence of the latter on the feedback coefficient is shown in Curve 4, Fig. 3. 

The ratio of the noise-storage coefficients for two and three memories with delayed feedback 
is equal to 

8 + 24m? + 3m4 

g 4(2 m*)(L1— ny 


and is shown as a function of the feedback co- 
efficient in Curve 5, Fig. 3. 

From a consideration of the curves given, 
it is possible to conclude that the ratio of the 
noise storage coefficients for two and three mem- 
ories sharply increase as the feedback coefficient 


Fig. 3. Noise storage coefficients for two and 
three memories and the ratio of these coefficients 
as functions of the feedback coefficient. 


approaches unity and even at m = 0.9 are 291 and 
142,600, respectively, which is a consequence of 
the strong correlation of the noise of the two 
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preceding memories with delayed feedback. This correlation increases as the feedback 
coefficient approaches unity. 


CONCLUSION 


The noise power at the output of one, two and three memories with delayed feedback 
and frequency characteristics in the feedback loop coinciding with those of an uncompensated 
resistance coupled amplifier in the high-frequency region have been calculated by the "slowly 
varying magnitudes" method. Applying the method described we may make calculations for 
delayed feedback memories utilizing filters having a different frequency characteristic, and, 
also, for systems involving frequency filters and delayed feedback memories. j 

The calculations show that noise is stored in the delayed feedback memories considered 
to a lesser level than with ideal filters with rectangular amplitude-frequency characteristics. 
This is explained by the fact that the power spectrum of the noise contracts each time it is 
circulated in the feedback loop, when the filter possesses non-uniform amplitude -frequency 
characteristics. 

Noise stored in the second, and especially, in the third memory behaves differently 
than in the first, and reaches an extremely high level, due to the strong correlation it ac- 
quires on passage through the first and second delayed feedback memories. 
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LINEAR CORRECTION OF TRANSIENTS 
IN SYSTEMS WITH FINITE PASSBAND 


P. Ya. Nudel'man 


The report discusses the problem of decreasing the transient phenomena caused by the 
passage of a pulse through a linear system with a finite passband. 


Let us examine an ideal filter 


1, |@|<, 
We shall assume that phase shift is zero, which does not diminish the general applicability 


of the discussion. 
If a pulse ~(t) is applied to the input of such a filter, on the basis of the Duhamel formula 


we obtain at the output 


T/2 ; 
NO te ee ere TU (2) 
—T/2 : 


Here (—T/2, T/2) is the duration interval of pulse p(t). All values of f(t) corresponding 
to instants of time outside the interval (—T/2, T/2) will be known as after-effect: 
Let us now determine the form of the impulse (t) of T duration having a minimal after- 


effect energy: 


co T/2 
Wae = W—o.c0— Wem, m= \ peod— \ Pode (3) 
—0 —T/2 


in passage through an ideal filter. 
The total energy and the energy in the interval (—T/2, T/2) of the function f(t) at the 


filter output are, respectively, 


co 


, = { Asinau—9 (t) ac] dt = 
Woman \[ | 


™ t—T 
—o —T/2 
Dip Sap (ore) : 
4 sinQ(t —7) sinQ (t —s) Aes 
= FF oiiew | Hate nated ardcas 
—Th—Tp =) 


T/2 T/2 F 


—T/2—T/2 
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dt dt ds. (5) 


T/: 

= \ sin Q(t — T) sede 
t—t t— 

—T/2 


os fc 
W -—T/2, T/2 = \ 
Tr 


In order to find the pulse yielding a minimum in Eq. (3) let us investigate the functional 


Tj, T/2 
‘ ee Ao 
U [9] = \ \ @ (t) p(s) = = — ) duds 
—T/2—T/2 
T/2 T/2 f , T/2 os 
| si = Q(t 
—i \ \ P(t) 9 (8) = ( — = a ane —9) dt dvds. (6) 


—T/2 —T/e Te 


It is necessary to find the function (henceforth designated as ((t) for which (see, for 
example, Ref. 1) 6 U[@]=0 with the smallest value of the parameter A =A Win: 


Introducing the arbitrary function ¢ (t) and the parameter ¢, we obtain the following 
equation for the determination of A pw jn and P(t): 


ZUG O +e Wl, =0. vs eT 7) 


Equation (7) converts to the form 


T/2 


( PT 7 Pe 


iv 


T/2 
1 sinQ(t—s) sin Q(s—t) = 
3 \ pa Fee dsdt 0. (8) 


Equation (8) is satisfied only by the eigenfunctions of a homogeneous integral equation of 
the Fredholm type with symmetrical positively defined kernel: 


eS ) E(x) ae (9) 
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ee 
et) 
ti 


with values of coinciding with the characteristic numbers of Eq. (9). 

This equation was derived by M.S. Gurevich [2] for the pulse carrying the maximum 
part of the total energy in the given bandwidth. 

Thus, the pulse with the maximum concentrated spectrum (in the energy sense) during 
passage through an ideal filter provides minimum after-effect (also in the energy sense). 

Let us note one property of @(t), flowing from the physical sense of Eq. (9), which may 
be regarded as the Duhamel formula [compare with Eq. (2)]. In passing through the system of 
Eq. (1) the pulse@(t) does not change its shape in the interval (—T/2, T/2) but only decreases 
its amplitude by Ayiy times and is analytically extended from the interval (—T/2, T/2) over 
the entire time axis (—™, ©). 

If, instead of an ideal filter, we examine another narrow band system K(w) with the same 
cutoff frequency 2 then it is evident that the pulse Y(t) must be replaced by pulse w(t) so that 


¥ (0) K (@) = O(a), (10) 


where 


D(o)= | owe Mar; (11) 
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T/2 


(0) = \ p(j)e dt, (12) 


—T/2 


In other words, the pulse must be "matched" with the system so that there is obtained at 
its output a spectrum, corresponding to the function f(t), with the best energy concentration 
in time. 

This may be assumed as the basis of another approach to the problem. Let it be nec- 
essary to find F(w) — the spectrum (at the output of a system with finite bandwidth) equal to 
zero at | w| > Q — such that the energy of the after-effect 


09 T/2 


Wae= \ Pa— \ Pat (13) 


—co —T/2 


is minimum; that is, we arrive at the problem discussed in Ref. 2, the only difference being 
that the time and frequency axes are interchanged. Here 


Q 
i Oe : \ F(o) edo. (14) 


—2 


Instead of a determination of p(t), we seek its Fourier transform © (w) in the frequency 
interval (—Q, ), which obviously is a unique definition of ® (w) over the entire frequency 
axis and, consequently, of O(t) . 

Investigation of the stability of the functional 


oo Q T/2 a 
Rese Sia Ie Pear cowie toot 2 - 
\ |e \ F (@)e ‘do| dt — \ lz \ F (@)e do| dt (15) 
—oo —2 —T/2 —2 
or, which is the same, 
% BOs Sins, (@ — y) 
ai ‘ Kv 2 ' 
se \ | F (o) |?do — 55 \ \ Jory F (o) FY) dody (16) 
== Oe 


leads to the equation 


DB (o) = 2 | ET Diy dy. (17) 


—2Q 


As stated above, we have designated by ® (w) that one of the admissible functions F(w) in 
which the functional (15) assumes a stable value. Comparing Eqs. (9) and (17), we arrive at 
the equality 


o(g)=00(F), —F<r<y. (18) 


which is at least valid with the simple characteristic number XA »jyn; & may be determined by 
means of the Rayleigh theorem. 

In Ref. 2 M.S. Gurevich performs an approximate solution of Eq. (9), According to 
Ref. 2, the shape of Q(t) with small values of T2415 almost rectangular and with an in- 
crease in 2 T/2 there is observed a considerable decrease in ((t) toward the extremities of 


the interval (—T/2, T/2) and with = >s 7m on the plot of @(t) there appear at the 
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extremities of the interval (—T/2, T/2) oscillations about the zero level. For values of 
QT/2 from the interval (7, = 1), O(t)is close to the function cos > tb 


In summation, it may be said that the pulse G(t) is a maximally frequency-compressed 
pulse. In passage through an ideal filter it is converted into a maximally time-compressed 
function repeating the shape of @(t) in the interval of the initial duration. The spectrum of the 
pulse, © (w), when Ayin is simple, is associated with the pulse by Eq. (18), that is, is its 
dual in the sense of the Fourier transform. A system with cutoff frequency © yields at the 
output a maximally time-concentrated function if the product of the spectrum of the pulse 
and the system transfer function coincides in the passband with © (w). 

Hence it is evident that compensation must reduce to selection of such a pulse with the 
given linear system (or of a compensating filter with a given code signal) that the spectrum 
of the signal at the output coincides with ® (Ww) within the passband of the system. 

Let us mention again that we have discussed only idealized linear systems with infinite 
attenuation in the stopband. Energy was chosen as the criterion in evaluating transient in- 
tensity. 
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CALCULATION OF CHANNEL CAPACITY 
IN TWO-PATH COMMUNICATION 


Lin Hai-Ts'uan 


The article discusses the capacity of a two-path communications channel with randomly 
varying parameters (propagation amplitude and propagation time) in the presence of Gaussian 
additive noise. A formula is derived for the capacity of the channel in the case of normally 
distributed random amplitudes and a uniform noise spectrum. 


INTRODUCTION 


Much attention is currently being devoted to the investigation of the theory of communica- 


tions channels with randomly varying parameters. These investigations include Refs. 1 and 2 
and others. 
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Reference 2 investigates the capacity of a two-path channel (frequently encountered in 
practice) and derives a simple formula for determining the capacity in the case of a uniform 
noise spectrum. Instead of a real two-path communications channel, the present paper ex- 
amines an idealized model of a channel in which the amplitudes of the two waves are con- 
sidered to be nonfluctuating. 


CALCULATION OF CAPACITY 


Let € (t) be the transmitted signal at the input of a channel having finite average power 
Pg and finite bandwidth (O, W), within the limits of which the propagation conditions for any 
of the spectral components of &(t) may be considered identical. In addition, it is considered 
that two-path propagation of radio waves occurs. 
= Under these conditions the signal 7 (t) at the channel output may be represented as 
ollows: 


9 (2) = a (4) & [tf — 1 (¢)] + ae (t) § [¢ — Te (¢)] + (2). (1) 


Here a4(t), ag(t), T(t) and To(t) are random functions, independent of &(t) and € (t), 
which describe amplitude and arrival time fluctuations of the corresponding waves; € (t) is 
additive noise, independent of € (t). 

As shown in Ref. 2 the calculation of the capacity of the two-path channel described by 
Eq. (1) in the general case is a complex nonlinear problem, 

The physical nature of the propagation of medium and short waves permits substantial 
simplification of Eq. (1). It may be replaced by an idealized model of the channel 


m (f) = a8 (6 — %) + as (t — te) + 6), (2) 


where a, and ag are fixed, nonrandom numbers equal to the root mean square values of the 
processes a(t) and ag(t); T, and Tj are random quantities whose distributions coincide with 
the one dimensional distributions of the processes T,(t) and To(t). In Ref. 2 there is derived 
a formula for determining the capacity of just such a channel. Below we examine a model 
which may be described by the relationship 


 (t) = 28 (¢ — 11) + yS(t—Te) + OM). (3) 
6 


Here x and y are random quantities whose distributions coincide with the one dimensional 
distributions of the processes a, (t) and ay(t). Such a model provides a closer approximation 
to the real communications channel. 

Let us assume that the transmitted signal & (t) is a random stationary process. Let the 
additive noise C(t) be a Gaussian random stationary process with zero mean. The mean 


noise power will be designated by Pn. 
The channel capacity of the model of Eq. (3) may be represented by the expression [3] 


Cc” = max MI(E, 1/2, y, M1; Te); (4) 
2 (t)} 


where I(é, 1/x, y, T 1, Tg) is the information rate; averaging is performed over allexemya 
Ty, Tg. If 7 and To are considered to be independent of x andy, then 


Cie ak \\P (x, y) dady\\ P (t1, Te) 1 (E, 1/2, Y, Tr, Te) dty dt. (5) 


As shown in Ref. 2 in the presence of additive noise with a uniform spectrum the 
capacity of the model of Eq. (2) 
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Co AX MT (é, Y/ tas T2) = 
{& (¢)} 


P Tye - Py} 
=F logs E + (a2 + a?) P. A ae (aj) Be + 2 (a2 + a?) =| ’ (6) 
Hence in the case of the model of Equation (3) 
Pe 
Gi a= \\ Pe, y)F log 4 + (a? + y?) Pa + 
; 'z 2 ' P, 
+V De eae Pe + 2 (2? + y") |acay = 
os mr \\? (x, y) O (2, y, Pe, Pn) dx dy, (7) 


where 


1 Als 
© (ey, Pe, Pa) ng | Lb G+ vA) pS + 


Bs ch Wt 2 P, 
sae 1+ — P42 +H) \e 


Let us assume that x and y are independent random quantities with normal distribution. 
In this case 


(x—d,)? (y—ay)* 
ae MS , SA Ger { a 207 P 
Coates \ \ 1g “1 ue 8 Oe yee Pere (8) 
In2 V 203, V 2150 
09 —0O 


where aj, ag, O and Oo are, respectively, means and standard deviations of the random 


quantities x andy. 
If the dispersion of the random amplitudes is small, then, expanding ® (x, y, Pa Pn) 


in a Taylor series in terms of x = aj, Y= ag, we obtain 
® (z, Yy; Ie Pn) = ® (a, a, Pe; Py) Sie (w@—a,) © (a, a2, Bes Pn) a 
a 4 2 2 
ae ts), ® (a,, ds, Po, Pn) + pea)? Ss Ou, dg, Pe, Pn) + 


o 
2(z a) (y a2) 355, (a, do, 1 Pn) + 
a2 
+ (y — 4)? Fa D(a, aa, Ie, Pa) +.... (9) 


Substituting Eq. (9) in Eq. (8) and integrating, we obtain 


W- 5 


rs F , Sy @ 
CY = ae (Q1, dz, Po, Pn) + a aa WDC. Wigs 1D, Be) — 
22 @ | 
ae re WGA, Oh, 225, Pye o. 7 : (10) 


The resulting series is an approximate expression for the capacity of the two-path 
channel with randomly varying parameters in the presence of additive noise with a uniform 
spectrum. Equation (10) is applicable in the medium- and short-wave ranges. 

Let us examine a few particular cases. 

1. In the case where fluctuation is absent (i.e. , G =O = 0) the expression for C'' is 
transformed into the expression for C' obtained in Ref. 2. 

2. When the ratio ErADe is small 


ake eee (11) 


This means that for small signal-noise power ratios the capacity of a two-path channel has 
the same value as a one-path channel with a transfer constant (for power) a, + ae + of +: of ‘ 


3. If the ratio P,/P,, is large, 


P 
“ 5 Cc . D 3 
F log @ p— With ai >a}, 


GC! ew n (12) 


12 
) py tir) ig * 3 3 
vi log as 7, wih ay << ay. 


The result coincides with that obtained in Ref. 2 for this case. This is evident since, 
with a large ratio PA P,, and small dispersions 0, and 09, the amplitude fluctuations are 
negligible. 

4, In the case where the second wave is absent (i.e., ag = 0» = 0) and the output signal 
is represented in the form 7 (t) = x€ (t —T,) + (t), we have 


Fo GP. 


CO" ~ 


In2 ee (ai 7 oi). (13) 


In conclusion I wish to thank V.I. Siforov and B.S. Tsybakov for their constant as- 
sistance in this work. 
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CAPACITY OF ONE-PATH AND TWO-PATH 
COMMUNICATIONS CHANNELS 
WITH POLARIZATION FADING 


Lin Hai-Ts'uan 


The report discusses the capacity of one-path and two-path communications channels 
with polarization fading. Formulas are derived for determining the capacity of such channels 
in the presence of Gaussian noise with uniform spectral density. It is shown that the decrease 
in capacity of a one-path channel due to polarization fading does not exceed 14.5%. 


INTRODUCTION 


The capacity of two-path communications channels with interference fading was dis- 
cussed in Refs. 1 and 2. However, in short-wave radio communications, in addition to in- 
terference fading, there exists polarization fading of the received signal, which also affects 
the capacity of short-wave communications channels. Hence it is of interest to calculate the 
capacity of communications channels with polarization fading. 


1. ONE-PATH COMMUNICATIONS CHANNEL 


Let us first examine a one-path communications channel. As is known, under the in- 
fluence of the earth's magnetic field a plane-polarized wave reaching the ionosphere is 
divided into two elliptically polarized waves. Along the propagation path of the radio waves, 
the direction of the electric field vector of each wave changes continually due to fluctuation 
variations in the concentration of electrons. The change in direction of the plane of polariza- 
tion will lead to fading of the signals. 

Let it be assumed that reception is achieved with a horizontal dipole and that the direc- 
tion of the electric field vector of the incoming wave relative to the dipole assumes various 
equiprobable values. Mathematically this is described by the following: 


1 
W (8) = { 2m] a el et? 


4 
0 for 6<0, 6>2zn, (1) 


0 A where @ is the angle between direction of the 
electric field vector of the incoming wave and 
Fig. 1. Relation between direc- the dipole Fig. (1): W(@) is the probability 
tions of the electric field vector E density of the angle 6. 
of the incoming wave and the hori- Let &t) be the transmitted signal with mean 
zontal dipole A. power P, and 7 (t), the received signal at the 
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channel output. Then 
1 (t) = 2 (¢— 1) + 6(4), (2) 


where C(t) is the additive noise, the mean power of which is designated as P,; x is a random 
quantity characterizing the component of the electric field vector in the dipole direction. 

If aj represents the transfer constant (for electric field strength), then in the given 
case x = a, cos 9. As is shown, if the angle 9 has a uniform distribution, then the probability 
density of the random quantity x has the form 


“ 


For a channel operating as described in Eq. (2) the capacity 


P (a) = for |#| <1, (3) 


n 


; P 
Cnea\ Pie) Plo (1 pe) de. 
\ (A VOg | int Hap | da (4) 
Substituting Eq. (3) in Eq. (4), we obtain 


, ; P 
Cy = \ a Flog( 4 ++ a? ) ae, 


nV a? — 2" a0 


(5) 


, 
—a, 


After integrating, we have 


5 : 4 aP., ‘ Pp 
Commi 108% (1+ IP +/ 1+a5* 


ale © 


This equation permits the calculation of the capacity of a one-path channel with polariza- 
tion fading in the case of Gaussian noise with a uniform spectrum. It follows from this that 
the capacity depends on the bandwidth occupied by the signal and the ratio of the received- 
signal power to the noise power at Po/2Py at the channel output. 

Let us examine a few particular cases. 

1. When the signal power is small in comparison with the noise power the capacity of 
the channel discussed is 


— 
“sl 
Ww 


4 aP. 
CieePlog(t -4 a ; 
n 


Hence it is seen that the capacity of the channel with polarization-fading coincides,in this 
case with the capacity of a one-path channel with a transfer constant (for power) ay Ys 
2. When the ratio P,/P,, is large 


ie 


CroxF log aj ak (8) 


It follows from this equation that the capacity of the channel discussed coincides with the 
corresponding quantity for a one-path channel having a power transfer constant ae, 

Now let us compare the capacity Cj, in the presence of polarization fading with the ca- 
pacity C; of a one-path channel in the absence of polarization fading. The average power of 
the received signal at the channel output (designated as P)) is assumed to be identical in 
both cases. 

From Shannon's formula 


Po 
C,=F log (1 + Pr ; (9) 


where P! = Past (see Ref. 1). In the given case the quantity x2, equal to at/2, is easily 
found from its probability density. Hence 


a2P ! ‘ 
Ck log (1 7 (10) bn 
Figure 2 shows the depenaene: of ratio a9 
Ci/Cy on Pe Pat where P}, = a7P_/2. It is * aa 1/0, 
easily seen that with a change in P}/P, a a ain WS at 


from 0 to ©, the ratio C}/C, decreases 

from 1 to 0.855 and then rises and approaches Fig. 2. Capacity Cw, C, asa 
unity. Hence it follows that the capacity of a 
one-path channel decreases due to polariza- 
tion fading by not more than 14.5% at any (Po/ Pa) 
signal-noise ratio. 


function of signal-noise ratio 


2. TWO-PATH COMMUNICATIONS CHANNEL 


The preceding section discussed the capacity of a one-path channel with polarization 
fading. If, due to single and double reflection from the ionosphere, there arrive at the 
receiving point two waves in each of which polarization fading occurs, then there is formed 
a two-path channel with polarization fading. Let us investigate the capacity of such a channel. 

As shown in Ref. 2, the capacity of a two-path channel with randomly varying amplitude 
and phase in the presence of white noise may be calculated by means of the equation 


” 


Cn= \\ P(x, y)\C (a, y)\dz dy, 


where 


, z P 
GC (x y) = Flog +1 + (a? + y?) xe a 


P? ] 
a V/ 1+ yy Bee ach ag) _ ie 


We shall assume that the quantities x and y are independent. Substituting the expressions 
for their probability densities in the expression for Cy we obtain 


4 te@— vps 42ers we bee 
y pe Se Ae +p di dy. (11) 


The equation permits the calculation of the capacity of a two-path channel in which 
magneto-ionic division occurs. 

Let us examine a few particular cases of Eq. (11). 

1. When the ratio P,/P, is small 


Te 


Cc 


Po (a + a). (12) 


’ F 
Cie 3 


The capacity of the channel discussed coincides with the corresponding quantity for the 


channel investigated by Shannon with a power transfer constant of _1 (a2 ae} 
2nd Qi 


2. In the same case where a1/a5 => 1 and pep i, 
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* iy 4 OR, (ime tiga Pay, Nhe Wes a5 
oa (1+ 2Ph +// maa Ales Py ae (13) 


In the limiting case with a+ 0, the quantity C}! goes over to Cr (Eq. (6)), which was 
to be expected. 


3. When the ratio P,/P,, is large 


P. 
Flogat-=— for a2 > a?, 
*) © n 2 
Chom (14) 


P 
F loga?—* for a@< a?. 
2 Py if Si 


From Eq. (14) and the corresponding equations of Refs. 1 and 2 it follows for this case 
that with a large signal-noise ratio the capacity does not depend on whether or not the signal 
amplitudes fluctuate. 


In conclusion I wish to express my thanks to V.I. Siforov and B.S. Tsybakov for their 
constant assistance and attention in this work. 
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THE PERTURBATION METHOD AND ITS RELATION 
TO THE EXACT DIFFRACTION METHOD 
IN SOLVING PROBLEMS OF INCOHERENT SCATTERING 


A.V. Shabel'nikov 


The exact solution of the problem of the diffraction of a plane electromagnetic wave 
at a dielectric sphere is expanded into a power series in the small parameter w = Ank,a. 
It is established that the solution obtained by the perturbation method is identical to the 
first term of the expansion. The limits of applicability of the perturbation method are 
evaluated. The perturbation method is used in solving the problem of the scattering of a 
plane electromagnetic wave at a dielectric sphere with constant and variable €¢, The solution 
with variable € under certain conditions, imposed in the form of a function € (r), is identical 
with the solution obtained in the scattering of radio waves at statistical discontinuities of 
the troposphere. 


INTRODUCTION 


In solving the problem of incoherent scattering of radio waves at discontinuities of the 
troposphere the perturbation method [1-4] is usually used. It is distinguished by simplicity 
and permits the calculation of the scatter field from volumes of any shape with any dependence 
of the dielectric constant on position. However this method is approximate and hence may be 
used only under definite conditions. 

It is the purpose of the present paper to determine the relation between the perturbation 
method and the exact diffraction method and to establish the limits of applicability of the 
perturbation method. For this purpose the exact solution of the problem of diffraction of a 
plane electromagnetic wave at a sphere is expanded into a series in terms of powers of 
a = Ankoa and it is established that the solution obtained by the perturbation method is 
identical with the first term of the expansion. 

In the conclusion the perturbation method is used in solving the problem of the scat- 
tering of a plane electromagnetic wave at a dielectric sphere with constant and variable e. 
The solution with constant € is identical with the exact solution of the diffraction of a plane 
electromagnetic wave at a sphere obtained by Jobst [5] on the same assumptions under 
which the perturbation method is applicable. 

The solution with variable ¢€ (r), where r is the distance from the coordinate origin to 
the reradiating point of the sphere with an appropriate selection of the function € (r) is 
identical with the solution obtained in the scattering of radio waves at statistical disconti- 
nuities of the troposphere by the use of the correlation function R(r). 


1. RELATION OF THE PERTURBATION METHOD TO THE EXACT DIFFRACTION METHOD 


Let us examine the propagation of a plane electromagnetic wave in a medium in which 
the dielectric constant is a function of position. In this case the equation satisfied by the 
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electric vector E will be [2] 


AE + heel = + *F ix,, (1) 


where 


e=e+ SAe; S<1; k=, &) = const; (2) 


J is the current density in the antenna; A€ in the general case is a function of position. 
Substituting Eq. (2) into Eq. (1), we obtain the discontinuous wave equation 


AE + het = — AeShh + eerily, (3) 


the solution of which may be written in the following form: 


= eihaR 


ree faye 2 
= Br) 5 he ae GF av @ 


» 
r 


E( 
where rg is the vector coordinate of the point of observation; 138 (Z5) is the field created 
by the transmitting antenna at the point of observation in the presence of a homogenous at- 
mosphere; V is the reradiating volume; R is the distance between the point of observation 
and a variable point of reradiation in the volume V (see Figure); kp is the wave number in 
a homogeneous medium. 

In solving the integral equation let us use the method of successive approximations [6], 
for which it is necessary to find the solution of Eq. (4) in the form of a series in terms of 
positive integral powers of S: 


> 


E (r,) = Ey (74) + SE, (ro) + S*Eg (74) +... + S?Bp (%) b... - (5) 


Substituting the series of Eq. (5) into Eq. (4), integrating term by term and equating coef- 
ficients of like powers of S on both sides of the resulting equality, we obtain the formulas 


for the sequential determination of Ep(¥Q) : é 
Eq (7) = Eo (70); (6) 
= 4 SOIES BIRD 
Bi.) = a \ Aek3E, (7) -— dV, (7) 
V 
> => 4 eee etkaR 
Bx (ro) = yey \ AekiBs () Q- ai, (8) 
Vv 
poe [ Ackie: (oh eS (9) 
Py Ames ) Brash R g 


where 


R= Vr ae ro 277 COs &} 


cos £ = cos © cos O, + sin O sin O,-cos (p— %); 7, 0, @ and 79, 9%, @o are the coordinates of the 
points of reradiation and the point of observation in the spherical coordinate system. 

We will note that Eq. (7) is the classical expression for the scatter field, which is us- 
ually derived in all theories of scattering in using the perturbation method. In order to es- 
tablish the relation between the perturbation method and the exact diffraction method it is 
necessary to precisely evaluate the integral of Eq. (7), wherein in this case A€ = const. 


iNehee 


Usually in scatter theory in evaluating integrals of this form the quantity R is ex- 
panded into a series, thereby taking advantage of the fact that the calculation of the field 
is being made in the wave zone. Hence in evaluating the perturbation method we must know 
the value of the integral in the vicinity of the scatter volume, that is, at points where rp 


is of the order of or even less than r. 
In this case the expansion of R into a series has no meaning and it is hence necessary 


to find a method permitting precise evaluation of an integral of the form of Eq. (7). 
Thus, it is necessary to evaluate the integral 


> Agk? > eikaR 
Ey =. = io ethir cos © oe -dV, (10) 
Vv 


where Eo is the amplitude of the plane wave. For this purpose let us expand the integrand 
into a series in terms of Bessel functions and Legendre polynomials [7] 


oO 


> Aeke = (7S 
B, = F,\ {2 i (2n + 1) Wa (Kar) Pn (cos @)} x 
V ‘n=0 
x {ike D} (2m +1) Fr (har) AL? (hyr.) Pn (cos 8)} AV, (tt) 


where 


ahr) = Vi cd a (har); 


PL nafe ae 
2 aaa 


HD (karg)= Vo se HO 4 (ker), 


while the function P_(cos €), in accordance with the addition theorem [8], may be written 
in the form of a series 


P,,(cos &) = P, (cos @) Pn (cos ®,) + 
- (n—m)! pm m 
+2 2 Teme Pr (cos ®) Pr’ (cos @,) cos m (@ — go). (12) 


Substituting Eq. (12) into Eq. (11) and performing the successive integrations for ~, © and 
r, we obtain the final expression for the radial component of the scatter field: 


wAnkoa? 


Ey = EO By D) i” (Qn + 1) h® (hear) P} (cos @y) cos @, X 


2ro 
(Higa) ses (haa) , (13) 


n=0 


f 2 « ‘Saini 
x te n (kya) — J 


are: 
2 


where Eo is the amplitude of the wave incident at the dielectric volume V; a is the radius of 
the dielectric sphere; ky is the wave number in the medium. 

Let us now examine the exact solution of the diffraction of a plane electromagnetic wave 
at a dielectric sphere of radius a. Let us represent the radial component of the scatter 
field outside the sphere as follows [10]: 


eve} 
Haas eu . 
Ey = — Pa > int (2n + 1) FAY (kor) Ph (cos Oy) cos Do, (14) 


n=0 


where Wis. (y) [e¥, (x))/ = o i, (x) (yt, (y)]’ 


HD W) eV DY — SV, CO LA DY | i 
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NOM 0% Nees 


> p= kya; 


the primes indicate differentiation with respect to the argument; x — y = a@ = Ankga; 


An = Ny — Ng; Ny is the refractive index of the medium; n, is the refractive indexe of the 
sphere. 


Let us designate 


[2 (x)]’ 
An (@) = ee 


whence Eq. (15) is rewritten as follows: 


(yy) — A, @) yt, (y) 


bn = By (@) = - 5 
(yh (y)\’ — A, (@) yn) (y) 


(16) 
Since x = y+ qa, the function Barly + @) may be expanded into a Taylor series: 
Bu (y +o) =Bn(y) + 77 Bn(y) +... 4-2 = BP (y + 72), (17) 


where y is a number lying between 0 and 1. 
Calculations show that 


r Bn (y) = 9, : (18) 
Bn (y) = — ify? [Pn (y) Pr (y)— (Pa (y))?) —2 (Wn (YP — y'P rn (y) En (y)}- 


After appropriate transpositions, Eq. (14) in the first approximation takes the form 


By = “RE Ey S i" (2n + 1) he? (karo) Pr (cos ®,) cos Po X 


x (rng 2 (hat) — J a (hd ol (hea)} . (19) 


a 
n—=- 1 plore 


Comparison of Eqs. (13) and (19) reveals their complete identity, which leads to the 
following important conclusion. The expression for the scatter field usually used in the 
theory of noncoherent scattering represents the first term of the expansion of the exact so- 
lution of the diffraction problem into a series of powers of & = Ankoa. Evaluation of the 
integrals (8), (9) etc. by the above method has permitted establishing the following identi- 
ties: 


ies = @ y\ pnt (2n 4-4) S Br (y) WNP (Far) Pn(cos Og) COS Gy, (20 
oT 
ag n=0 
Ep = — Fone BS) i?* (2n + ire <3 (y) RD (kar) Pn(cOS Oy) COS Gy: (21) 


n=0 


Equations (19), (20) and (21) show that the series (5) is an exact solution of the problem 
posed, wherein the rate of convergence of this series is determined by the rate of con- 
vergence of series (17). Series (17) permits the determination of the limits of applicability 
of the perturbation method. For this purpose it is necessary to satisfy the following in- 
equalities: 


oP ay + 70) <a (y) (22) 
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” 


or 
a By (y+ 12) 
> a <i. 23 
PiegieG) = co) 

Evaluation of the quantity 41 = [Bn (y+712)/Bn(y)| shows that it does not exceed unity for any 

values of the argument. This indicates that inequality (23) may be rewritten as follows: 


Inequality (24) defines the limits of applicability of the perturbation method. Chernov [9] 
also arrives at such an inequality in applying the perturbation method to an examination of 
the scattering of sound waves at fluctuations in the density and velocity of the medium. 

Let us determine the number of terms of series (5) which must be taken in order to 
obtain the solution of integral equation (4) with an error not exceeding a given value of Uy. 
It will be determined from the following equality: 


p 
sr Be (y+ 7) = Bre (25) 


The procedure for establishing the limits of applicability of the perturbation method may 
easily be extended to dielectric volumes of any shape. For this purpose let us examine an 
arbitrary volume V' with a dielectric constant differing by 4¢« from the surrounding me- 
dium. Let us inscribe around this volume a sphere with radius a, where 2a is the distance 
between the two outermost points of the volume V'. The sphere volume V is 


V=aV'+V", (26) 


where V" is the volume required to bring V' up to the volume of the sphere. From Eq. (26) 
it follows that 

sae SS eik2R ae eS 
| Aeng#, (7) G- aV < EL), (27) 
&, 
that is, condition (24) is also applicable for the volume V'. 


1 
4m& 


2. APPLICATION OF THE PERTURBATION METHOD IN SOLVING THE PROBLEM OF 
SCATTERING OF A PLANE ELECTROMAGNETIC WAVE AT A DIELECTRIC 
SPHERE WITH CONSTANT AND VARIABLE e€ 


A. Let us demonstrate the advantage of the perturbation method over the exact diffrac- 
tion method in the problem of scattering of a plane electromagnetic wave by a dielectric 
sphere with constant «. The coordinate systems and geometric relations used in this case 
are shown in the figure. 


Poly Yo 


200; Yo, 


Geometric relations in the scattering of radio waves at a dielectric sphere. 
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It was shown in Section 1 that in the first approximation the field scattered by the 
sphere may be represented as follows: 


(28) 


where E(t) is the incident plane wave propagated in the positive z direction: 


=> 


oe Oe: 
sO (73) ee eres 


(29) 
Substituting Eq. (29) into Eq. (28), the field scattered by the sphere may be represented as 
follows: 


E, (r)) = Aek®E,L, (30) 
where 
4 ( oiklz+R) 
ECR R 
Vv 


(31) 
Let us evaluate the integral of L in the spherical system of coordinates. We will con- 
sider that the dimensions of the sphere are much smaller than the distance to the point of 
observation (a < ro), that is, we shall examine the field within the wave zone of the sphere. 
In this case the integral of Eq. (31) may be rewritten as 


iherye —ihy [Bat My (1—% | 

| ee ae ai e 3 i ae (2 aE 
4mro 

Vv 


(32) 
or in spherical coordinates 


0 
iker, ike2rsin — cos 
i = etelo it? 2 
4sro 


(33) 
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é 
where cos € is the cosine of the angle between the lines ( Sane 


(3) 
Denoting 2 sin = by b and considering that 
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sz)» Gp ande@, 9. 
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\ sin eae \ eibkar cos? dy — Ty (dkar), 


0 0 


jo (Ohar) = VJ, (Ober); 
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bkar) = iy wes sin bker, 


[4 
~~ 
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the integral of Eq. (33) may be finally written in the form 


etal 2 in Bhar 
L= "7, \ bky rdr, 


0 


where a is the radius of the sphere. 
After integrating in Eq. (34) we obtain 


pikaro sin bkoa ai if 
pe ( Z cos baa | ben? | ; ? : 


To bkoa 


The magnetic vector of the scatter field is related to the electric vector in the following 
manner: 


ii, = Aek3L Ves (rk, I, (36) 


where r is the unit vector directed from the coordinate origin to the point of observation. 
The power (per unit solid angle) scattered by the sphere in the direction of the ob- 
servation point is [10] 


Pg = Soy (37) 


where S;* is the projection of the vector S = > [E,Hy] in the direction of r*, while the 


total power scattered by the sphere in all directions, 


a ( ( Pasin Od0dg. (38) 


« 


0 0 


Let us define the scatter cross-section Q, as the ratio of total energy scattered per second 
to the power flux density of the incident wave, that is, 


ee 
tic S, ’ (39) 
where 9, = + Vc, £2. Substituting into Eq. (38) the expressions for E, andi, from Eqs. 


(30) and (36) and integrating with respect to 0 and@, we obtain for the scatter cross- 
section 


DO 1 1 sin 4k3a cos 4kg ie 
s = 2matks Ae? [ = = = 
& ace 4 (2k2a)* 8 (2k2a)4 4 (2kea)® | 8 (2ksa)4 4 | (40) 


Let us examine the two extreme cases: 


aks <1, (41a) 
Case (41a) corresponds to Rayleigh scattering, case (41b) corresponds to the scattering 
which occurs in long-distance tropospheric propagation of microwaves (a is of the order of 
tens of meters, A is of the order of tens of centimeters). 

With ak, < 1 the sine and cosine in square brackets in Eq. (40) may be expanded into 
a series in powers of 4k,a and limited to the first four terms of the expansion. After ap- 
propriate transposition fn (40) is rewritten as 


CG = ma’As? Be (4) ; 


v / 


(42) 
We have thus derived the well-known Rayleigh formula for the scattering of electromagnetic 
waves by small particles. 

With aky > 1 the first four terms in square brackets in Eq. (40) may be neglected and 
only the last term considered. In this case Eq. (40) takes the form 


(43) 
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Using for the scattering cross-section an expression in the form of a series obtained 
by the usual diffraction method, Jobst [5] obtained closed forms for the scattering cross- 
section: 


ee) a“ 
for the condition 
phe<p<I (45) 
and 
Qe =p ne Bet (S") (46) 
for the condition 
phe<l <p, (47) 


where 0 = 27a/X, wherein inequalities (45) and (47) are equivalent to condition (24). Com- 
parison of Eqs. (42) and (43) with Eqs. (44) and (46) shows their complete identity, which 
permits the following conclusion: being highly descriptive and mathematically simple, upon 
fulfillment of condition (24) the perturbation method provides solutions coinciding with the 
solutions obtained by the exact diffraction method. 

B. Let us examine the same problem as in section A, but let us now consider the di- 
electric constant of the sphere as dependent on r. On the basis of Eq. (28) the expression 
for the electric vector will have the following form: 


Dae es 
ee ro sl (48) 
where 
i — \ e(r)sin (bkyr) rdr. (49) 
0 


Substituting in Eq. (49) different functions of € (r), we obtain various dependences of E, on 
a, A and @®. For example, let us examine the relationship, widely used in scatter theory, 


r 


e(r) =Aee !, (50) 
where 


l<a (51) 


is a parameter which in scatter theory is considered as a correlation radius. Substituting 
Eq. (50) in Eq. (49) and_integrating, taking inequality (51) into account, we obtain the fol- 
lowing expressions for E, and Hj: 


2Aekel3 eiksro > 


ge E, (52) 
(1 + b?A31?)? ro 
— 2AEk213 thier > > 
H,=Ve, Cees Si 53 
es Dipmerrmersrterre (Att a 


where Vé,. is the mean characteristic conductivity of the medium. The power per unit 
solid angle scattered by the sphere in the direction of the observation point is the real part 


of ae [E, Hi] multiplied by ro. Then the scatter cross-section Q. is given by 
‘ omAerhdls ( ‘ | Uh 
ugh foie 42 (1+ 4h212)3 ) 
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For the special cases: 


(1) hal <1, 


Q. = 812A? Gor (55) 
(2) klS1, 
2 2m *% 
gee rey (08) 


From Eqs. (55) and (56) it follows that scattering at a dielectric sphere of radius a with 
variable € under condition (50) may be regarded as scattering at a sphere of radius / with 
constant ¢€. 


3. RELATION BETWEEN THE CORRELATION FUNCTION R(r) 
USED IN SCATTER THEORY AND THE FUNCTION € (r) 


It is of interest to determine that dependence of the dielectric constant of a sphere 
€ (vr) such that at the receiving point there is created a field equal to the field which would 
result from the scattering of radio waves at random discontinuities of the troposphere. 

The complex power density at the receiver, resulting from scatter within the volume V, 
will be [11] 


r 


digi kosiny)2 ,_ VAe? j 
shih =s V e2k5 | at fa |B (7) sin (bhar) rdr, (57) 
0 


where R(r) is the correlation function of the dielectric constant of air; V is the scatter vol- 
ume. The power density resulting from scattering of a plane electromagnetic wave at a di- 
electric sphere with variable € (r) is 


ke sin X¥)2 Ag? 
b2k2 


g(r) sin (bkgr) rdr) (58) 
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a ryan 1 — 2 
5 [LA] = 3 VeES| 
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Comparing Eqs. (57) and (58), we obtain the following equation for determining the relation 
between R(r) and € (r): 


bial ( R(r) sin (bkyr) rdr = ( e(r) sin (bkgr) rdr|’. (99) 
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Extracting the square roots and differentiating Eq. (59) with respect to r, we obtain the 
following relation between € (r) and the correlation function R(r): 


g(r) = R(r) y, . eee : (60) 
) 16x | R (r) sin (bar) rar 
0 


From Eq. (60) we may draw the following interesting conclusion. 
The scattering of radio waves at statistical discontinuities of the troposphere may be 
regarded as scattering at a sphere at the center of which the dielectric constant is V ae 


and decreases toward its extremities according to the law expressed in Eq. (60), while its 
dimensions correspond to the dimension / of the correlation variation. In the special case 
where € (r) is given by Eq. (50) the scattering of radio waves at statistical discontinuities 

of the troposphere may reduce to scattering at a sphere of radius / with fixed dielectric con- 
stant [11]. In other words, the problem of the scattering of radio waves at statistical dis- 
continuities of the troposphere may be reduced to the problem of the scattering of radio 
waves at a sphere whose dielectric constant is a regular function associated with the 
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correlation function by Eq. (60). 


In conclusion the author expresses his thanks to B.A. Vvedenskiy, A.V. Sokolov 


and N.A. Armand for review and discussion of the present work and their valuable com- 
ments and advice. 
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GRAPHICAL METHOD 
OF DESIGNING APLANATIC ANTENNAS 


N.G. Ponomarev 


The paper describes a simple method of designing reflector antennas and lens an- 


tennas satisfying the Abbe sine condition. This method is based on geometrical optics and 
is similar to the well-known method of designing mirrors from a diagram of special shape. 
The results of an experimental investigation of a two-reflector aplanatic antenna are given. 


In antenna engineering both reflector and lens antennas are used, wherein beam scan- 


ning is achieved by displacement of the primary source of radiation over a certain curve. 
As is known, undistorted beam scanning is possible upon satisfying the ''sine condition, "' 
which in optics provides undistorted images of objects of finite extent. By analogy with op- 
tical systems, antennas in which the sine condition is fulfilled are known as aplanatic an- 
tennas. 


Reflector and lens systems in which the sine condition is fulfilled have two reflecting 


(refracting surfaces, since, in addition to the sine condition, they must satisfy the ''co- 
phasality condition" ensuring the transformation of a beam of rays emanating from a point 
into a parallel beam. 
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Equations for the curves of the central section of an aplanatic reflector system 
are known [1], but they are complex and calculation with them is extremely laborious. 
Profiles of a double-refracting aplanatic lens may be plotted as the result of integration 
of the differential equation [2]. However, the analytic solution of this equation is unknown. 

Below we present a simple graphical method of designing both reflector and lens 
aplanatic systems. It is based on methods of geometrical optics and is similar to the 
method of plotting the curve of the central section of a mirror from given diagrams of a 
primary source and a mirror [3]. 

In order to fulfill the cophasality condition it is necessary that all the rays leaving the 
antenna be parallel to its axis. The sine condition is fulfilled if all the incident rays from 
the focus and their corresponding emergent rays intersect at a circle whose center lies at 
the focus of the system[2]. The radius of this circle is known as the focal length. 

In order to plot the profile of a reflector system satisfying both conditions let us draw 
a circle of radius f with its center at the coordinate origin (see Fig. 1). The coordinate 
origin is taken as the focus and line Fx as the axis of the system. We plot a system of in- 
cident rays from the focus and their corresponding emergent rays parallel to the axis and 
intersecting the incident rays at the circle. At the outermost ray and the corresponding 
emergent ray let us select arbitrary points A and B and join them by a straight line. These 
points are the extreme points of the profile of the reflector system*. Selecting points A and 
B, as shown in Fig. 1, and joining them by a straight line, we thereby define the path of the 
extreme ray. In order for the ray to have the shape of the broken line FABD the tangents to 
the surfaces of the reflectors at points A and B must be perpendicular to the bisectors of 
angles a and B (see Fig. 1). Plotting these tangents, we find the position of the elementary 
plane reflectors at the given points. Extending the tangents to the reflector profiles at points 
A and B to intersection at points A; and B, with the following incident rays from the focus 
and their corresponding emergent rays, and joining points A; and B, by a straight line, we 
obtain the path of the adjacent ray. Drawing at points A, and By straight-line segments, the 
normals to which coincide with the bisectors of angles @, and 81, we obtain the next ele- 
ments of the reflector profiles. Subsequent elements of the profile are obtained by continu- 
ing the plotting in the same manner. 

As the result of plotting we obtain the broken lines AA, A2A3...Ay and BB] B2B3...By 
approximating the profiles of the reflectors of a system satisfying the cophasality and sine 
conditions. It is evident that the smaller the 
angle A@, the more accurate the approxima- 
tion. The accuracy of the graphical plotting 
may be increased by other means, for ex- 
ample, by the method used in the calculation 
of lenses with a given amplitude-phase dis- 
tribution at the aperture [5]. Since the system 
is symmetrical relative to the focal axis, 
plotting may be limited to the profile of only 
one side of the axis. The second half is ob- 
tained by reflection relative to the axis. 

Figure 2 shows the plot of the profiles of 


F 
a two-reflector system for the case where 
Fig. 1. Plotting the profile of a two- point A lies outside the circle of radius f and 
reflector aplanatic antenna. point B lies within it. 


On the basis of graphical plotting we 
may obtain recursion formulas expressing the coordinates of points A, _ ; and Bier 
the two-reflector system these formulas have the form 


—(k—1)Ap4+ 8, 
tg [max —(# — 1) Ag] + etg | Sax ( va p+8, ) fi 
x = LA, ’ 
coun are 2 [Pmae —kAQ] + ctg| Pmax — (*— 1) Ag + ] 
5 Urmax / 2 


*The selection of initial points A and B is somewhat arbitrary. It is not difficult to see that 
both points may lie within a circle of radius f or outside it or that one of them (A or B) 
may lie within it and the other outside it. One or another choice of initial points determines 
the shape of the profile of the reflectors and their relative position. 
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Fig. 2. Plot of the profile of a two- Fig. 3. Plot of the profile of an aplanatic 
reflector aplanatic system (point A lens. 

is outside the circle of radius f and 

point B is within it). 


In deriving Eqs. (1) - (5) it was assumed that (),,,, is divided into identical angles AQ. 
Let us discuss the graphical plot of the refractive profiles of an aplanatic double- 
refracting lens. For this purpose, as in plotting a two-reflector system, let us plot a fam- 
ily of the rays incident from the focus and the corresponding rays emergent from the lens — 
intersecting at a circle of radius f (see Fig. 3). Joining points A and B by a straight line, 
let us assume that the outermost ray has the shape of broken line FABD. In addition, if the 
refractive index of the lens is n > 1, point A is chosen on a ray incident from the focus within 
the circle, while B is chosen on an emergent ray outside the circle*. Since the refractive 
index of the lens is considered known, by means of the law of refraction we can determine 
the position of the elementary refractive surfaces (at points A and B) ensuringthe given shape 
of a ray in the form of the broken line FABD. It is not difficult to determine that the angles of 
incidence B WN and Bp at points A and B are, respectively, 


we —nsind, 6 

By anche); +ncosa, ’ (6) 
n sina 
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where @, is the large angle between lines FA and AB and @ is the small angle between lines 


AB and BD. 
The angles of incidence Bq and Bp, and consequently also the normal to the elementary 


refractive surfaces at points A and B, can be plotted graphically. Let us discuss the plot of 


*It can be shown by graphical plotting that for an aplanatic dielectric lens (n > 1) point A 
must lie within the circle of radius f, while B must lie outside it. For a lens withn< 1 


points A and B must lie within this circle. 
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the normal to the refractive surface at point A. Let us plot a unit circle (see Fig. 4) 

with its center at point A and on line AB let us mark off segment AN', the length of which 
is equal to n. From the point N' we drop a 
perpendicular to the line FA and draw a line 
from N' to M', which is obtained as a result 
of the intersection of the unit circle with the 
extension of the line FA. As is seen from the 
figure and from Eq. (6), the angle N'M'P is 
equal to the angle of incidence 8 q at point A, 
while line M'N' is parallel to the normal to 
the profile of the refractive surface. Thus, 
the position of the elementary refractive sur- 
face at point A has been determined. Plotting 
of the normal at point B and at all subsequent 
. and B,, Bo, Bs, . = 


Fig. 4. Graphical plot of the normal points Aj, Ag, Ag, . : 
MN and the angle of incidence Ba at is performed in the same manner as at point A. 
point A. The entire procedure for plotting the profiles 


of the refractive surfaces of a lens is the 
same as in the case of a two-reflector system (see Fig. 3). 
As the result of the plotting we obtain the approximate shape of the profiles of the re- 
fractive surfaces of an aplanatic lens. The recursion formulas for the coordinates of the 
points of the lens profile have the form 


te 1Pmax — (* — 1) Aq] + ctg (a, —Ba,_,) 
tg (Omax = kAg) <r ctg (a), = Pa) 


, (8) 
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Disp = fsin (@ max — kAg), (11) 
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Plotting the profiles both of the two-reflector and lens antennas was performed by be- 
ginning with the outermost points. However, such plotting may be performed, by assuming 
that the initial points lie on the axis, and being given the thickness of the lens or the dis- 
tance between reflectors. Such plotting is preferred for lenses, for a lens of minimum 
thickness may be derived in this manner. 

Several types of aplanatic systems were developed by means of this method. A few pro- 
files of reflector and lens antennas are given in Fig. 5-8. As is seen from these figures, it 
may develop that one of the reflectors considerably shields the other (see Fig. 6). However 
in actual construction of a two-reflector system such shielding may be avoided by various 
means. For example, if for a given antenna the sine condition is satisfied in the horizontal 
plane, it is necessary to displace the reflectors in the vertical plane (see Fig. 9). Then 
there will exist no shielding between reflectors. 

The above method was used in designing a reflector antenna the experimental model of 
which is shown in Fig. 9. Scanning of the beam is achieved by the rotation of a parabolic 
horn, which for this purpose is mounted on a special bracket. The ratio of the focal length f 
to the aperture D was 0.8. Figure 10 shows experimental radiation patterns which were 
plotted by scanning of the pattern maximum by +50° from the antenna axis (in Fig. 10 the 
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angular values are referred to the width of the pattern in the central position). As seen 
from Fig. 10 the experimental model of the antenna (with {/D = 0.8) ensures undistorted 
scanning over a wide angular sector. The accuracy of the plotting of the reflector profiles 
was adequate for practical purposes. 

The author is indebted to L.S. Benenson for his interest in this work and his valuable 
advice. 


Fig. 5. Profiles of the reflectors of an aplanatic system with 
negligible shielding of the large reflector by the small. 


Fig. 6. Profiles of the reflectors of an aplanatic system in which, due to considerable 
shielding of one reflector by the other, it is convenient to displace the reflectors vertically. 
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Fig. 7. Profiles of an aplanatic lens. 


Fig. 8. Profiles of a "thin" aplanatic lens. 


Fig. 9. Experimental model of an antenna. 
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Fig. 10. Radiation pattern of two-reflector system. 
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SUPPRESSION OF CURRENTS EXCITED 
IN METAL SHIELDS BY DIFFRACTION ANTENNAS 


OF FINITE DIMENSIONS 


O.N. Tereshin, A. Ye. Sokolov 


The article discusses a method of suppressing electrical currents excited by a diffrac- 
tion antenna at a metal surface. Suppression is achieved by means of a purely reactive sur- 
face impedance. The method of suppression may be used for the decoupling of two diffraction 
antennas as well as for lowering the back radiation level and decreasing the effect of metal 
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shield dimensions on the radiation pattern of diffraction antennas. Experimental data 
show that the degree of suppression of surface currents is determined by the maximum 
rate of change of the purely reactive surface impedance which is attained in practical 
realization of the latter. 


INTRODUCTION 


Reference [1] discussed the possibility of creating a decoupling arrangement between 
two slot antennas of unlimited length. It was shown that in this two dimensional case the 
degree of decoupling is determined by the maximum attainable rate of change of the purely 
reactive surface impedance. Reference 2 pointed out that in the creation of decouplers 
for antennas of finite length the attainable degree of decoupling is limited by the level of 
the currents surrounding a decoupling structure built on the basis of the conclusions ob- 
tained from an examination of the two dimensional problem. , 

This paper presents the result of theoretical and experimental investigations of de- 
couplers for annular-slot antennas of finite dimensions. Results are given for an experi- 
mental check of the operation of decouplers constructed in accordance with the conclusions 
of Ref. 1 as well as of the effect of decouplers on the radiation patterns of slot antennas. 


1. STATEMENT OF THE PROBLEM 


Let us examine the problem of a decoupling arrangement for diffraction antennas of 
finite dimensions. 

An impedance surface is given in the plane z = 0. 

In the region S at the coordinate origin there is located a diffraction antenna (Fig. 1) 
in which the current distribution is symmetrical with 
respect to ~. In this case the boundary conditions at the 
impedance surface with z = 0 may be written in the form 


a j ma) 
ic 0) wit , (1) 


where E,., Eg» H,,., H,, are the tangential components of 
the electric and magnetic vectors of the field; Z(r) is the 
distribution function of the surface impedance. 

In obtaining a purely reactive surface impedance by 
means of a ribbed structure Eq. (1) is applicable in the 
case where the width of the channels is much smaller than 
the wavelength and the width of the ribs is much smaller 
than the channel width. 

Let us assume that the field of the primary sources 

Rigel given in region S has the free-space properties of TM 
waves. 

As shown in Ref. 1, the degree of limitation of cur- 
rent spread may be characterized by the degree of diminution of the tangential component of 
the magnetic field at the surface z = 0 in the radial direction. The solution is similar to the 
calculations performed in Ref. 1; hence we shall discuss only the principal stages. As be- 
fore, we are given the field in the form 


Ld i, BD) SOOO) (2) 
where 
Zr, 2) = Dee Vee 22, (3 Be (3) 


Here Z, and Zo are real functions of the arguments r and z. 
For the special case of TM waves from the primary source we obtain from Maxwell's 
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equations 


z i OH, 
Dr we Os ? (4) 
where € is the absolute dielectric constant of the medium. 
Substituting Eqs. (4) and (2) into Eq. (1), we obtain 
(ie a7,” We OLie a x 
BNE ae eg os mle DE 7: eMC Sa (9) 
The condition for a purely reactive surface impedance may be written in the form 
OZ, Nepal 
= = OD Wwaitlanea—t (0) (6) 


The tangential component Ho of the magnetic field must satisfy the homogeneous wave 
equation 


Ol HH. OL 
 ) (- *) ae aa = 2a fe 2H. (in (7) 
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Further, as in Ref. 1, in order to determine the class of functions satisfying the 
boundary conditions (1) and the wave equation (7) it suffices to consider the region adjacent 
to the plane z = 0 where the functions Z; and Z9 may be written as asymptotic expansions 
in powers of z in the form 


AC, 2)=2O+-2,(07)2- 2,02, (8) 
Lat, 2). = Lar) © Zp (rye -- Zp (r) 22. 


Substituting the asymptotic expansions (8) into the wave equation (7) and assuming z to be 
sufficiently small and then separately equating the real and imaginary parts to zero, we ar- 
rive at the equations 


2 LGD NG 4 ADO. ere é (ie 
Zy eh 25 V Ge :) ee) r ae a 2Z, (7) — (2 a = : @) 
BZ? i) om 5 ae: dZ) dZ\ 
Chee Ty dr 22; (Te Ae js = 0. (10) 


W " 
For the case of diffraction antennas it is necessary to assume in Eq. (8) that Z, and Zo 
are equal to zero. 

Using the equations derived in Ref. 3, it is easily shown that at those points of the 
surface z = 0 lying within the region S where the primary diffraction antenna is located it 
is necessary to choose the minus sign in Eq. (9) in order to fulfill the radiation conditions. 

With the limitations presented, the conditionforapurely reactive surface impedance is 


written in the form 


278 ily dL) 9 (%) Gon 


2 i 
dp 2a 


while the required distribution of surface impedance is found from the equation 


ay. (11) 
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If the coordinate of the observation point (r) in Eqs. (11) and (12) tends to infinity, then 
these expressions asymptotically become the expressions derived in Ref. 1 for the two 
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dimensional problem. 

Thus, the distributions of purely reactive surface impedance which have been de- 
rived for one and the same variations in the tangential component of the magnetic field 
about the surface z = 0 differ in the case of a symmetrical antenna of finite dimensions 
and the two dimensional case only at small distances from the antenna. 


3. EXPERIMENTAL RESULTS 


(a) Measurement of the distribution of the amplitude and phase of the tangential com- 
ponent of the magnetic field along the structure. Experimental investigation of the distribu- 


tion of the amplitude and phase of the surface currents excited by a slot antenna was per- 
formed by the usual method on the prototype shown in Fig. 2. The required purely reactive 
surface impedance was attained by means of a ribbed 
structure with period b = 0.156. and rib thickness 
Z ae A= 0.03 X. The length of the attenuating structure 
} was 1.56X. In the calculation of the structure, the 
function for the change in amplitude of the tangential 
y component of magnetic field H, along the y-axis was 


YYW 
MMe given in the form 


/ 


| Ha] = = Pu, (13) 
Fig. 2. 1) Slot antenna; 2) sup- 
porting structure; 3) measuring where B is the attenuation constant. Then the re- 
probe. quired variation in the purely reactive surface im- 


pedance, in accordance with Ref. 1, is defined by 


% t Anya el a 
Z{y) = = V tv C — B= k (14) 


Results of the experiment are shown in Fig. 3 and Table 1. Figure 3 gives the 
calculated (a) and experimental (b) distributions of amplitude and phase (for B = 0.475 k) 
of the field along the attenuating structure for various attenuation constants. The 
rectangles represent the channel depths of the ribbed structure for the case where 


B= 0.475 ke 
Table 1 lists the values of field attenuation 
,, atten, db caused by the entire suppressing structure. 


An investigation of the suppression of surface 
currents excited by a primary antenna in the form 
of an open-ended circular waveguide (Eg; mode) 
within a metal shield was performed on the model 
shown in Fig. 4. 

The required purely reactive surface impedance 
was attained by means of a ribbed-disc structure with 
period b = 0.12A and A=0.012. 

The variation of the amplitude of the tangential 
component of the magnetic field along the radial 
axis was given in the form 


—Br 


|#,|=—-. (15) 


Fig. 4. 1) Transmitting antenna; 2) suppres- 
sing structure; 3) receiving (slot) antenna. 
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Table 1 Table 2 


Measured 
attenuation, db 


The required variation of the surface impedance, according to Eq. (12), is deter- 
mined from 


: Yee ; : 
Lele) Sip ie meant eR eas ah ae (16) 


The measurement of the amplitude of the surface currents was performed on the basis of 
the degree of excitation of a receiving antenna (slot) placed a fixed distance away from the 
transmitting antenna. A suppressing structure of length 0.9 was introduced into the zone 
between the receiving and transmitting antennas. 

Results of the tests are given in Table 2 (only the additional attenuation attributable to 
the presence of the suppressing structure is listed). 

Let us evaluate the results. 

Examination of the results shows that in both cases there is adequate correspondence 
between the design and experimental values of field attenuation and the variation in amplitude 
of the field along the structure for small values of the attenuation constant 8. The deviation 
in results for large values of the attenuation constant may be explained as follows. With an 
increase in the attenuation constant there is also required a higher rate of change in the 
purely reactive surface impedance [see Eqs. (14), (16)], which condition cannot be realized 
by means of a ribbed structure with finite channel width. Actually, due to the finite channel 
width and wall thickness the required continuous distribution of purely reactive surface 
impedance is approximated in steps. Hence, in practice the variation in the purely reactive 
impedance is determined by the first 3-4 channels, after which the impedance ceases to in- 
crease due to inaccuracies in channel depth (see Fig. 3). With this in mind, it is convenient, 
in order to increase the attenuation, to periodically repeat over the entire structure the 
variation in surface impedance as produced by the first channels. 

If we extend over the entire structure the variation reflected in the first three channels 
for the case of a slot antenna and & = 0.475 k, then the attenuation increases considerably 
and reaches 33 db (Curve din Fig. 3). For the dise structure (repeating over the entire 
structure the variation reflected in the first two channels) the following results are ob- 
tained: with 8 = 0.48 k the attenuation rose to 25 db, with 8 = 0.64 k the attenuation rose to 
33 db. 

With given dimensions of the suppressing structure the degree of attenuation may be in- 
creased by using a ribbed structure with small period(b ~ 0.025, and less). In testing for 
the case of a disc with purely reactive surface impedance, such a structure was made in the 
form of dielectric rings of different height wrapped in brass foil. The structure was also 
made of such rings in close contact with one another. 

Measurements showed that with such a structure of size 0.1, the resulting decoupling 
was 10 db, which corresponds to B = 2.4 k. It is obvious that this is not the maximum value 
attainable. The use of ribbed structures with still smaller periods will result in considerably 
greater attenuation of the field. 

We note also that in periodic repetition of the segment with greatest rate of change in 
purely reactive surface impedance Eqs. (14) and (16) become inapplicable to the segment at 
a sudden reverse jump in surface impedance (transition to the next period) as soon as at 
this point aZ/dr —-o, andas phown in Ref. 1, Eqs. (14) and (16) are applicable with ac- 


curacy up to terms of the order Apes 


However, the curve for the change in amplitude of the field with periodic repetition of 
segments with sudden change in the purely reactive surface impedance (Fig. 3, Curve d) 
shows that the variation in amplitude of the field is quite close to that calculated for a 


195 


continuous change. This means that within the limits of each period of change in impedance 
Eqs. (14) and (16) remain applicable and at each point of jump in the surface impedance the 
coefficient of reflection is small. 

Thus, in the final analysis the rate of change in the purely reactive surface impedance 
can be determined. 

(b) Use of a suppressing structure in eliminating the influence of a metallic surface 
(beyond the limiting structure) on the radiation pattern of a slot antenna. For this purpose 
the radiation patterns of two models which differ in the presence or absence of a metal 
shield beyond the suppressing structure were plotted. The suppressing structure had a 
length of 0.1, and provided attenuation on the order of 10 db. 

Figure 5 shows the radiation patterns of a narrow slot antenna cut out of the center 
of a shield with length 6 ro (the schematic representation of each of the experimental models 
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is shown beneath its corresponding radiation pattern in each figure). Figure 6 shows the 
radiation patterns of the same slot antenna with an attenuating structure at the center of 

each arm, designed for two frequencies differing by 20% (a for frequency fg and b for fre- 
quency f; = 0.8fp). Figure 7 shows the radiation pattern for the same frequencies in the 

case where the outer shield portions are omitted. The patterns in Figure 8 are plotted 

from a model with a short shield (L = 0.63A 9) for the case of the omission of the attenuating 
structure at the end of the shield (wider pattern) and for the case of its use (narrower pat- 
tern). In comparing the radiation patterns for frequency fy in Figures 6 and 7, it is evident 
that elimination of the outer portions of the shield does not distort the radiation pattern. This 
indicates the absence of currents in these portions. 

The radiation pattern in this case is created both by the radiation from currents flowing 
along the shield and by the radiation from the suppressing structure (compare with Fig. 5), 
which is confirmed by the change in radiation pattern at a frequency f, differing from fg by 
20%. At frequency f, a slight distortion is observed upon eliminating the shield sections be- 
yond the suppressing structure. This is explained by a decrease in the cut-off properties of 
the ribbed structure because of the decrease in the drops in purely reactive surface impedance 
along the attenuating structure with a decrease in frequency. 

Installation of the suppressing structure at the edge of a small shield (Fig. 8) leads toa 
decrease in the back-radiation level. 
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CONCLUSIONS 


1. In practice, an adequate realization of the required purely reactive surface im- 
pedance results in an excellent agreement of the surface current attenuation with calcu- 
lated values. 

2. For the purpose of providing maximum attenuation, optimum results are obtained 
by the use of periodically repeated segments with sharply variable purely reactive surface 
impedance. 

3. A suppressing structure in the form of periodically repeated segments of ribbed 
structures with a small period permits considerably greater decoupling than is obtainable 
with known decoupling devices of the "cut-off slot" type [4] with identical dimensions 
(27 db in comparison with 9 db). An even greater degree of decoupling is provided by de- 
creasing the period of the ribbed structure. 

4. The use of a suppressing structure in the form of periodically repeated segments 
with sharply variable purely reactive surface impedance along its length permits elimi- 
nating the influence of a metallic surface beyond the limiting structure of the radiation pat- 
tern of a slot (diffraction) antenna. 

5. The use of limiting structures at the edge of the shield of a slot antenna reduces the 
back-radiation level. 

6. A ribbed structure with purely reactive surface impedance loses little of its lim- 
iting properties with a change of 20% in frequency. The changes in radiation patterns at 
different frequencies are apparently associated with the change in nature of the radiation 
of the limiting structure as well as with the change in current distribution in the segment 
from the antenna to the limiting structure. 

It is to be expected that by an appropriate choice of the radiation of the limiting struc- 
ture it will be possible to establish optimum variations in the purely reactive surface im- 
pedance for a limiting structure with a view to small distortions of the radiation pattern 
due to a shield of finite dimensions. 
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TRANSITION CALCULATIONS 
FOR A SYMMETRICAL MAGNETIC WAVE 
IN A CIRCULAR WAVEGUIDE 


B. Z. Katsenelenbaum, Z. A. Malina 


The paper discusses methods of numerical calculation of a symmetrical waveguide 
transition between two waveguides of different diameter for a symmetrical magnetic wave. 


INTRODUCTION 


In the transition of an H,, wave in a circular waveguide from a waveguide of one diam- 
eter to a waveguide of another diameter there arises reflection of the H,, wave and para- 
sitic waves H,, (i > 1) in both directions. In a number of cases the occurrence of parasitic 
waves substantially disturbs the operating conditions of the waveguide system. A decrease 
in the amplitude of the parasitic waves over a wide range of frequencies by increasing the 
length of the tapered transition requires the use of extremely long (more than | m) tapers. 
Below we present several methods for the calculation of tapered transitions with genera- 
trices of complex shape. The use of such tapered transitions permits obtaining small 
transformation with a relatively small length of transition. The calculation is based on the 
system of differential equations for wave amplitudes within a nonuniform waveguide as 
obtained in Ref. | and discussed in Refs. 2 and 3. 


The presence of the so-called critical section for a parasitic wave introduces a com- 
plication in the mathematical organization of the problem. The method developed in Ref. 4 
is used in the analysis of this case. The examples given in this paper illustrate the partic- 
ular characteristics of the proposed methods of calculation. 


1. SELECTION OF GENERATRIX 


In choosing the shape of a tapered section we shall proceed from the approximate anal- 
ogy (see Refs. 2 and 5) existing between the problem discussed and the problem of deter- 
mining the optimum variation in the characteristic impedance of a section of long line 
matching two lines with different characteristic impedances. 


The coefficient of reflection R from a section of long line with characteristic imped- 
ance W(z) is 


where L is the length of the section and € =z/L. 
On the other hand, if in the system of equations (12) of Ref. 1 or, what is the same, 
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Eks. (39)-(40) of Ref. 3 we neglect (except for the exponent) the difference of wave 
numbers h, and h, of the Hoi and H02 waves from the free-space wave number (in Sec- 
tion 2 the validity of this approximation is shown by an example), then the amplitude of 
the forward H o2 Wave will be 


1 
dina es 
D(s) = ch ca dé. - 


The relative power losses are equal to| D! 2. Here the parameter o is the phase excursion 
difference of both waves over the entire length of the transition. 


‘4 


G= V(r alts Be (3) 
0 


€ is anew variable, defined as 


= = (1 =a); (4) 


where) yj. — \ l,dz; C is a constant equal to 


0} 


DS 
2uyu 


(en L (5) 


2 
in 


is the p-th root of the equation J 1 (4) = 9. With the same degree of accuracy as 
exists in Eq. (2), the relation between € and z may be written in the form 


Doe aee 

z(t) = ——\a" dk, (6) 
av 0 Se 

y § 
where 

1 

a= \@ (SB) ches (7) 
0 


a — is the transition radius considered as a function of variable €. 

Comparing Eqs. (1) and (2), we see that in Eq. (2) D and a are related in the same 
manner as R and Win (1). The parameter 2kL is replaced in Eq. (2) by the parameter Oo, 
which is defined by 


ho (8) 


ao decreases as frequency increases. Further, the amplitudes of forward parasitic waves 
created at the transition will increase, since the difference in phase constants of the var- 
ious modes decreases, as does the mutual compensation of parasitic waves which arise 
in the various waveguide sections. We will note that in Eq. (1) the wavelength differs 
from that in Eq. (2), for the expression for the amplitude of the reflected wave includes 
the sum of the phases of both waves rather than their difference. The functions a (€) will 
be chosen in the same manner as were W(é) in the theory of long lines. Given a(&), from 
Eqs. (6) and (7) we can find z (€), that is, we can determine the shape of the generatrix in 
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parametric form. The amplitude of the parasitic wave is determined from Eqs. (8) and 
(2). The functions given in Ref. 6 will be used in choosing a(&). Let it be noted that deter- 
minations of the true optimum [7] for the integral of Eq. (2) are of no significance in our 
case, for Eq. (2) itself, for the problem under discussion, is approximate. 

The form of the function a(€) depends on the maximum permissible amplitude of the 
parasitic wave. Assume, for example, that |D| must be of the order of (3-5) - 1072 (the 
amplitude of the Hg; wave is assumed to be unity). In this case we may use the function 
W(é) defined by Eq. (8) in Ref. 6. Let us assume further that with € = 0 and = 1 the 
waveguide radius assumes the assigned values of a (0) and a(L). Then 


a(t) =a(0) a q (< -2:838 sin oxk) " 


where q = a(L)/a(0). At the ends of the generatrix defined by this function there are 
breaks equal to 


da Gey 
rT. = 0,364 Ing L 


aos (10) 


At the wide end the break is smaller than at the narrow end. 
For the variation defined by Eq. (9) the function |D(o )|has the form 


| sin 2 
2 An? — 0,364 o2 
| Di(s) |= | 6 Inge (11) 


SG 
2 


The plot of |D@)|/|Inq|is shown in Fig. 1. With o>cmin =37, /D@)| < 0,048{Inq| . 
Knowing q and denoting the minimum wavelength within the band as Amin, from Eqs. (7) 
and (8) we may determine the minimum length of the transition. Let us note that a tran- 
sition providing sufficiently low loss for the Ho2 wave will provide [within the limits of Eq. 
(1)] even lower loss for the Ho3 wave, etc. 


WDi\tng\ 


a(z), MM 


re) pari 
O° On bn 4n on 6n 2 ne) eT) 
Fig. 1 Fig. 2 


Let us consider an example. If a(0) =9 mm, a(L)= 30 mm and Amin=6.6 mm, then 
the transition length will be 193.7 mm. At all X>2 min [within the limits of applicability 
of Eq. (2)], we will find| Dia)|< 0.06. Fora taper such losses would be permissible at 
lengths 8 - 9 times greater. The plot of the function a(z) is shown in Fig. 2. The breaks 
at the ends are equal to 5.4 and1.6°. In order to achieve a transition with minimum loss 
x is necessary to use another function from Ref. 6 or other works on the theory of long 
ines. 
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2. NUMERICAL METHOD OF DETERMINING THE TRANSFORMATION COEFFICIENT 
IN THE ABSENCE OF A CRITICAL CROSS-SECTION 


Let us now calculate the transformation coefficient without admitting the simplifica- 
tions introduced in Eq. (2). In calculating the amplitude of the forward wave| D| with the 
chosen generatrix it is convenient to use the system of equations obtained from Eqs. (39) 
~(40) of Ref. 3 by converting to a new independent variable € and new functions b(é) and 


d(é): 


b(&) = Bee, d(&) = De-*e, (12) 


The quantities b and d represent the amplitudes of the forward and backward parasitic 
waves. The amplitude of the Ho1 wave is considered to be equal to unity over the entire 
transition. 

Differentiation with respect to € will be indicated by a dot. The system of equations 
then takes the form 


cae Be ei. ies Caths Phi) 5 
d + ih,d an Th = : swage a ee 
baited SE fea Bott 
The boundary conditions for this system have the form 
d(0)=0, b(1)=0. (13b) 


Direct solution of the problem in Eqs. (13) leads to a number of difficulties, for Eq. 
(13b) includes conditions at different ends of the interval. Equation (13a) is a linear system 
and hence may be reduced to two Cauchy problems [8]. In the first problem the system of 
equations coincides with Eq. (13a) and in the second it differs from Eq. (13a) by the absence 
of the right-handmembers. The initial conditions of both problems are 


d(ty)=1, b(1)=0. (14) 


The solution of the system of Eqs. (13) is given in terms of the solutions d(1), b(1) and 
d(2), b(2) of these Cauchy problems by 


es d) 0) 49 (15a 
d(&)=d 5) a0) B), ) 
yore el Oise) 15b 
b(§) =b" (&) To. (15b) 


In accordance with Eq. (15) the amplitude of the forward parasitic wave Ip|= |a(a)| is 


d (0) 
ao : 16 
A) (0) ( ) 


|p| =|1 


Upon separation of the real and imaginary members, the given systems of equations 
lead to systems of four equations which have been solved by us on the BESM [high-speed 
electronic computer] by the Runge-Kutta method. Figure 3 gives the results of a calcula- 
tion for the transition whose generatrix is shown in Fig. 2. The solid line is plotted from 
the results obtained on the BESM by use of Eqs. (15) and (16) and the broken line is plotted 
from results of calculation with Eq. (11). Comparison of the curves shows that the assump- 
tion made in Section 1 is valid. This confirms the appropriateness of the method chosen 
for selecting the shape of the generatrix. Over a wide range (even at the wavelength at 
which the critical cross-section arises) the losses will be but a few decibels higher than in 
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Eq. (11). However, it will 

Di/\ing| not reflect in detail the 

| | shape of the curve for the 
dependence of |D|on A in 
Eq. (11). 

The dash-dot curve 
is plotted from the ampli- 
tude of the backward para- 
sitic wave |b(0)|calculated 
from Eq. (15b). As is to be 
expected from general con- 
siderations, the backward 
wave under these conditions 
is extremely small. 
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3. NUMERICAL METHOD OF DETERMINING THE TRANSFORMATION COEFFICIENT 
IN THE PRESENCE OF THE CRITICAL CROSS-SECTION 


Let us assume at a certain cross-section z =2 that the wave number of the H02 wave 
vanishes: h2 (z) = 0. The method of direct numerical calculation from Eq. (13), as used in 
the preceding section, is not applicable here, for the coefficients in (13a) have singulari- 
ties at z =z. For the sake of definiteness let us discuss a tapered transition. The conver- 
sion losses in the transition are determined by the amplitude of the parasitic wave reflected 
into the wide waveguide, |B(0) | . From Eqs. (11) - (13) of Ref. 4 for B(0) we may 
arrive at the expression 


Zz 


~ 7 L 
B (0) = Vas oe ee +) ¢ Path One oid) dg 4 
he (0) 2Vigh, % 


0 


ca Sie Ee L aay: 
ye hh hy dina a EE a neat ( y+ =) \ eS dina V’ -ihV Al 17 
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Here zy, is any point for which zy, < %. Let us choose zy, extremely close to the critical 
cross-section. The functions V and dV/dz are expressed as Airy functions (see Refs. 3 
and 4). The evaluation of the first two integrals in Eq. (17) will reduce to the solution of a 
system of differential equations. Let us denote these integrals by Jj and J2, respectively, 
and introduce the functions 


UW Cm ee (18) 


These functions satisfy an inhomogeneous system of equations in the variable € differ- 
ing from (13a) only in the absence of terms containing hg. The initial conditions of the 
system have the form 


TAO) =O, a) 0. (19) 


As in system (13), we calculated the first two components in Eq. (17) by the Runge-Kutta 
method. 

The integrand in the third component of Eq. (17) decreases rapidly in the region z > Z; 
hence the third integral was calculated by the Gaussian quadrature formula with eight 
abscissas. 

As an example we calculated a series of generatrices for various values of a(0) and 
a(L). The curves for these are given in Fig. 4(CurvesI -IV). The table lists the values of 
| B(0) lobtained for these transitions by the method described. The calculations were per- 
formed for X =8 mm. 
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Figure 4 


As is seenfrom the table, in certain cases the form of the generatrix as obtained by the method 
described in Section 1 may prove unsatisfactory in the presence of the critical cross-section. 
Analysis shows that the first integral in Equation (17) has greater effect on |B( 0)| than 
other components. Hence it is sometimes necessary to choose a(z) so that this integral 
and not integral (2) is as small as possible. 
For this purpose let us introduce the function F(a), defined by the condition 


dF = @ uth 
4 2V Iyhy’ 


(20) 


F(a) may be determined easily from Equation (20), but due to its unwieldiness we shall not 
write it out. Let us further choose a(&) so that F(a) is proportional to the function chosen 
[optimum for Eq. (1)]In W(€). Here the variable € is determined as in Eqs. (3) and (4) by 
replacing L with zy. We will note that the expression for F(a) includes frequency (wave 
numbers hj and hg depend on k) and we choose the most effective frequency of the band. 
For other frequencies of the band a transition of modified form will not, generally speak- 
ing, be better than the transition obtained in Section 1. 

Having determined F(a) and F(€), let us determine the dependence of a on €. The de- 
pendence of z on € is determined from the condition 


dz the a (§) dé 
a T: Lay +h; ‘ (21) 
ZL 2 
jap 7 8) 4 


Knowing a(z), we can, as in the preceding case, find the length zj,. The form of the gen- 
eratrix for z>zy, affects only the value of the third component in Eq. (17) and the quantity 
| B(0) | for wavelengths shorter than that for which the optimum a(z)* (for our purposes) 
was determined. 

In Fig. 4 Curves II-V are given for the same conditions, while Curve V is derived by 
the method described. As seen from the figure, with the exception of the region adjacent 
to the critical cross-section, for the wavelength chosen (A = 8 mm) it almost completely 
repeats curve II]. Hence, a transition with this form of generatrix yields for other wave- 
lengths approximately the same result as a transition with generatrix III. The value of 
) B(0)Jobtained for curve V was 3.4%. 

It is obvious that the method proposed here, beginning with determination of F(a) from 
Eq. (20), may also be used for designing transitions in the absence of a critical cross- 
section. However, satisfactory results may also be obtained by the simpler method pro- 
posed in Section 1. 


CONCLUSION 
The calculation of losses for a chosen form of a(z) in the absence of a critical cross- 


section for each value of A requires approximately 5 minutes on the BESM, while in the 
presence of a critical cross-section it requires 10 minutes. The design of wideband 


* In practical calculations we chose zy, so that the condition |?! <WN |s|S>1, was 
fulfilled, where t and N are determined from Ref. 3. 
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waveguide transitions with extremely low conversion loss (on the order of 40-950 db) and 
relatively short length probably requires the testing of several forms of a(z) and perhaps 
also consideration of the change in amplitude of the Ho wave along the transition, and 
hence will require several hours of machine time. 

We take this opportunity to express our thanks to co-workers of the theoretical- 
design group of IRE AN SSSR for their assistance in the calculations. 
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SYNTHESIS OF STEPPED DIRECTIONAL COUPLERS 


A. L. Fel'dshteyn 


The paper discusses the theory of directional couplers in the form of a chain of sec- 
tions of coupled lines. The system possesses theoretically ideal directivity and matching 
within an infinite bandwidth and optimum (Chebyshev or maximally linear) frequency S 
response in the division of power between the branches. Also discussed are several 


properties of a system having an infinite number of steps. Several laborious i es 
culations are tabulated. design cal 
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INTRODUCTION 


A section of uniform coupled lines may serve as a directional coupler within a wide range 
of frequencies [1, 2]. In a cascade connection of several such sections the bandwidth of the 
system is increased [3]. It is of interest to examine the theory of the synthesis of optimum 
cascade connection of coupled lines. Criteria of optimality may differ; hence we shall dis- 
cuss two cases: a directional coupler with minimum number of sections (Chebyshev or iso- 
extremal characteristic) and a coupler with maximally a linear characteristic. The physically 
attainable properties of directional systems are also discussed, since they are intimately 
related to the synthesis problem. 

Below we present the principal assumptions of this theory and several of its applications. 


1. SINGLE-STEP DIRECTIONAL COUPLER 


Let us examine a uniform symmetrically coupled line (Figure 1) the equations of which, 
according to Reference 1, have the form* 


a U, =U, cosml + 7 (pI, + rJ,4) sin ml, 
(Se U,=U,cosml +7 (rl, + pl,) sin ml, 
1 
2 LT a, 1 =1,cos ml + j (= — +) sin ml, @) 
—_— 
Figure 1 I,=I,cosml +] ee — sin ml, 


\ / 


where p and w are the characteristic impedances and r and v are the current and voltage 
coupling impedances respectively. 
For such a line, as is known [1], the following relationships hold 


Peale Fe e° = lie 5 i) == lon. (2) 


a” ? 


Let us assume that the load resistances at all outputs of the eight-terminal network are 
identical, R = 1 ohm, and let us introduce the condition of ideal directivity and matching of the 
system 


mrvu=R=1. (3) 


Under these conditions the wave transmission matrix of a single-step coupled line has the 
form 


abcd 
i en / Sites a 
aj eel 
mno p 
cos 9 + jpsin® 0 0 —jrsin@ 
o~ 0 cos8-+jpsinO0 —yjrsin®@ 0 (4) 
i 0 jrsin®@ cos ® — jpsin@ 0 
jrsin® 0 0 cos 9 — jpsin OJ 
where r and p are normalized (r ~ r/R and p ~ P/R). From Equation (4) it is easy to find 
the desired scattering matrix elements of one step: 
Si = 0, Sein => 0, (5) 


*The notation is changed slightly. 
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1 S jr sin ® (6) 
cos@+jpsin® ’ 12°" “cosO+ jp sin® * 


Sig — 


2. MULTI-STEP DIRECTIONAL COUPLER 


A cascade connection of any number of uniform coupled lines with balanced coupling (rv = 1), 
discussed in the preceding section, will be referred to as a multi-step directional coupler. Such 
a system is ideally directive and matched over an infinite bandwidth. Thus, confronted with the 
problem of the synthesis of an optimum multi-step directional coupler, the only parameter re- 
quiring optimization is the power division factor between branches 2 and 3 of the system: 


ro 


(7) 


2% | Syof 
1 pol” 


In the synthesis of this frequency function it is necessary to keep in mind the following. 
1. The power division factor is numerically equal to the square of the absolute value of 
element M(T41) of the wave transmission matrix [T] of the system*. 


2. The power division factor is a polynomial in cos2O**:; 


1] = a + a, cos? 6 + a, cost + ... + een cos?” 6, (8) 
where 
1 
8= ml = 2n—- > (9) 


is the length of the step; @j are real constants. 
Due to the unitarity of the matrix [S] 


| Si Poe PS Pees eu) 


and, consequently, the quantity 7 fully defines the system. Let it be required that the function 
7 (@) have the form 


ne 5 om2 (CSO 11 
n=6 er (), oO 


where 8 is a constant; h and S are amplitude and scale factors; 


T,(£) = cosnarc cos x (12) 


is a Chebyshev polynomial of the first kind and of the nth order. The approximate shape of 
the frequency characteristic corresponding to Eq. (11) is shown in Figure 2. 
The relation between the parameters of a multi-step system follows from an evident condi- 
tion: with 9 =07 =0. Then from Equations (11) and (12) we obtain an expression defining 
the attainable passband of a multi-step system: 


s=—- (13) 


4 
ch — arch & 
n h 


Thus, the bandwidth of constant power division (the passband") in a directional coupler is 
fully determined by the number of sections n, the nominal level of power division B2 and the 
allowable deviation from the given level h2. 


*The matrix elements [T] of the system are designated by italic capitals in the same sequence 
as given in Equation (4) for a single step. 
**See Appendix 1. 


p12 S46 676 9 we 


Figure 3 


After the scale coefficient S has been determined from Equation (13), the band overlap 
factor X is found from the equation 


_ Ay  m—arecosS 
K Ma ~arccos SS? (14) 


where 9 and Aq are the cutoff wavelengths of the system with a given tolerance. Equations 
(13) and (14) are illustrated by the plots in Figure 3 for n = 2 and 3. These curves bound the 
physically attainable bandwidth of a directional coupler with given parameters n, B and h. 
From analysis of Equation (11) it follows that the length of a single step in the system is 
defined by 
7 1 take 
mo (15) 


In a number of cases it is of interest to consider a number of cases it is of interest to 
consider a multi-step directional coupler with a maximally linear power division frequency 
characteristic (Figure 4). Let us express this characteristic in the form 


1 = B2 —B2 cos" O (16) 


(a special case of this relationship - the simplest single-step (n = 1) system - was discussed in 
Section 1). 

Calculation of the coupling impedances of the coupler steps in both of the types mentioned 
may be performed by the method of undetermined coefficients, as was done in References 4 and 
5 for stepped transitions. The design formulas are derived by comparing the coefficients of 
the same powers of sin? © in Equation (11) or (16) and in the expression for this function as 
determined by direct multiplication of the matrices [T] of the steps of the directional coupler. 
The explanation is given in Appendix 2. This same appendix contains tables of values of r 
calculated by this method for the case of a two-step ideal directional coupler with a Chebyshev 
power division characteristic. 


3. LIMITING MULTI-STEP SYSTEM 


It was shown above that a multi-step directional coupler has a periodic frequency charac- 
teristic (Figures 2 and 4) and, consequently, that the passband of such a coupler is always 
limited at the low- and high-frequency ends. There may arise cases in which a constant 


Figure 4 Figure 5 


207 


transition loss within the coupler is required within an unlimited range of frequencies. An 
excellent approximation to such a characteristic is obtained in the case of the limiting multi- 
step system, in other words, in the case of coupled tapered lines. Let us use the limiting 
relation of the first approximation for the power-division constant of loosely-coupled non- 


uniform lines [6] 
] * 
n= Sip = jm|r (ye dy, 


Q) 


(17) 


where r(y) is the coupling distribution function; y =/-x is a coordinate read from the origin 
(Figure 5). Assuming that the coupling changes exponentially 


r(y) = re” 


and integrating in Equation (17), we have 
(18) 


Uis=fots Sam (A nie Jen cos mt. 


A typical characteristic plotted from this equation is shown in Figure 6. The broken line 
represents the limit curve (with / +o). It is evident from the curves that the passband of a 
directional coupler in the form of tapered lines is 
limited only at low frequencies. At this point we may 
compare stepped transitions and tapered transitions: 
while the first have a periodic characteristic and thei1 
passband is limited on both sides, the passband of the 
latter is limited only on the low-frequency side. The 
comparison may be extended as follows: with a given 
finite passband a stepped system will always be shorte 
than a tapered system, as is the case with ordinary nonuniform lines [6]. 


Figure 6 


CONC LUSIONS 


Modern synthesis of circuits is based on frequency characteristics in the form of a Cheby- 
shev polynomial, Zolotarev fractions and the so-called maximally linear function. In recent 
years this method has found wide application in microwave techniques. Particular interest 
attaches to the case where synthesis is applied directly, without equivalent circuits — by an 
appropriate substitution for a frequency variable ("frequency conversion"). In this respect, 
devices consisting of uniform line sections of identical length are extremely convenient and 
are successfully applied in the synthesis of matching transformers and filters [7, 8]. 

In the present paper this method was extended to directional couplers. These devices, 
possessing (with a proper choice of parameters) ideal directivity and matching over an un- 
limited range of frequencies, have a polynomial frequency characteristic of power division 
between branches. This polynomial relationship permits evolving a synthesis for directional 
couplers with Chebyshev and maximally-linear characteristics which in many respects is 
similar to the synthesis of stepped transitions and stepped filters. 


APPENDIX 1 
It follows from Reference 9 that 
Sy  BI—AM 
a i (1) 


but it is evident from Equation (4) that for our case 


I1=0, A=F, i) 
consequently , 

Sie 

Ses eS iil 

Sis (III) 
and 


| Sie |? + | Sis? = 4. 
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Thus, the element M of the wave transmission matrix wholly defines the power division be- 
tween branches 2 and 3. Let us now investigate the nature of the frequency dependence of this 
element. It is evident from Equation (4) that in multiplying any number of such matrices each 
element of the resulting matrix is the sum of the products of the following terms: 


7P, sin ®, cos@-+ jp; sin 0, 


(iV) 
that is, the possible combinations 
sin” ®, cos™ 0, sin” 6 cos™@. (Vv) 
We are interested in the modulus of Equation (11); if 
S 
a =a+]y, 
then 
= es 
(Sanaa ae 
that is, terms of type (V) acquire the form of even functions 
sin?” @, cos?” @, sin?" @ cos?” @. (VI) 
All even powers of sines and cosines may, after elementary conversion, be expressed 
in terms of cos26 or sin20. Thus, the power division constant is a polynomial in cos20: 
1 = % + a cos? 6 + a4 costO +... +4, cos?”"O. (VII) 
APPENDIX 2 
Multiplying the matrices of two steps, we find 
M =asin?0-+ 7ksin @ cos 0, (VIII) 
where the following designations are used: 
& = pilo — Par; 
k= Ty + To. ay 
Hence, 
H =|M |? = sin? 0 (k? — (k? — a2) sin? 0}. 
(XI) 
On the other hand, it is necessary that 
; iC) : » 1 ems 
n = B—WT2 (“) Se) eh ee sin? Oh. (XII) 
Comparing the coefficients of the same powers of sin2© in Equations (XI) and (XII), we ob- 
tain the system of equations 
2 4 ; 
{ + raja ah (= — a) a 
ln tr)? — (nar — nes = 4 2 pabis) 
Considering that —_ 
a= M1 + a Pz V1 + res (XV) 
we obtain after conversion 
(XVI) 
A (rire)? + 2B (ryr2) — BY = 0, 
where ; } vi 
A= 408 (=~); (XVII) 
Beeps (XVIII) 
St 
Solution of Equation (XVI) yields 
(XIX) 


naa ee V wach yrza). 
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where 


Bali ul 1 
deus te (Se 

Allowance for | My: | Overlap 
transition loss 2D a PUNT coefficient 

5, db rl r2 A2/A1 

0.1 S290 0.2311 2.3461 

0.2 1.3515 0.2609 2.7790 

0.3 LS Od 0.2911 Sealant 

0.5 VE S8T96 0.3235 3.6897 

LO SSA; 0.3970 4.8500 


We present a table of values of r;.and rg for a two-step Chebyshev directional coupler with 
equal power division between branches 2 and 3 (79 = 1, transition loss Lg = 3 db). 
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ATTENUATION MEASUREMENTS IN WAVEGUIDES 


Yu. N. Kazantsev 


This paper proposes a method for the measurement of attenuation in waveguide elements 
by determining the standing-wave ratio in the waveguide in front of a cavity formed by the piex 
to be measured, a diaphragm and a shorting plunger, at resonance. Measurement setups are 
presented for the Hg; wave in circular waveguides and formulas are given for error determin: 
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1. MEASUREMENT PROCEDURE 


The resonant cavity method is usually used in measuring attenuation in relatively short 
aveguide sections. This method consists essentially of the measurement of the Q of the cavity 
ormed by the waveguide section under investigation and two metal end plates. The Q is then 
recalculated taking into account the effect of losses in the end portions on the attenuation in 
the section. 

Methods of Q determination at microwave frequencies may be divided into two classes: (a) 
determining the Q from the duration of free oscillations; (b) determining the Q from the width of 
the resonance curve. 

The method of determining the Q from the duration of free oscillations, described in Reference 
I, is suitable for systems with extremely high Q; with a decrease in Q the measurement error 
increases considerably. 

Methods of determining Q from the width of the resonance curve find their widest application 
in microwave techniques [2, 3]. However, in a number of cases they prove unsuitable owing 
to the high stability requirements for the signal generator. 

The method described below permits the measurement of an extremely wide range of attenua- 
tion values with the use of simple equipment. In addition, the method does not place high stability 
requirements on the signal generator. Attenuation within the section is determined from the 
relative power absorbed in the cavity at resonance. 

Let us determine the relationship between the value of this power and the parameters of 
the cavity resonator. For this purpose we shall examine the cavity formed by diaphragm 1 and 
a shorting plunger 2 placed within a waveguide 3 which is matched with the signal generator 
(Figure 1). 

The diaphragm must pass the given wave 
mode without exciting other, unwanted modes. 
Under conditions near resonance the 

field within the cavity may be regarded as 

the sum of two opposed traveling waves of 
equal amplitude. The expression for the Figure 1. Diaphragm and plunger 
amplitude of these waves was derived in in waveguide. 

Reference 4 by the method of superposition 

of waves multi-reflected from the end faces. At the resonant frequency, after transformation, 
this expression takes the form 


RARVVAVTAVAVagaa 


SAAS SSS AAS SAS) 
iY, 
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as ! 

u "i (1) 
where A is the amplitude of the traveling-wave field to the right of the diaphragm in the case 
where, instead of a shorting plunger, there is connected a matching load; 7 gp K 1 is a quantity 
indicating that part of the power which is lost by the wave within the cavity (the quantity 7) gp in- 
cludes both the thermal losses and radiation losses within the waveguide). 

Before proceeding to the equation defining the power absorbed in the cavity, let us explain 
the physical meaning of certain quantities. If we designate the power of the incident wave as Py and 
the transfer constant of the diaphragm was 1, then the product Pj” 1 represents the power of 
the wave passing through the diaphragm in traveling-wave operation (diaphragm within a waveguide 
matched with the generator and load). The ratio u/A indicates the number of times by which 
the wave amplitude within the cavity has increased as the result of multiple reflections from 
the end faces; the difference 7 g,-71 defines the thermal losses in the cavity. 

Now it is easily seen that the power P absorbed in the cavity is defined by 


P=P\m & y (Map m1) (2) 
Substituting Equation (1) in Equation (2), we obtain 
pebep eens (3) 


2 


‘dp 
Assuming in Equation (3) that P=P ,, let us determine the condition for the total absorption 
of the power of the wave reaching the cavity from the waveguide: 


sa) oye (4) 


From Equation (3) we can also easily express the thermal loss 


20d 


1 = Ngp— th (5) 


in terms of the diaphragm transfer constant and the relative power p = P/P, absorbed in 
the cavity: 


n=m|(+—1)+—2V1—9]. (6) 


Equation (6) can be written in simpler form by expressing the quantity p in terms of the 
standing-wave ratio in the waveguide in front of the cavity. Since p = 1-k2 (k is the modulus of 
the reflection coefficient of the cavity) and p = (1+k) / (1-k), then 


= 1? (7) 
for n>, [plus sign in Equation (6)] and 
n=m/P (7) 
for <n [minus sign in Equation (6)]. 
From Equations (7) it is evident that for the determination of the quantity 7 it suffices to 
know the transfer constant of the diaphragm in traveling-wave operation as well as the standing- 


wave ratio in front of the cavity at resonance. The quantity 7 includes losses in the side walls 
of the cavity 7,4, losses in the plunger 7,1 and the diaphragm +9 of the wave within the cavity | 


Y= Ner — Nr + ro. (8) 


The principal part of the setup for attenuation measurement is the cavity, which consists 
of the two end faces of the diaphragm and plunger and a waveguide section of one or two half- 


AB 


wavelengths. The cavity portion may include a test piece with length of aoa where n is an 


integer and \p is the waveguide wavelength. 
The number 7 is determined for two cases: (a) the piece is included (7 p) and (b) the piece 
absent (70); then the losses in the piece are defined as the difference of 7 p and 7 9*: 


1 piece = Nz — No- (9) 


From Equations (7) and (9) we can write the following expressions for the attenuation within 
the piece (in decibels): 


N piece = 2,17 41 (P —Po) With ny < Ho <p, 


Oech Oy —f A, 

piece = wu is with 1 > Ns > No, (10) 
2.47 1 (ep>—t ; 

1] piece = lee vith No< 1 <p, 


Po 


where p and Pq are the standing-wave ratios in front of the cavity with and without the connecti« 
of the piece, respectively. These expressions, as in the case of the initial Equation (1), are 
valid under the condition that 7 gp = Npiece + NO +74 «1, which is easily satisfied. 

From Equations (10) it is easy to derive the equations for calculation of the total relative er 
in measurement: 


Anpiece Am hogs au Ne 

= : T - j 4% ae aa 
Npiecet;, Mh. 0 sate PiaPe with 1 < No < Np, 
Anpi A A oAo 

piece Am , — Podp pApo i * 
As wi 

"piece " p(Po— p) * £o(Po— P) th 91> > No, (11) 
Anpiece An PoAp Ady : 

= ie arial with aa 
Npiece ™ eh Rea) No <1 < Np. 


We will note that the case 7g < 71 <n RB corresponds to the maximum range of measurable 
attenuation at one and the same maximum standing-wave ratio. 

Below we describe two variants of the setup for attenuation measurement of the Ho 1 wave 
in circular waveguides. 


*If the length of the test piece differs considerably from n AB there is introduced into the 


measurements an additional error, which is extremely small under the condition that n > 1. 
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2. MEASUREMENT SETUP 


Variant 1. Indication is achieved by means of a directional coupler within a circular 
waveguide. 
The block diagram of the measurement setup is shown in Figure 2a. The signal from 
_klystron generator 1, modulated by a sawtooth voltage, is applied through decoupling attenuators 
3 and 3' and calibrated attenuator 2 to an Ho 1 wave exciter in circular waveguide 4 [5, 6]. 


10 


a, he V3 4 wv) 7 b 
Figure 2. Block diagram of measurement setup. 


a) First variant; b) second variant 


Passing through helical filter 5, a backward-wave directional coupler in circular waveguide 
6 [7, 8] and telescoping junction 7, the Ho; wave reaches the diaphragm with an aperture system 8. 
The diaphragm excites the cavity of test piece 9. Part of the power of the backward wave 
from detector 11 is applied to oscillograph 12. In tuning the cavity by means of plunger 10 through 
the klystron range (reproduced on the oscillograph screen) a sharp dip appears (Figure 3a). The 
standing-wave ratio in the circular waveguide in front of the diaphragm is easily determined by” 
comparing (by means of the calibrated attenuator 2) the height of the response with the level of 


the dip. * 


Figure 3. Oscilloscope pattern in tuning. 


Variant 2. Indication is achieved by means of a directional coupler in a rectangular waveguide. 
The block diagram of the measurement setup for this variant is shown in Figure 2b. In this 
case the backward-wave directional coupler is placed within the rectangular waveguide before the 


Ho 1 wave exciter. 


*In practice the klystron signal is conveniently modulated by rectangular pulses. Then at 
resonance, due to parasitic frequency modulation, only part of the pulse "collapses" (Figure 3b). 
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This variant is simpler in construction, for it does not require a directional coupler within 
a circular waveguide, but the accuracy of measurements, as shown below, is lower than in 
variant 1. 


3. ATTENUATION MEASUREMENTS IN WAVEGUIDE ELEMENTS 


The above method has been used in attenuation measurements in straight pieces as well as 
in curved pieces. Attenuation of the Hg; wave has been determined in periodic (annular and 
helical) waveguides with dielectric and metallic sheaths [9, 10]. Attenuation of the Hp; wave in 
curved sections of these waveguides has also been investigated [11]. 

Let us dwell in greater detail on the process of measurement. As has been stated, attenua 
within the cavity due to the two kinds of dissipation of the wave is determined from the standing: 
wave ratio in the waveguide in front of the diaphragm in resonant operation. The standing-wave 
ratio is measured by comparing the height of the pulse on the oscillograph screen with the level 
of the dip. Due to the presence of reflections from the H,, wave exciter and attenuators, and 
also due to the fact that the backward-wave coupler has finite directivity, the height of the 
pulse and of the dip vary upon displacement of the diaphragm and the cavity by means of the 
telescopic junction. 

Let us examine the influence of the above factors on the results of measurement of the 
standing-wave ratio (SWR) or the coefficient of reflection from the cavity for both variants of 
the measurement setup. Since the SWR is easily expressed in terms of the reflection coefficien 
subsequent discussion will deal with the reflection coefficient k. 

In the circuit shown in Figure 2a the power of the wave at the detector on the directional 
coupler will vary with displacement of the diaphragm from 


A’ = B (kh? + 6) (12’) 
A” — B (k? <3 61), (12”) 
where Bp is a constant for the given setup, while 6 4 depends on the coefficients of reflection 
from the cavity and the exciter from the direction of the circular waveguide as well as on the 


directivity coefficient of the coupler. 
The relationships in Equations (12) are valid when 


kip 1, bn <k,. (13) 


where kp is the modulus of the coefficient of reflection from the exciter and ky is the power 
directivity of the coupler*. 

If it is considered that in the absence of resonance practically all the power is reflected 
from the diaphragm and, consequently, the absolute value of the reflection coefficient is equal 
S ae then in the case of resonance the absolute value of the corresponding reflection coef- 

icient kis 


to 


——— 


Iie V/ aa (14) 


; “ies 
a 
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where A' and A" are the maximum and minimum values of the backward-wave power at resonan 
at the detector in the directional coupler; ) and Aj are the corresponding values in the absenc 
of resonance. 

From Equation (14) it is seen that in order to determine the reflection coefficient it is 
sufficient to compare the values of A' and A" (dip level) and Af and § (pulse level). As was 
stated above, this comparison is performed by means of a precisely calibrated attenuator. 

In the circuit shown in Figure 2b the coefficient of reflection from the diaphragm within 
the circular waveguide at resonance is determined by comparison of the reflected waves in 
the rectangular waveguide up to the Hg; exciter both at resonance and off resonance. 

Upon displacement of the diaphragm, the wave power at the detector in the directional 
coupler will vary from 


*By kp is meant the absolute value of the total reflection coefficient resulting from the 
presence of the exciter and from incomplete matching of the decoupling attenuator (3') with the 


line; ky takes into account the coefficient of reflection from the power-matching termination of 
the directional coupler. 
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to C’ = D (hk + 6.) (15’) 
CY = D(k? — 6,), (157) 
_ where Dis a constant for the given setup; 6 9 depends on the coefficient of reflection from the 


cavity and the decoupling attenuator as well as on the directivity of the coupler and the charac- 
teristics of the Hg; wave exciter. 


Equations (15) are valid when 


sy< oe (by Bs << ile Te << sick, ka <<! (16) 


where S11, S99 ands 9 are the scattering matrix elements of the exciter-filter junction [12]; 
_ kg is the modulus of the coefficient of reflection from the decoupling attenuator. 
If it is considered that in the absence of resonance practically all the power is reflected from 
the diaphragm, then on the basis of Equations (15) we can write the following expression for 
the modulus of the coefficient of reflection from the cavity at resonance: 


i (GEOR oi ~ 
jo vee ; (4 7) 
V are 


where C' and C'' are the maximum and minimum values of backward-wave power at the detector 
in the directional coupler at resonance; 6 and Co are the corresponding values in the absence 
of resonance. 

Comparing the values of C', C"', Q and Ci by means of the calibrated attenuator, on the 
basis of Equation (17) we can easily determine She reflection coefficient. 

Since it is more difficult to satisfy conditions (16) than the analogous conditions (13), the 
accuracy in determining the reflection coefficient for the second variant will be lower than for the 
first variant. 


4, SELECTION OF DIAPHRAGM TRANSFER CONSTANT AND 
LENGTH OF TEST PIECE 


From Equations (7), (9) and (10) it is evident that the SWR (p and p g) depend on the losses 
in the test piece and the end faces as well as on the transfer constant of the diaphragm. The 
minimum value of SWR must satisfy conditions (13) or (16) and, for convenience in measurement, 
the maximum value must not exceed 5-6. Let us discuss the sequence of operations in the selectior 
of the diaphragm transfer constant and the length of the test piece, for example, for the case 
no0<"1<Np- First the cavity is assembled without the test piece and, gradually increasing 
the diaphragm transfer constant, the SWR before the cavity is measured by the method described 
above. A typical dependence of SWR Po on the diaphragm transfer constant for the case 7 9 = 0. 002 
is shown in Figure 4. From this curve it is seen that with an increase in the transfer constant 
the SWR at first decreases and, upon reaching unity, begins to increase. 
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Figure 4. Standing-wave ratio as a function 
of diaphragm transfer constant. 


The transfer constant 7; must be so chosen that the value of pg is the maximum permissible 
value (the right handbranch of the curve in Figure 4). 

The length of the test piece is determined by the condition that the SWR in front of the cavity 
with the test piece inserted is of the same order as without the test piece (in gradually increasing 
the length of the test piece, for the chosen transfer constant 7 4 the SWR will decrease at first and, | 
upon reaching unity, begin to increase). 

It frequently occurs in practice that the length of the test piece is given and cannot be varied. | 
In this case the diaphragm transfer constant and initial attenuation at the end faces 7 g are so chose 
that the SWR, not exceeding 5-6, satisfies conditions (13) or (16). 
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5. THE DIAPHRAGM AND ITS CALIBRATION 


The diaphragm is a thin waveguide partition with one (Figure 5a) or several (Figure oy, ©) 
apertures. The number of circular apertures or slits is so chosen that upon incidence of the Hg: 
wave within the cavity no asymmetrical waves are excited. However, such diaphragms may, in 
addition to the Hg; wave, give rise to symmetrical magnetic waves of higher order. Hence, 
diaphragms a and b may be used in the case where within the test waveguide there may be propa- 
gated a single symmetrical magnetic Hj, wave (Ka < 7.02), while diaphragm c may even be used 
in the case where within the waveguide, in addition to the Hg, wave, there may be propagateda | 
symmetrical magnetic wave of the second order, that is, Hgo (Ka < 10.2). The coupling aperturi 
are located on the radius where currents induced by the Hgg wave are equal to zero. Here Kis 
the free-space propagation constant and a is the waveguide radius. 


Figure 5. Various diaphragms. 


Calibration of a diaphragm consists of determining the transfer constant 7, within a range 
of frequencies. The block diagram of the setup for diaphragm calibration is given in Figure 6. 


Figure 6. Block diagram of setup for diaphragm calibration. 


The signal from klystron generator 1, modulated by rectangular pulses, is applied through 
decoupling attenuators 3 and 3' and calibrated attenuator 2 to an Ho; wave exciter in circular 
waveguide 4. The Ho; wave, passing through helical filter 5, reaches diaphragm 6, on both 
sides of which are located telescoping junctions 7 which permit eliminating reflection from the 
generator and receiving ends of the setup. Then (through helical filter 8, second exciter 9 and 
decoupling attenuator 10) the signal reaches detector 11 and, finally, measuring amplifier 12. 

The transfer constant of the diaphragm is determined by comparing the signal level at 
detector 11 in the presence of the diaphragm with the signal level at detector 11 with the dia- 
phragm removed. It is easily seen that reflections from both exciters introduce errors in 
determining the transfer constant. 

In order to eliminate this error telescoping junctions 7 are used. With changes in distance 
from the diaphragm location to the Ho, wave exciter the amplitude of the wave reaching detector 
11 varies from 


aed A ; 
Umax = i i to “min= TB (18’) 


in the absence of the diaphragm and from 


AY AV 


Wd max= > - +O ie 
1 — hy — ky + hy hey dmin Pa ae a Rian 


(18”) 


in the presence of the diaphragm, where A is a constant, ky « 1 and ky « 1 are the absolute 
values of the coefficient of reflection from the exciters*. 


*k; and ky represent the absolute values of the total reflection coefficients from the directio 


of the circular waveguide as caused by the presence of the Hg; exciter and incomplete matching 
of the decoupling attenuators 3' and 10 with the line. 
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On the basis of Equations (18), ignoring second order terms, we can write the following 
xpression for the transfer constant of the diaphragm: 


1 /“dmax , “dmin \? 
=> — | —— ——————— . 19 
: 4 ( Umax ij a) i) 
min min 


| Thus, in order to determine the transfer constant it is sufficient to compare the amplitudes 
of ug ax and ug min we Umax ~ Umin> this comparison is performed by means of the cali- 
brated ¢ attenuator. 


a 


CONCLUSION 


The proposed method for the measurement of attenuation permits considerable simplification 
of the measurement setup. In using this method it is not necessary to measure the width of the 
resonance curve, which eliminates the need for high signal generator stability, high linearity of 
the oscillograph sweep, etc. Although the setups described above dealt with the testing of cir- 
cular waveguide elements, the proposed method can be used in measuring losses in waveguides 
of other configurations. 
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SOME PROBLEMS OF MICROWAVE CHANNEL MATCHING 
IN RADIO SPECTROSCOPES 


V.I. Muromtsev and A.K. Piskunov 


The paper discusses special aspects of matching the microwave channel of a radio spectro- 
scope with a reflex resonator. Analyses are given for possible nonlinear distortions of the 
paramagnetic absorption signal which arise under certain conditions in the matching of the 
resonator with the waveguide system. The influence of the dispersion signal on the shape of the 
absorption signal is evaluated. 


INTRODUCTION 


Analysis of the operation of various types of radio spectroscope has been presented in 
References 1-3. In Reference 1 the chief causes of limited sensitivity of spectroscopes were 
examined in detail and the optimum conditions for matching of the microwave channel with the 
resonator were determined. The theory of simple and superheterodyne spectroscopes and the 
influence of dispersion on the observed shape of the paramagnetic absorption signal in a re-entrant 
resonator were first discussed in Reference 3. 

The present paper presents a detailed analysis of some specific aspects of matching of 
the microwave system of spectroscopes with the reflex resonator, possible nonlinear distortions 
of the signal shape, and evaluation of the influence of the dispersion signal on the shape of 
the absorption signal. 

In order to increase spectroscope sensitivity it is necessary to match the microwave wave- 
guide system with the resonator. The optimum matching condition depends not only on the type 
of spectroscope and the detector characteristics, but also on the internal noise of the generator, 
the noise figure of the amplifier and the noise temperature of the detector or converter. 

It is most important to know the optimum matching conditions from the point of view of 
obtaining the maximum signal-to-noise ratio. As a first rough approximation it is often useful 
to know the optimum matching condition at which the maximum signal is obtained at the output of 
the system. Calculation of the matching conditions is based on the assumption that the waveguide 
system with a resonator is described by an equivalent lossless long line with characteristic 
resistance Rog loaded by a lumped impedance Z,. 

The resonator is represented as an Beulpeient: tuned circuit and is characterized by an equiva- 
lent loss resistance r, Qo and resonant frequency fo. The coupling of the resonator with the 
system is represented by an ideal transformer with transformation ratio n. 

The lumped impedance loading the long line is defined as the input resistance of the equiva- 
lent circuit shown in Figure 1. 

The reflection coefficient Tis defined by 
vet be oy, mae 


1 
Vine Zi +R, ; ( ) 


iP 


where Uyegs is the amplitude of the reflected 
wave and Uj;n, is the amplitude of the incident 
wave. 

During observation of the absorption signal 
the resonator is tuned to the generator frequency. 
In this case the load impedance Zy, of the long Figure 1. Equivalent circuit of waveguide 
line will be resistive and equal to system with resonator. 
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Z,= RL =-,. (2) 


: oe tuning of the resonator to the generator frequency the reflection coefficient Tis de- 
ined by 


ope 2 
r— Ron 


Mp FE grt (3) 


Paramagnetic absorption of the microwave power causes additional losses in the resonator, 
which is the same as an increase in the equivalent resistance r, characterizing the loss, by Ar. 
- As the result of the change in loss there is a change in amplitude of the reflected wave. Fron 

_ Equations (1) and (3) it follows that 


r— Ryn? 


Orep = ro Ron Vine (4) 


or 
Uref =f (7, n). (5) 
Assuming that the paramagnetic losses are small, the magnitude of change in amplitude of 
the reflected wave may be written in the form 


a 
AU ret = $2 Ar 
or 
Ron? 
As follows from Equation (6), at one and the same value of paramagnetic absorption (Ar) the 


change in amplitude of the reflected wave depends on the transformation ratio n¢. 
The magnitude of the signal AU,..¢ will be a maximum when 


AU yet =2 ArVince. (6) 


Ron? 
aaa (7) 


and is given by 


4 A 
AU et = > Vine — - (8) 


In practice it is more convenient_to express AUyer and the optimum coupling condition (7) as 
functions of Ip. Replacing Ron? in Equation (6) by the reflection coefficient [p, we obtain 
from Equation (3) f, 
4 A 
AU ret = = Vine(1 —T5) —. (9) 
Equation (7), defining the condition at which the change in amplitude of the reflected wave will 
be maximum is written in the form 


Ty =0. (10) 


Condition (10) defines the optimum coupling in observing the absorption signal in a reflection 
spectroscope with linear detection. 


1. MATCHING IN SQUARE-LAW DETECTION 


It is known that the matching condition at which the maximum signal is obtained at the output 
depends on the detector characteristic. In the general case the choice of the matching condition 
depends also on the level of the incident power at the resonator. 

Let us determine the condition of optimum matching in square-law detection in a circuit of 


the reflex type in observing an absorption signal. 
In the spectroscopic case the change in the amplitude of the microwave carrier represents 
the electron spin resonance signal. The voltage applied to the detector is 


Vine = Uref + AU) cos at, (11) 
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where Uyeg is the carrier amplitude and AU, is the change in amplitude of the microwave 
voltage caused by paramagnetic absorption. The expression for the signal component after the 
square-law detector has the form 


I, = 4 [(Uye¢ 4+- AU.) cos wt]? = CHUL INU (12) 


Thus, the value of the signal at the output is proportional to the amplitude of the carrier Uyef 
and the value of the signal AU, at the input. Here, for an increase in the transfer constant of the 
detector at the same signal AU,, it is necessary to increase the carrier amplitude. Spectroscope; 
in which in order to increase the detector transfer constant, the carrier level is artifically 
increased are described, for example, in Reference 2. For a spectroscope with a reflex 
resonator, from Equations (1) and (9), we have 


ref = lr Vine, 


AU AU oi = (le) Eee ae 


Substituting Equation (13) in Equation (12), we obtain the dependence of the signal at the output 
on the reflection coefficient Tp: 


ie = const P3(t— Ve eae ; (14) 
Analysis of Equation (14) shows that with 
. { 
DS re (15) 


the signal at the output will be maximum. A plus sign before Tp indicates coupling below critical 
the negative sign indicates coupling above critical. In square-law conversion the signal at the 
output is proportional to the square of the voltage Utne: that is, the incident power at the resonat« 
In spectroscopes employing crystal detectors the optimum matching condition (15) is applicable 
at powers not greater than a fraction of a milliwatt. At high powers Equations (14) and (15) are 
not applicable. 

We will note that with To! less than optimum the decrease in signalatthe output is caused 
by a decrease in the transfer constant of the detector, and at values of |[p| greater than opti- 
mum; adecrease in signal at the output is associated with a decrease in the absorption signal 
AUe, which is caused by a decrease in the microwave power entering the resonator. Matching 
ensuring maximum spectroscopic sensitivity differs from condition (15), which was pointed out 
in Reference 3. 

Results of the investigation performed may be briefly formulated as follows: 

1. The optimum matching condition in square-law detection is defined by Equation (15). 

2, For optimum square-law conversion the optimum matching conditions do not depend 
on the level of power incident at the resonator. 

3. The signal at the output of the square-law detector is proportional to the incident power 
at the resonator. 


2. MATCHINGIN A SIMPL# REFLECTION SPECTROSCOPE 
AT MEDIUM AND HIGH POWER LEVELS 


In the operation of a simple spectroscope with nonsaturable specimens, in order to increase 
sensitivity it is necessary to increase the power level of the generator. At powers of the order 
of several milliwatts, depending on Ip), the detector may operate both in the square-law region 
and in the linear region. Let us discuss the manner in which in this case it is necessary to 
select the reflection coefficient [p . Over a wide range of input voltages the detector charac- 
teristic has a relatively complex shape. Analytical determination of the optimum matching 
conditions is an extremely laborious undertaking. Hence we shall limit the discussion to a 
qualitative examination of the matching conditions at which maximum spectroscope sensitivity 
is obtained. 

In the case under discussion the selection of the reflection coefficient Tp is based on the 
following considerations. In the square-law region of the detector characteristic the transfer 
constant is proportional to the carrier amplitude. With transition to the intermediate region 
the increase in transfer constant is retarded and with further increase in amplitude (into the 
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linear region) the transfer constant is practically independent of carrier amplitude. | 
However, in linear operation with an increase in voltage applied to the detector there is 
a sharp increase in anomalous noise of the detector, particularly in the audio-frequency region. 
At high levels of incident power at the detector even a small deviation of Tp from zero may 
lead to transition from square-law detection to linear detection. Consequently, in contrast to 
the preceding case, the optimum matching condition is determined, on the one hand, by a 
decrease in transfer constant upon a decrease in carrier level and, on the other hand, by an 
increase in anomalous noise of the detector with an increase in the applied voltage. The choice 
of Ip is determined by proper choice of detection conditions. Optimum detection conditions 
depend on the noise figure of the amplifier and on the individual properties of the crystal detector. 
For a simple spectroscope the frequency and amplitude noise of the generator may be disregarded 
Experiments reveai that for spectroscopes with modulation of the magnetic field no higher 
than 1 ke the optimum operating conditions lie between the square-law and linear regions. 
Thus, in a simple spectroscope operating at various power levels, in order to ensure opti- 
ce detection conditions it is necessary to select appropriate values of reflection coefficient 
p: 


3. DISTORTION OF ABSORPTION SIGNAL SHAPE 
IN THE MICROWAVE SYSTEM 


Experience in operation with a simple reflex circuit (particularly with the superheterodyne 
circuit) shows that with the reflection coefficient IT close to zero, the shape of the absorption signal 
is distorted due to the presence of the dispersion signal, even with small deviation of the generator 
frequency from the natural frequency of the resonator. A detailed discussion of this problem 
is presented below. For large absorption signals distortion of shape is observed even with 
precise tuning of the generator to the resonant frequency of the resonator. These distortions 
are not associated with the influence of dispersion. 

Experimentally observed distortions of the electron spin resonance signal shape of diphenyl- 
picryl-hydrazyl depending upon matching conditions with precise tuning of the generator frequency 
are shown in Figures 2-4. The signals were observed on a superheterodyne radio spectroscope 
similar to that described in Reference 3. 

As is seen from these figures, at large signal levels the shape of the absorption curve may 
be considerably distorted. Let us discuss the cause of the observed distortion. The matching 
condition Ip = 0 corresponds to the case of critical coupling. Departure of the reflection 
coefficient from zero may occur in the above-critical direction and in the below-critical 
direction. With the reflection coefficient close to zero the amplitude of the reflected wave is 
defined by the approximate expression 


af RVing (16) 
Ure = (0 — Ron) 


With matching conditions below critical (Ryn?/r < 1) and at critical (Ron?/r = 1) Equation (16) 
is conveniently written in the form 


Ue 
Useg = AR, a ’ (17) 


where 
AR, =r—RwS 0. 


Paramagnetic absorption causes an increase in resistance ARy, by Ar and, consequently, a cor- 
responding increase in the amplitude of the reflected wave 


U: 
Uref = (AR, + Ar) —¢ (18) 


The increase in the amplitude of the reflected wave is the signal. The magnitude of the signal 
is proportional to Ar and, consequently, with matching conditions below critical and at critical 
the shape of the signal is not distorted. 

With matching conditions above critical the paramagnetic absorption will cause a change 
in amplitude of the reflected wave which is defined by the expression 


U.: 
U vet ee | AR, Ar — ’ (19) 


where 


As follows from Equation (19), with an increase Ar the amplitude of Uyef at first decreases, 
at ARp = Ar becomes equal to zero, and with a further increase in Ar begins to increase. 

Thus, in the case of above-critical coupling the proportionality between the paramagnetic 
absorption of microwave power and the change in amplitude of the reflected wave may be 
disturbed and the shape of the signal distorted. In order that the signal shape may not be 


Figure 2. Absorption signal of diphenyl- 
picryl-hydrazyl at below-critical coupling. 


Figure 3. Absorption signal of diphenyl- 
picryl-hydrazyl at above-critical 
coupling with ARp <Ar,ay- 


Figure 4. Absorption signal of diphenyl- 
picryl-hydrazyl at above-critical 
coupling with ARp2 >A R,, but ARR9<Armax. 


AR, = |r — Ron? |. 
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distorted at above-critical coupling it is necessary to increase the mismatch (ARp) or to 
decrease the size of the speciment. 

Let it be noted that, in accordance with Equations (18) and (19), in the transition from 
below-critical matching to above-critical matching the signal phase is reversed. Equations 
(18) and (19) clearly define the observed signals. The signal in Figure 2 corresponds to the 
case of below-critical matching. The signals in Figures 3 and 4 correspond to the case of 
above-critical matching, wherein Tp} for the signal in Figure 3. In both cases (Figures 3 
and 4) ARB< Atmax:- 


4. AMPLITUDE OF REFLECTED WAVE AS A FUNCTION OF DETUNING 
Aw UNDER VARIOUS MATCHING CONDITIONS 


The amplitude of the reflected wave is defined by the expression 


Z, — Ron® 
L oe 
Uret = Zy + Ro® ine) (20) 


where 


Zy=r +i(ol aa): 
Letting X = w L- =~ , the modulus of the reflection coefficient may be written in the form 


= (r = Ron*)?-+ X2 - 
Me (FF Row PE a 


As is known [4], for frequencies near resonance 


x 2A@ Ao 
Te = Qo Oo a A@o ’ (22) 


where wo is the resonant frequency, Qg is the Q of the unloaded resonator, AW g is one-half 
the total 3-db bandwidth; Aw = w - wo is the detuning from resonance. Applying Equation (22) 
the modulus of the reflection coefficient is conveniently expressed in the form 


rare eh ae 
[T= iY Ratio hae (23) 
(1 ae ) +( el 
It follows from Equation (3) that 
Ret.» MV (24) 
r AS ee 


Substituting Equation (24) in Equation (23), we obtain the final expression for dependence of the 
modulus of the reflection coefficient |I'| on detuning under various coupling conditions: 


/ (ae) + (42) 


[P|= Sara eC (29) 
V eres 


Equation (25) is of primary importance in evaluating 
the influence of the dispersion signal and the frequency 
stability of the generator under various matching con- 


ditions in the reflex circuit of the spectroscope. 2 16 42 08 04 0 04 08 42 16 2 
Curves for the dependence of |T'| on the relative Figure 5. Dependence of the 
detuning y = Aw/Awo at values of Tp = 0, 0.1 and 0.3 modulus of the reflection co- 


are given in Figure 5. efficient on detuning at various 
As is seen from Figure 5, at a given value of de- values of Tp. 

tuning y the change in value of IT | and, consequently, 

in the amplitude of the reflected wave will be maximum with [Tp | = 0. This means that the 

value of noise caused by frequency instability of the generator at Ip = 0 will be maximum. If 

the dispersion of frequency instability is small in comparison with the bandwidth of the resonato1 
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then an increase in Tp | to 0.1 - 0.2 substantially decreases the frequency noise. 
5. MATCHING WITH || CLOSE TO ZERO 


In the observation of the absorption signal in a superheterodyne circuit the value of the 
reflection coefficient Tp is chosen close to zero. 

Let us examine the modulus of the reflection coefficient as a function of detuning y at 
values of Ip lying within the range from zero to 0.3. 

We shall limit the discussion to cases of small detuning: 4 


=| ~ |< 0,2. (26) 
In this case 
Ges dora Comme | 
and Equation (25) may be written in the from | 
Ty =bVe+y, (28) 
where ae oF 
b= - Dey Os eae : 


Curves for the dependence of |I| on detuning under various matching conditions are given in 
Figure 6. 

It is seen from the curves that at the same value of|I'| and the same relative detuning y 
the value of the modulus of the reflection coefficient at above-critical coupling is smaller than 
at below-critical coupling. 


Figure 6. Behavior of the modulus of the re- 
flection coefficient at small values of detuning 
y and at above-critical and below-critical coupling. 


1 - Cirtical matching (I'p=0); 2 - below-critical 
matching (Ip=0.1); 3 - above-critical matching 
(Tp=0. 1). 


Let ['p represent the reflection coefficient at below-critical coupling, Ta the reflection 
coefficient at above-critical coupling, and 


2 2\T 


fe Pele As Teal ah 


= (29) 


The ratio of the modulus of the reflection coefficient at below-critical coupling to the modulus 
of the reflection coefficient at above-critical coupling may be written in the form 


Meee V Deas 


(30) 


Since, in accordance with Equations (29), a, > ap, then from Equation (30) it follows that 


|Tg|>|Tal- (31) 
It can be shown that for y 4 0 


d| Tp | Vat y 


ea re (32) 
dy 
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from which it is seen that 


a|\T rea We 
= | = Ie (33) 


Inequality (33) indicates that at above-critical coupling frequency instability of the generator 
has less effect than at below-critical coupling. 

Substituting into Equation (32) the values of a, and ap at various I, and various detuning 
y, we may obtain a quantitative determination of the decrease in slope and, consequently, the 
decrease in the influence of frequency noise in the transition from below-critical matching to 
above-critical matching. The expression for the slope of the reflection curve d| T'| /dy at small 
detuning values, utilizing Equation (28), has the form 


aq{|V\_ janie hie 
“ay + yaTe 4 


It follows from Equation (34) that with a decrease in |T| to zero the value of d| T| /dy ap- 
proaches +1/2. With |I)| equal to zero, the value of d |T| /dy = +1/2 for all values of y 
satisfying Equation (22). With Ty = 0.1 and y = 0.05 the value of d |T| /dy = 0.15. With 
Tp = 0.3 and y = 0.05 the value or. |dI'| /dy = 0.075. These examples show that at small 
detuning (y = 0.05) the slope of the reflection curve depends substantially on | Tp | . he 
curve for the dependence of the slope of the reflection curve d|T| /dy on the value of | Tp! 
at a relative detuning of y = 0.02 for below-critical coupling is shown in Figure 7. With 


| KH! < 0.2 for above-critical coupling the derivative of the reflection curve is practically 
the same as the curve in Figure 7. 


ail 
ay 
3 
04 
93 ==. = ae 
Figure 7. Dependence of slope of 
92 reflection curve d |T | /dy on the 
value of lg |I',| with relative 
Bf detuning y = 0.02. 
lgIT! 
T° 10” ig? 1g" f 


6. INFLUENCE OF THE DISPERSION SIGNAL ON THE SHAPE 
OF THE ABSORPTION SIGNAL 


In observing the absorption signal it is desirable that the change in amplitude of the re- 
flected wave, caused by paramagnetic absorption, be maximum and that the influence of para- 
magnetic dispersion be minimum. It was previously shown that if the chosen value of detuning 
y = 0, then the value of the absorption signal AU, will be maximum under the condition Ip = 0. 
The equations 


T,=0, y=0, (35) 


defining the optimum coupling and tuning in observing the absorption signal, were obtained with- 
out consideration of the influence of paramagnetic dispersion. The influence of dispersion is 
more evident at detuning y £ 0. It is of interest to examine the influence of dispersion on the 
shape of the absorption signal at detuning y = 0. 

Let us discuss this problem in the case of a reflex spectroscope without automatic frequency 
control. Let us assume that in the absence of paramagnetic resonance y = 0. Paramagnetic 
dispersion causes a change in the energy of the magnetic field within the resonator, which is 
equivalent to a change in the resonator inductance L by the value AL. 

Assume that paramagnetic absorption causes an increase in equivalent loss resistance 
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by a value of AL. This results in changes in Ip and y: 


Pa ets Oe 
AT) = (Rye A” 


AL i 
Ay = Qy 5 = SB AL. (36) 


The result is that the value of the reflection coefficient [ is changed. Taking into account 
Equations (25), (35), (36) we obtain for AT the following approximate expression: 


AT ~ + V (2aT 5)? + Ay? = = V Ar? + (@ AL). (37) 


It follows from Equation (37) that the shape of the reflected signal under condition (35) is a 
complex function of the absorption and dispersion signals. 

In a spectroscope with a reflex resonator and without automatic frequency control of the 
operating cavity we shall compare the contribution of dispersion to the observed paramagnetic 
resonance signal under condition (35) with the contribution of the absorption signal. We shall 
show that with an increase in IT | at tuning of y = 0 the influence of the dispersion signal is 
substantially decreased. Let us examine the matching at which 


IT >| |Tninl, y= 0, : (38) 
wherein 
iMmin| >| Ay], |Uminl >| AT, |. (39) 


Since the magnitudes of |Ayl and |ATp| areusually quite small, the value of | Tyin! may be 
considerably less than 0.01. In the case defined by Equation (38). 


4 
Disregarding quantities of second order in Equation (40), we find that 
AT ~ AT). 


Thus, with sufficiently large | I,| the influence of the dispersion signal on the shape of the 
absorption signal may be disregarded. 


7. MATCHING IN SQUARE-LAW DETECTION IN THE PRESENCE 
OR FREQUENCY INSTABILITY OF THE GENERATOR 


It has been shown that in square-law detection the signal magnitude after the detector will 
be a maximum when the reflection coefficient Ty =41/V3. If it is not possible to disregard 
the influence of frequency noise within the spectroscope, then the two cases of matching are 
not equivalent. 

It can be shown that the slope of the reflection curve with Tp = 1/V3 and y < 0.1 is defined 
by the approximate expression 


The slope of the reflection curve at below-critical matching is approximately 14 times greater 
than at above-critical matching. Usually the ratio of the dispersion of frequency instability 

to the one-half bandwidth of the resonator is much less than unity. Hence we may approximately 
consider that the dispersion of amplitude instability at the detector input is proportional to 

the slope of the reflection curve. This indicates that in the case of square-law detection at 
above-critical matching the frequency noise will be 14 times less than at below-critical matching. 


226 


CONCLUSION 


In a spectroscope with a reflex resonator a decrease in IT | leads to an increase in the 
slope of the reflection curve at small values of detuning y. 

The increase in slope results in an increase in the effect of the dispersion signal on the 
shape of the absorption signal and there is also an increase in noise caused by frequency in- 
stability of the microwave oscillator. 

With IT] = 0 and precise tuning at y = 0the shape of the electron spin resonance signal is 
a complex function of the absorption and dispersion signal. In this case the noise caused by 
frequency instability is maximum. 

The increase in frequency noise and the influence of the dispersion signal with a decrease 
in ITpl explains why superheterodyne radio spectroscopes do not need i-f amplifiers with ex- 
tremely high gain, for an increase in i-f gain in ordinary superheterodyne spectroscopes implie: 
a corresponding decrease in |Ip|. 

With an increase in incident power at the resonator, in order to maintain the required 
operating conditions of the frequency converter in an ordinary superheterodyne spectroscope it 
is necessary to decrease iTp |. This leads to an additional increase in the influence of frequency 
instability. 

We wish to thank A.A. Manenkov for his valuable advice. 
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LARGE-SIGNAL OPERATION 
OF BACKWARD-WAVE AMPLIFIER 


é 


G.N. Rapoport 


On the basis of the numerical solution of nonlinear equations this paper investigates the 
amplification of signals of finite amplitude by a single-stage O-type BWT. It is shown that 
nonlinear effects in the BWT are much greater than in an ordinary TWT and practically 
exclude application of small-signal theory to the analysis of a BWT power amplifier. The 
paper also presents a brief discussion of the operating conditions in which self-oscillations 
are locked by a periodic signal at the BWT input. 


SINGLE-STAGE AMPLIFIER 


Operation of the O-type BWT as an amplifier was first reported by R. Kompfner [1]. 
The small-signal theory of BWT amplifiers has been elaborated in References 2 and 3. It 
has been established that up to the self-excitation threshold the BWT is a regenerative 
amplifier with relatively narrow bandwidth, which may be shifted to wider limits by varying 
the accelerating voltage. The small-signal power gain in the middle of the G band is determinec 
from the approximate equation 
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where Ig is the electron beam current and I, is the critical current corresponding to the 
self-excitation threshold. In approaching the self-excitation threshold the gain rises, 
the bandwidth Af decreases and the produce AfV G remains almost constant. 

The present paper discusses the calculation of the gain and efficiency of a backward- 
wave amplifier taking into consideration nonlinear effects and establishes the limits of ap- 
plicability of linear theory (small-signal theory). 

As is known, for conventional TWT's substantial departures from linear theory occur only 
at output powers greater than 0.4 - 0.6 of saturation power. Hence small-signal theory is 
applicable in practical TWT calculations. 

In conjunction with the regenerative nature of the amplification in a BWT it is to be ex- 
pected that nonlinear effects will be of greater significance, particularly as beam current 
approaches the critical value. 

The calculations presented below are based on an approximate model and the system of 
BWT equations discussed in References 4 and 5. For the sake of simplification the problem 
of the influence of the space charge was not considered. In this case the system of BWT 
equations has the form 


dF — I we = 
ee = iz(a, x) 
aS = Ie da, (2) 


2 (a, 2) Az (a, x)\* lek) BE 
(1-8) noun Gy, a 


az 


where F(x) is the dimensionless amplitude of the high-frequency field, wherein 


F (a) = — F (a); SS = ory (4) 


P is the energy flux of the field; C is the gain parameter; Ip is the beam current; V is the 
accelerating voltage; v is the initial electron velocity; u =u (w) is the backward-wave phase 
velocity in a "cold" retarding system at frequency w; 


ss v @ 
GS eG Ge 
uU v 


s is the electron path; 6 is the loss parameter; Z(@,x) is the electron phase relative to the 
wave in a cold system at point X; @ = -Z(a@, O) is the initial electron phase. 
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Figure 1. Distribution function of field Figure 2. Distribution function of field 
amplitudes for various amplitudes at amplitudes for various amplitudes at 
tube output. Parameter b = 1.385. tube output. Parameter b = 1.465. 
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Equation (2) describes the excitation of the retarding system by the electron stream and 
Equation (3) describes the movement of electrons within the high-frequency field. If the wave 
dispersion u(w) in the retarding system is known, then the frequency w is uniquely associated 
with parameter 


which, together with the amplitude of the field at the tube output F(0), defines the initial 
conditions for the system of Equations (2)-(3). 

With the frequency and amplitude of the wave at the output F(0) so given, the amplitude 
of the wave at any point s > 0 is determined by the solution of the system of Equations (2)-(3). 

For tubes with length of interaction space L the amplitude of the signal at the input will 
be F(xy), where x; = a Ie 

The system of integro-differential equations in Equations (2)-(3) leads approximately to 
a system of 34 ordinary (nonlinear) differential equations of first order by replacing the 
continuous initial-phase distribution of electrons with a discrete distribution for 16 equally 
spaced initial phases. Integration was performed on the "Ural"! electronic computer by the 
Runge-Kutta method of the fourth order. 

The following values of the parameters were adopted in the calculations: C = 0.05, 6=0.1 
("'cold" losses are approximately 0.27 db per length of retarded wave). 

Figures 1 and 2 show the distribution functions of field amplitudes for various amplitudes 
of F (0) at the output, the values of which are indicated on the curves. For Figure 1 the 
parameter b = 1.385, which corresponds to a signal frequency equal to the self-oscillation 
frequency of a BWT the length of which is brought to the threshold of self-excitation. For 
Figure 2, b=1.465, which corresponds to a lower signal frequency equal to the natural 
frequency of self-oscillation of the BWT with a current ratioI./I, = 1.15. 

Values of x are plotted on the lower abscissa scale and the ratio of operating current to 
starting current ([9/I,) is plotted on the upper abscissa scale. Similar calculations were 
performed for a number of values of the parameter b from 1 to 2. In the region x < 2.055, 
where I,p/Ix < 1, the resulting solutions of F(x) describe the BWT amplifier conditions. With 
I)/I, > 1, where x >2.055, the resulting solutions of F(x) correspond formally to conditions 
wherein self-oscillation is locked by an external signal. 

At small output amplitudes (less than 0.3) the functions F approximate the functions 
calculated from small-signal theory. A characteristic nonlinear effect is the dependence of the 
position and value of the minimum of F(x) on the amplitude of F(0) at the output. The envelope 
of minima of F(x) is represented by the broker#line in Figures 1 and 2. 

Let us discuss the principal characteristics of amplifier operation. Figures 3 and 4 show 
the dependence of gain G (in decibels) on the parameter b for various excitation powers Pjy. 


0 ai 4) 
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Figure 3. Gain Gas a function of parameter Figure 4. Gain G as a function of parameter 


b for various excitation powers Pjn; 16/li- = b for various excitation powers P;,; nee = 
0.8; numbers on curves are values of 1.0; numbers on curves are values of 
Pjn/ClIoV (db). Pip/CIoV (db). 


From Figures 3 and 4 we can calculate the dependence of gain on the signal frequency f 
at a fixed accelerating voltage V = const or the dependence of gain on V with f = const. Let 
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f; and fg represent frequencies corresponding to values of by and by with V = const. Then, 
as is easily seen, 


(mas = = ; (by b,), (5) 
ftp 
gr 


with f = const 
a1 = 2 (by — by), (8) 


where fT is the mean frequency, V is the mean accelerating voltage, Uph and Upp are the phase 
and group velocities of the retarded (backward) wave. 

It is seen from Figures 3 and 4 that with an increase in excitation power the maximum 
gain of the BWT decreases, the maximum of the gain curve is displaced in the direction of 
lower frequencies and the bandwidth increases. From Figures 3 and 4 and from Equation (5) 
it is seen that,with the given values of I,/I, and relative excitation power P, ,/CIgV>in order 
to increase the bandwidth it is desirable to decrease the dispersion of the system and to increa: 
the coupling of the system with the electron beam as well as the value of the beam current I). 

For signals with frequencies above or equal to the BWT self-oscillation frequency fy at the 
self-excitation threshold the gain decreases with an increase in signal power; however, at 
frequencies below f, a nonlinear increase in gain is possible with an increase in excitation 
power. An explanation of a similar effect in the TWT with large values of the parameter b 
is given in Reference 6. 

Figure 5 shows the gain G (in decibels) as a function of the normalized excitation power 
(also expressed in decibels) for signal frequencies corresponding to maximum gain. For 
these same frequencies Figure 6 shows the output power Pot Clan as a function of the 
excitation power. The broken curves in Figure 6 represent the difference (Pate in Clg: 
For comparison in Figure 6 the arrow indicates the maximum BWT power with the given 
parameters C, 6, Ig, V in oscillator operation (the tube provides such power if the length 
of the retarding system is 40% greater than critical). 

FromFigures 5 and 6 it is seen that with 
I)/Ik = 0.8 the region of linear gain is limited 
to extremely small powers, not exceeding 10-1! 
of the saturation power (efficiency less than 1% 
In approaching the self-excitation threshold the 
"linear' region is further diminished. It follow 
from this that small-signal theory is practicall: 
inapplicable for BWT power amplifier calculati 
With I,/I, = 1 small-signal theory leads to in- 
finite gain, whereas even at relatively low ex- 
citation power Pj, = 0.01 CIV the actual gain 
is less than 17 db. M.S. Neyman has discussec 


WP We le a Fin , db, the possibility of the synchronization of an O-ty 
Iv backward-wave tube by an external periodic 
Figure 5. Gain G as a function of signal. In Reference 7 he examines methods 
normalized excitation power for of synchronization and presents a qualitatave 
signal frequencies corresponding characteristic of expectable dependences of the 
to maximum gain. I0/I, = 0.8 locking bandwidth on tube parameters and sign 
and I)/Ix = 1.0. power. 


Turning to an examination of locking condi 
from the data of the calculations performed, we note that the subject of the stability of these 
conditions is not strictly within the province of this paper, for we proceed from assumptions c 
the stationary nature and single frequency nature of the oscillations. Nevertheless, with the 
information at hand we can easily confirm the existence of stable locking conditions. 

The heavy line in Figure 2 represents the field distribution function F(x) for the self- 
oscillating state of a tube with length Xy, = 2.15 corresponding to Lidk =1.15. In this state 


the frequency of self-oscillation fy =f, corresponds to b = 1.465 and the amplitude at the 
output Fy (0) = 0.927. It is evident that a shorter tube, for example, xy, = 2.12 d./I, = 1.08) 
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for which the frequency of self-oscillation 
fo =f2 > f,*, may be synchronized by a signal 
of frequency f1 if the amplitude of the signal 
corresponds to the distribution function of the 125 
self-oscillating state and is Fj (x{,) = 0. 06 (the 
point with coordinates x} and F = Fy (x{) in 10 
Figure 2, indicated by a dot). Since the self- 
oscillating state of a tube of length Xz, is stable, 
the locking conditions of a tube of shorter length 
x} are also stable, which conditions are char- 
acterized by the same field distribution and 05 
current distribution as the self-oscillating state 
(heavy line in Figure 2). 0,25 
The conditions of synchronization of a tube 
by signals of frequency f, but with large input 


amplitudes must also be stable. These condi- Ls 005 Of = OG 2 OS a db 

tions are represented by the descending branches 0 

of the distribution function which are located Figure 6. Output power P9,;4/CIoV as a 

above the line of the self-oscillating condition function of excitation power for signal 

F, and correspond to large output amplitudes. frequencies corresponding to maximum 
From Figure 2 is seen that the lines of gain. Broken curves in Figure 6 cor- 

stable conditions intersect the ascending respond to difference (Pout-Pin)/CIoV. 


branches of the distribution function with output 
amplitudes F(0) = 0.278-0.741 smaller than 
for the self-oscillating state at the same frequency. Consequently, the ascending branches of 
the distribution function over the examined segment 2.055 < xy < 2.15 correspond to unstable 
conditions. 
If in conjunction with the above discussion we assume that within the limits of stability 
of the locking conditions the BWT output amplitude must not be less than the amplitude of free 
self-oscillations F,, then we may show that the change in parameter b within the locking 
F 
bandwidth is ~ 0.8 oa where F, is the amplitude of the synchrenizing signal at the BWT 
a 

input. Considering Equation (5), for a BWT with the above parameters we obtain the fol- 
lowing evaluation of the width of the region Af of frequency locking: 

a > F . (7) 


: Ugr 


It also follows from the above remarks that the self-oscillations of a BWT with finite 
output amplitude may be locked by as small a signal as is desired at the tube input if the 
frequency of the signal is close to the frequency of the self-oscillations. 

The locking bandwidth of the tube, all other conditions being equal, will increase in inverse 
proportion to the dispersion of the retarding system. Let it be noted that the results obtained 
are in qualitative agreement with the conclusions presented in Reference 7. 

A further problem in the theory is the examination of bi-frequency (asynchronous) con- 
ditions of a BWT and the determination of the limits of stability of synchronous operation. 
However, the results obtained also indicate, in principle, the possibility of use of a syn- 
chronized tube as a narrowband "amplifier" of modulated signals with a correspondingly 
limited spectrum. 


* With an increase in tube length or beam current the frequency of self-oscillation 
decreases [4]. 
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AMPLIFICATION OF SPACE-CHARGE WAVES 
IN THE PASSAGE OF ELECTRON BEAMS 
WITHIN MEDIA WITH INDUCTIVE SUSCEPTANCE 


L.V. Kasatkin 


On the basis of analysis of a one-dimensional system, this paper discusses the mechanism 
of interaction leading to the amplification of space-charge waves in the presence of inductive - 
resistive reaction of the medium beyond the natural passband. An approximate calculation 
of the nonlinearity of the bunching process is presented which permits evaluation of the gain 
characteristics. 

The question of the amplification of space-charge waves in the passage of electron beams 
through a system with inductive susceptances or conductances relative to the beam has been 
discussed in a number of works. 

In Reference 1 there is described a retarding system of the ''double-comb" type and as 
the result of solution of the field equations it is shown that beyond the passband of this system, 
when the susceptance referred to the beam is inductive, the occurrence of a rising wave in 
the electron beam is possible. In Reference 2 the amplification characteristics of space- 
charge waves are determined (by a linear approximation) for different types of susceptance 
of the drift walls for the cases of cylindrical and plane symmetry. An extremely interesting 
application of this principle of amplification is described in Reference 3 where the so-called 
"inductive drift" (achieved in the form of a retarding system whose passband lies within a 
frequency region above that of the operating region) served as the active amplifying decoupler 
ina "hybrid" TWT. By using the amplification of the waves within the beam during its in- 
ductive drift passage it was possible to shorten the sections of the retarding waveguides of 
the buncher and of the collector and thereby to insure high stability of the amplifier at con- 
siderable gain and bandwidth. 

The amplification mechanism of space-charge waves is not difficult to trace in the ex- 
ainination of a one-dimensional model in which the electron beam freely passes along thin 
channels in a medium having an inductive susceptance y = -jBy, = -j/wLy and also, in the 
general case, a conductance 0. The parameters of the medium are so chosen that in the 


absence of an electron beam no electromagnetic waves may be propagated in the medium. 

In the passage of an electron beam with variable components of current density ig and 
charge density 0g through the system discussed there arises an induced charge with density 
Pc and current with density 


; 1 
io = OE, + 7-\ Eudt, (1) 


Here the coefficient of induction Lj characterizes the specific inductive reactance of the 
medium. During harmonic excitation the coefficient L, and the negative dielectric constant of 
the medium €¢ (€¢ < 0) are related by the relationship L; = -1/w2e,. The total charge density 
in the system 


Oe Paina: (2) 


and the total current density 


. 5 5 1 1 LE, 1 
i, =rot H ig + cE +e TG t \ Edt. (3) 


In determining the charge induced in the medium we shall use the equation 


div he — Px 
: (4) 


and the discontinuity equation of beam convection current 


ie dp 
div tg Las = (9) (5) 
Applying the divergence operator to Equation (3), taking into account Equations (2), (4), 
(5) and the time function eJ“t, we obtain 


p= — mq (1+ ——-|*. (6) 


7@&  ~—«WE 


Depending on parameters o, €, L there arise in the medium induced charges, differing in 
phase and magnitude, the electric field of which interacts with the electron beam. The 
nature of this interaction may be determined by proceeding from the wave equation for the 


one dimensional system discussed: a 
PE, OE, 1 ‘ 0i 
a2 Wire Be Sapa lverye © (7) 


Here i is the current density, equal to the sum of the densities of the convection current 
within the beam and the current in the medium. Expressing Pg by Pg in accordance with 
Equation (6), expressing Pg and i by Ez, (as was done in References 2 and 4) and assuming 
that all variable quantities depend on z and t according to the expression ej(“t-$z), we obtain 
instead of Equation (7) 


r -(€ / 8) Br 
E, (8?— ope + jous 4 wi : oe B, ' ==i(0) (8) 
(8. — 8) (1 ] WS we 


where By = p/ug = (M0/€ ug)1/2 is the wave number in the plasma, € is the dielectric 
constant in vacuum, Ig is the constant component of beam current density, ug is the con- 
stant velocity of the electrons, and Be = w/ug. 

Equation (8) is a fourth-degree equation in terms of f and thus defines four possible 
x waves: two decaying waves in the medium, not associated with the beam, and two space- 
charge waves. The propagation factor of the space-charge waves in the beam is defined by 


the expression 
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B Sie a a ‘9 
a (9) 
where 
ar use! | eres (10) 
Sy ee we oe |? x >’ 


€4 is the equivalent complex permittivity which includes the conductance and susceptance of 
the medium. 
Substituting Equation (10 in Equation (9) it is seen that in the general case, with inductive- 
resistive reaction on the part of the medium, amplification of the space-charge waves in the 
beam is possible, that is, 
B= Ge tBo(p+ ia), (11) 


where 


(12) 


alsa GT 20s y 
/ ae 6 \2Th : 
[(O=se) +e) | 


The maximum values of p and q occur at magnitudes of By/ we defined by the expressions 


(13) 


By, 3 / we 
( we i A ae Vou i) 
By, 5 / ws 
ah = {+ —_., 15 
(= aux 1 73 i 
In which case 
081 


Pmax > max ~ (e/aey (16) 
Figures 1 and w show families of curves for p =f (By /€) and q =f (By /w€) with various 
values of the parameter o/we. In the case of small losses in the medium there occurs con- 
siderable amplification of the space-charge waves in the beam with By /we slightly larger 
than unity. With 0< B,/we< 1 significant amplification is possible only with considerable 
losses in the medium, while with o/we = 0 in this region of values of By/we there is no 
amplification. An increase in the losses in the medium decreases the maximum value of 
amplification (with o/we ~0 q>~), but the region of the parameters By / we in which amplifica- 
tion is possible is enlarged. The parameter p (determining the change in the characteristic 
plasma frequency) with o/we = 0 lies within the limits 1 < p< », when0<B,/we< 1 and 
p = 0 with By /we > 1. Thus, in the absence of losses in the medium, amplification of space- 
charge waves in the beam is possible for By /we > 1 with an undisturbed velocity of the 
amplified wave Vy, =Up. Losses lead to a decrease in py,4y and an increase in the region 
of frequencies within which p is other than zero. 
For a physical picture of the interaction in the presence of inductive reaction of the 
medium let us examine the special case of the absence of loss (a/we = 0). In this case the 
induced charge in the medium is defined, in accordance with Equation (6), by the expression 


3%; | 1 fy 
Pe Pd Fi jas | : (17) 


Figure 3 shows the variation of po/pg =f(Br/we). With By /we <1 there is induced 
in the inductive medium a charge of the same sign as in the beam, and with By /we 1! 
Po >. In addition, the electric field intensity EG resulting from the induced charge is 
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Figure 1 Parameter p (change in plasma frequency) as a func- 
tion of the characteristics of the medium B;/we and 0/Bg (Be=We) 


added in phase with the space-charge intensity E, in the beam, with the result that there is 
an increase in the longitudinal repulsion of electrons and amplification of the space-charge 
waves proves impossible. With By/ We >1 the sign of the charge induced in the medium 
changes. The space-charge forces in the beam and the forces created by the induced charge 
EY are opposite in phase (see Figure 4). With |,|>|pgl the resulting intensity of the 
electric field acting in the beam leads to the aeteleration of retarded electrons and to the 
retardation of accelerated electrons relative to the center of the bunch, resulting in an 
improvement of electron bunching in the beam and ensuring amplification of the space-charge 
waves. Upon fulfillment of the condition By /wge = 1 resonance exists in the medium: the 
displacement currents and the induction currents in the medium become equal in value. In 
this case the charge induced in the medium approaches ~ and changes sign in the vicinity 

of W = Wo. 

The a rapiifieation of space-charge waves in the passage of the beam through a system 
with inductive reaction is determined by the improvement of bunching of the electrons in the 
beam along the transit path and, as a result, an increase in the induced current in the medium. 
Power is not transmitted along the system. In this connection, in order to excite waves in 
the beam it is necessary to use a buncher (of the TWT delay-section type or a resonator), 
causing, in the general case, velocity and density modulation of the beam. For amplification 
in the inductive drift a device similar to the buncher is used, wherein the amplified variable 
component of current in the beam excites electromagnetic oscillations. 

Let us conduct an approximate examination of the nonlinear effect of bunching in the 
inductive drift. In connection with the amplification of the variable components of current 
and velocity in the beam nonlinear effects are more intense here than in a drift with ideally 
conducting walls. 

If the amplitudes of the variable components of current and electron velocity at the in- 
put into the inductive drift are ijg and uyo respectively, then the variable component of 
velocity z of an electron flying into the drift at moment tg is defined by the expression 
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Figure 2. Gain parameter q as a function of the characteristics 
of the medium B; /we and 0/B, (Bg = €). 


‘ P oq 
Wy (5) = Wy Chip gs Sin Of, -F Thy ee sh /',, G3 COs fy. (18) 
0 
Electron transit time within the drift (z =s) is 
; ' : ds 
LOL Ugg, é 19 
) Ug 4 Uy Ch I, G2 SiN Oly =F M9 7 P sh ky, gq COs Oly ( ) 
U 0 
Upon satisfying the conditions 
rics Ray ahah ee ity oy ei yo o 
ae <1, FH shhpgs<! (20) 
we obtain from Equation (19) 
0 = o(t, —-t,) ~ ot —X sin (wt, — yp), (21) 


where T = s/uo, 


Danan 
ie iS aa) 


W490 4 (22) 


sh k, qs 
Uy On 1 


y = arcle 


i ia@iae aay None (aay Si 12) */2 
\ rr, sh Ip qs) 2] BG hhy qs — 1)| (23) 


where X is the bunching parameter. 
As is seen from Equation (22), the bunching parameter for electrons within the inductive 
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Figure 4. Characteristics of the longitudinal 
Figure 3. Charge 0, induced in the medium _ distribution of charge density within the beam 


as a function of the characteristics By /we Pq, space-charge field strength Eg and induced- 
of the medium in the absence of losses charge field strength E' (with B;/we < 1) and 
(a/we = 0). El! (with By /we > h, 


drift increases rapidly with an increase in the transit length s (almost exponentially for 
large values of kyqs). The amplitude of the first harmonic of the current within the drift tube 
at a distance s from the buncher is determined from the expression 


ae 3 
ies siege 


» (X) + 21, cos spl, (X), (24) 


Jo(X) and J;(X) are the zeroth-order and first-order Bessel functions. 

The curves in Figure 5 show itm/2Io = f(s) calculated from Equations (22) and (23) at 
various values of Wyq/. In these calculations it was assumed that u,0/ug = 0.01 and 
ij0/Io = 0.05. The region of application of these ratios is limited by the extremely rigid 
conditions (20). The bunching characteristics at large values of k.qz may easily be determined 
in the special case where ijg = 0 and the excitation within the drift is achieved by velocity 
modulation of the electrons. In this case the electron transit angle in the drift 


dz 
w (ty— i.) = 0\ (25) 


Ug > Up ch keg, q2 SiN Oly 


and with m = uj9/ug « 1, which is always attainable in practice, 


A 1--+mchk,, qs sin wlo 
© = ot ar In | ! (26) 


®, 1 + msin Of 


The curves in Figure 6 show wt, =f (wto), determined from Equation (26) with m ch 
kpgs = 0.05, 0.1 and 0.15 and Wygq/w = 0.72. The broken curves show wt, = f(wtg) as cal- 
culated from Equation (21) for m ch kpqs = 0.15. Comparison of these curves confirms the 
possibility of using Equations (2 1) -(23) for approximate evaluation of the first harmonic of 
the current. 

From an examination of the curves for ij, = f(s) it follows that, all other conditions 
being equal, the nature of the frequency dependence of the gain K = iy (f, s)/ij9 may be 
varied by the choice of the length of drift with inductive reaction. Thus, with a short length 
of drift, when large values of kyq correspond to the optimum bunching condition, the gain 
decreases rapidly with an increase in By /we, that is, with a decrease in frequency relative 
to resonance (Ww =wg). At such a drift length, when the bunching optimum is attained at a 
frequency 1 < wo, the gain characteristic will have its maximum at a frequency close to 
w . In the region of frequencies wg >w>w 1, due to overbunching (X>Xom) amplification 
within the inductive drift will decrease. All other conditions being equal, with a constant 
length of inductive drift the nature of the frequency dependence of the gain is determined 
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Figure 5. Change in the first harmonic 

of beam current along a drift with induc- 

tive reaction at various linear gains kpq 
(in nepers). 


nla 


Figure 6. Dependence of the electron 
arrival phase wt; on the departure phase 
wtg with various values of m ch (kpqs). af 


by ui10(f) and i10(f), that is the choice of operating conditions and the characteristics of 

the exciter section. For example, if ujg and ij0 increase with a decrease in frequency, the 
we may obtain sufficiently small variations in gain within the inductive drift over a wide 
range of frequencies. 
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CALCULATION OF A MAGNETIC PERIODIC 
FOCUSING DEVICE FOR TWT 


A.L. Igritskiy 


This paper discusses a method of calculating a magnetic system for periodic focusing 
of an electron beam in a TWT which permits consideration of the influence of the material 
and the dimensions of the pole pieces as well as of a finite number of magnets on the in- 
duction of the resulting periodic field. The method also permits consideration of the influence 
of disturbances of the periodicity and the ends of the device on the induction. 

The calculation was carried out by the method of successive approximations. The 
method described ensures better agreement between calculated and experimental results. 


INTRODUCTION 


The use of periodic magnetic fields for focusing of the electron beam in traveling-wave 
tubes (TWT) permits a substantial reduction in the size and weight of the focusing device and 
also eliminates the need for external sources of power. The diagram of a device for periodic 
focusing is shown in Figure 1. Permanent 
ring magnets M (individually magnetized in 
the axial direction) are so stacked that the 
facing polarities of adjacent rings are iden- 
tical. Between the permanent magnets there 
are pole pieces S of a material possessing 
a high saturation induction (e.g., Armco 
steel). The pole pieces contribute to field 
symmetry and greater field induction along 
the axis. 


Figure 1. a - Longitudinal cross-section of Ll 

magnetic focusing device: M - magnet; S - 

pole pieces; 1, 2, 3, 4, paths for fluxes v1, SS 

~2, v3 and wy, respectively; b - distribution d, 
of magnetic potential over radius R = Rj. b 


Design calculations of a periodic focusing device are intended to determine the optimum 
shape and dimensions of the principal elements of the magnetic system: the permanent mag- 
nets and pole pieces. A method of calculation of periodic focusing devices was discussed 
in References 1 and 2. 

As investigations have shown, calculations performed by this method yield values for 
induction which are 20-40% larger than experimental values. This difference in values is 
explained chiefly by the fact that the method of calculation described in References 1 and 2 
does not consider the effect of the reluctance of the pole pieces on the value of induction, 
although this effect may be quite significant. 


1. CALCULATION OF THE MAGNETIC SYSTEM 
IN THE FIRST APPROXI MATION 
In calculation of the magnetic system in the first approximation we disregard the reluctance 

of the pole pieces. It is seen from Figure 1 that the pole pieces form a series of cylinders of 
identical diameter wherein the magnetic potentials of adjacent cylinders are opposite and 

equal to +g while in the gap between them at radius r RL the potential varies linearly 

(see Figure 1b). The distribution of magnetic potential in the inner region between pole 

pieces in this case may be determined by solution of the Laplace equation; it has the form [1] 


- 205 ea ile { $ nq) sin (n5,) sin (= z) (1) 
ie COS ‘ BS) #)} 4 
@M(z;r) = > CLE: (1—co 1 te) 
n=1 Jet ia R,) 
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d = 
y=7. (2) 


In Equations (1) and (2) Lis the period of the magnetic field, Ry is the inner radius of the 


pole. pieces, d 1 is the distance between pole pieces, and Joe r) is a zeroth-order Bessel 


function of the first kind with imaginary argument. ; ; 
The component of the magnetic field induction along the z-axis may be defined as the 


gradient of the magnetic potential and, in absolute value is 


OD Be — COs i S: . q 2 
Bees a AS 4D Uo (1 — cos nt) sin (215) JI,(i Bae r eos ( nw 2) . (3) 


902 a 2n nas 
n=1 ES (i ee Rr) z 


The magnetic field along the z-axis (with r = 0) will vary in accordance with 


4M, (1 — cos na) sin (ns51) 2n70 ) 
- cos Zo. , 
Le (i ae 1) NMS, ( I8, (4) 


The total magnetic flux in the plane z = 0 between pole pieces in the inner region 


x 


Vy == | Be 2Qurdr — 20,P3. (5). 


Here Pj is the permeance for flux 41: 


2n% 
LE Ry) sin nto, 2(1 —cos nm) p . (6) 
2nn \ NING, 1 29 


m=1 Jo (i Sas Ry) 


og -iS,{i 


Py =o 
Ib 
Fee) is a first-order Bessel function of first kind with imaginary argument. 
Tne leakage flux existing beyond the limits of the ring magnets (with R > Rg) also has 
considerable effect on the value of the maximum attainable field along the axis. 
The leakage flux between pole pieces in the outer region is calculated in the same 
manner as the effective flux in the inner region and for it we derive the expression [1] 


P. = 2D,P». (7) 


Here Po is the permeance for the leakage flux: 


Phe oN 
co p7(1) (i a ) 
eS eee) sin (3069) 2 (1 —cosnt) pp (8) 
2 [Py \ Se 2 
n=1 Os Re) a ‘" 
Og == a a (9) 


In Equations (8) and (9) Rg is the outer radius of the magnet, dg is the thickness of a 
single magnet, H(})and H(}) are Hankel functions of the first kind and of the zeroth and first 
orders, respectively. 

In addition to the magnetic fluxes ~1 and v2 in the focusing system there exist the magnetic 
flux #3 between the end faces of the annular pole pieces and the flux w4 in the gap between the 
ring of a pole piece and the inner surface of the magnet. 

The paths of these fluxes are shown in Figure 1 by 3 and 4. In accordance with Figure 1 
the permeances for fluxes wg and #4 may be calculated from 


Pe ae (11) 


Here R3 is the outer radius of the annular pole piece, R4 is the inner radius of the permanent 
ring magnet. 
The total permeance between pole pieces 
4 
Py = 5) Pi. (12) 


i=1 


If the difference of magnetic potential between pole pieces is 2g, then each magnet 
will create a magnetic flux 


ip ee Pe (13) 


The induction in the permanent magnet is 


2D,Po 


B ee EATS Re TTENS 
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(14) 


Let us assume that the magnetic field strength is constant over the length of the magnet. 
Then from the definition of the field strength as the gradient of the potential we can write 


Had = 2p, (15) 


where Hy is the field strength in the permanent magnet. 
From Equations (14) and (15) we obtain the relationship between induction and field 
strength in the permanent magnet: 


Be = dzPo 
Hy, 1 (R; — Ri) ‘ (16) 


On the other hand, By and Hy are related by the demagnetization curve B =f(H). Fig- 
ure 2 shows the demagnetization curve for anisotropic ferroxdur. If all the dimensions of the 
magnetic system are given, we can calculate the 
right hand side of Equation (16). Extending a straight 
line from the coordinate origin to the curve 


doP, 
o% = arctg — (Re — BB) : (17) 


we find the operating point of the permanent mace 
net (see Figure 2). Knowing the values of Hy and 
Bm at the operating point, we can ascertain how 
close the operating conditions of the magnets are to 
optimum conditions. Equations (15) and (4) permit 
plotting the distribution of the magnetic field along j 
the axis of the periodic focusing device. ofp att APD a 

As the result of Calculation of the magnetic Hf, 08 
system in the first approximation we also deter- 
mine the fluxes #1, v2, w3 and ~4, which will be 
used in calculation of the magnetic circuit in the 
second approximation. 


Figure 2. Demagnetization curve for 
anisotropic ferroxdur B = f(H). 


2. CALCULATION OF THE MAGNETIC SYSTEM 
IN THE SECOND APPROXIMATION 


In a similar calculation of the magnetic circuit let us consider the reluctance of the 
material of the pole pieces. Figure 3 shows the distribution of magnetic fluxes in the pole 
piece of a periodic device. It is not difficult to conclude from the figure that there must 
exist a certain separating radius Rp. The lines of force which proceed from the magnet at 
R2 Rp will form the leakage flux 


tp: = Bra (R, — Rp), (18) 
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and the lines of force proceeding from the magnet 
at R< Ry form the flux 
Pe = P+ Ys + the = Bin (Rp — Ri). (19) 
From Equations (18) and (19) the magnitude of 
the separating radius is 


Ree /% Ro + bee (20) 
he + Dp 


In Equation (20) fluxes Jp and pp are known 
from the calculation in the first approximation. 

At the various portions of the pole piece there 
flow different fluxes; hence let us divide it into 
four portions ba, be, ad andde. First let us 


Figure 3. Distribution of magnetic determine the drop in magnetic potential at the 
fluxes in the pole piece of a periodic individual portions. 
focusing device. Portion be. In this portion there flows the 


leakage flux, the value of which varies'from zero 
at R= Rp to vg at R= Rg. Through the section of radius R, as is seen from Figure 3, there 
flows a flux ~2 minus the flux into the pole piece above the section of radius R: 


Ye = Ye — 8, n(R3— R?). (21) 


The induction in section R of the pole piece 
bs — B,, 1 (R2 — R?) 
pes 2 eee (22) 
onR —P 


Here B,, is the induction in the permanent magnet, which is known from the calculation in 
the first approximation; A» is the width of the pole piece (see Figure 1). 

Using the magnetization curve B = f(H) of the pole piece material*, from the value of 
BR we can find the corresponding value of HR. Given several values of R in portion bc, we 
can calculate for them the values of Br and plot the curve for Hp = fpe(R). 

The magnetic potential drop across the section be equals 

Rz R, 
Ute = \ Hp dR = \ foe(R)aR. (23) 


Rp 12 


This integral can be evaluated by one of the approximate methods, since Hp = fpc(R) is 
given in the form of a table or plot. 
Calculation of the drop in magnetic potential at portions ba and ad is quite similar. 
Portion ba. At portion ba there flows an internal flux, the value of which varies from 
zero at R = Ry to vp =%,+ 9+ 4 at R=Ry. Through the section of radius R there flows 
2 flux ~p minus the flux entering the pole piece from the magnet below the section of radius 


Vp aa vp a By, 1 (Re— 4 . (24) 


The induction in the section of radius R 


Dp SF IB nie te Ri) 


B = 
B aRAp 


(25) 


*If the pole pieces are not annealed after fabrication, in order to take into account the 
possible effects of hardening, itisdesirable to plot the magnetization curve of the material 
from specimens having the same width as the pole pieces. 
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By calculating the value of Br for several values of R and by using the magnetization 
curve we can plot the curve of Hp = fpa(R) for portion ba 


Rp oe 
Ong = \ Hn OR Rh Toa(R) dR 56 
\ aybosle) aR, (26) 


Portion ad. At portion ad there flows an internal flux, the value of which varies from 
pp =%1+ dg + v4 at R = Ry to ¥1 + Yg at R=Rg. Through the section of radius R there 


flows a flux ~p minus the flux branching into the gap between pole pieces along path w4 at 
the portion between R4 and R (see Figure 3): 


vp =, — Bn (R?2— R?). 


(27) 
Here B is the induction in the gap between pole pieces along the path wy: 
2 Onan 
B= wy eee aes 0 An, (28) 
In the absolute system of units uo = 1 and the induction in the gap along path w4 is 
numerically equal to the field strength in the magnet. 
The induction in the pole piece section of radius R 
pg ae Te 
Ba Dp m mt ( 4 ) ; (29) 


TRAD 


By calculating the values of Bp for several values of R and using the magnetization 
curve we can plot Hp = faq(R) for portion ad. The drop in magnetic potential at portion ad 


Rg 
Una = \ faa (R) (dR). (30) 


Rs 


Portion de. At portion de there flows through the pole piece a flux which varies in 
value from #7 + wg at point d to wg at the end face of the pole piece (point e). Through the 
circumferential cross-section at a distance z from the central plane of the gap there pass 
flux {3 and the part of flux ~1 which leaves the inner surface of the pole piece cylinder 
and proceeds to a point with coordinate z: 

be = Ya + \ By (Ra A 2aR, da. (31) 


0 


Here By(R1;2Z) is the radial component of induction at the inner surface of the pole piece 
cylinder. It can be determined as the gradient of magnetic potential from Equation (1); 
in absolute value it will be 


B, (7; 2) = Wo la. 
: tPA Oa 
—=id eee ge 3 
re S 2 ( L '/sinnas, 4D, (1 — cos nat) xa (aun 2) : (32) 
= Bi » (/ Qn NTS, ib \L 
a, Je (3 panes hs) 
iB 


Substituting By into Equation (31) and integrating, we obtain 


Dz = Ws 35 
3 2 aie , 
—iJ eit 
(oe) u i(i iE 1) 


sin nt31 4D, (1 — cos nz) [ ane (208 2)| 33 
Joli PBiRi\al Be n aT. 2) ; 


The induction in the pole piece cross-section with coordinate z is 


Lae ee 34 
B= a (34) 
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Calculating the values of Bz for several values of z and using the magnetization curve 
for the pole piece material, we can plot the curve of Hz = fge(z) for portion de 


U ae = (35) 


wl Pee | 
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Q 
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If the drop in magnetic potential at portion de is considerable, this will lead to a change 
in the distribution of magnetic potential at radius R = Rj (see Figure 1b) and it may prove 
convenient to derive a new expression for ®(z; r) instead of Equation (1), which was 
derived for a trapezoidal potential distribution with R = Rj. 

For an individual magnet of the periodic system we can provide the equivalent conversion 
circuit [3, 4] shown in Figure 4. Here Rj takes into account the air resistance along the path 

of internal flux: 


20, 4 


R, res _ 
1p, Pep Paes (36) 
Re, Ree Ra represents the (magnetic) resistance of 
the path of internal flux in steel. This‘resistance 
we shall relate to Jp and it will be 
Rey Re, 


Figure 4. Equivalent conversion diagram of an 
individual magnet in the periodic focusing device. 


ari Ura ar U ae (37) 


Ro Dp 


Here Upg, Ugg and Uge are the drops in magnetic potential at the corresponding portions, 
the values of which were defined above. 
R2 is the resistance in air to the leakage flux wo: 
Xi 


R= = 5, (38) 


Resistance Re2 represents the resistance to the leakage flux wo in steel: 


Vy. 
=e (39) 


Reg a 
Using the equivalent circuit, let us determine the total resistance of the magnetic 
circuit 


el e2 
Pa ale eee (40) 


Here Rei = Rj + 2Ro1 is the equivalent resistance for the internal magnetic flux and Reg 
is the same for the leakage flux. 
On the basis of Ohm's magnetic law the flux created by the magnet 


2 A yds 
Ab ==0B 7 a0 (Ro Ae Ry i (41) 


From Equation (41) we obtain the ratio between the induction B,, and the field strength 
Hy in the magnet in the second approximation: 


== $ Otay (42) 
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On the other hand, Bm and Hy are related by the demagnetization curve B = f(H) 
shown in Figure 2. On the plot of B = f(H) we extend a straight line through the coordinate 
origin at an angle of 


ds, 


Oe = arc t 
= 8 wR RRs” (43) 


and find the operating point of the permanent magnet in the second approximation. 
The mmf created by the magnet: 


UO, tliat (44) 
The induction of the magnetic field at the axis can be determined in the second approxi- 

mation from Equation (4) if, instead of the difference of potentials between pole pieces 2p, 

we substitute the drop in magnetic potential at the resistance R;, which, as is easily seen, is 


(45) 


If the calculated value of Bzo in the secondapproximation differs but little from the value 
of Bzo9 obtained in the first approximation, then the calculation can be considered concluded. 
If there is significant departure in these values, the calculation must be repeated with the 
new values of mmf and fluxes obtained in the second approximation. 


3. INFLUENCE OF THE ENDS OF THE FOCUSING DEVICE 


In the case of calculating the reluctance of the pole pieces the magnetic circuits of the 
individual permanent magnets prove to be interrelated and for the entire periodic focusing 
system we may present the equivalent conversion diagram shown in Figure 5. Calculations 
from this conversion diagram permit evaluating the mutual influence of the magnets and also 
permit considering the influence of the ends of the focusing system on the distribution of 
the field within it. 


¥! ie Yelm +9 baap V200%2) bya bon- Yoien2) bony Yorn bon Ylow Donst 


” Ondw! aneg Gngoy, 2 G, a, 
Qn Yin) Anes Wns) nd Wins) N+3 292 Wi aq. 9) 2 Witney Yoon) entt 


| 
1p 
Figure 5. Equivalent conversion diagram of the focusing system. 


The number of magnets in the focusing device must be even, for in this case the 
opposed directions of the fields of the individual magnets will result in mutual compensa- 
tion and along the axis of the focusing device 

ae 
\ H,dz =0 


—oo 


without the occurrence of additional undersirable tails of the magnetic field at the ends 
of the device. For the sake of simplicity we shall consider the periodic focusing device 
to be symmetrical with respect to the edges. By virtue of this symmetry we can consider 
that the mmf's of magnets equally distant from the ends are 


U, a Von} U, = UVon—13 eee Uy = Ven—K- ans Un == Ont: 
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The magnetic fluxes are 
Prin) a Pinta) Won) = Worngay? Vocn) = Won tay" 


Thus it suffices to perform the calculations for only half of the system, for example, 
to the right of the plane of symmetry PP. The currents and mmf's to the left are found 
from the conditions of symmetry. Calculation from the circuit shown in Figure 5 is 
conveniently performed by the method of successive approximations. On the first ap- 
proximation let us consider that the mmf's of the permanent magnets and the resistances 
of the flux paths for iron are identical and known. The values of the mmf's and resistances 
are determined from Equations (44), (37) and (39). It is required to determine the mag- 
netic fluxes flowing along the various sections of the magnetic circuit. 

For the lower contour with number n+ 1 and subsequent contours, on the basis of 
Kirchoff's second law [4], we can write the equations 


R 
meni + Whenpay a + Punta) =U, 
Pa 
home ~ + Winsett + Punts) =U, cts (46) 
sori daetsateess ais ata ee ones 
capt 4D npiyhts + Punter —- =U, 


UCT AY EC On OMM Te Chatce Oe "Om creat aires ey “Sahih ON OO 


die 
Pian) —5— a+ vp, conti t Prien) Ria = Ue 


Applying Kirchoff's first law to nodes an+1, ant+1, etc, we obtain 


, Piin-+41) 
REITE Nae fr ee 


; ; 4 
Vingay = Punta) — Punta = Punta — F Pan? 


Rt fs co So Se aE ee ee ae (47) 
Pinte = Panty — Pape = 
1 
= Vinten — Prepay 1? =P gepay Ph Pugepay Punt? 
i can) mm Pren-ta) a Py (on) re ved ca 
4 
= Bren) — Prana Fo — Punts) + Prt) — ZF Yum: 
The last two equations are applicable with k and n even. 
Substituting the fluxes from Equation (47) into Equation (46), we have 
4 R 
> Punta) (Rex + Ry) + Prints) att U, 
( i 4 R 
Panta) (— +R) — x Puntoki + Pints =- =U, 
te | R 7 
Prn+n an + AR) Pring Fo + Panter — (48) 


= (0. 


Bao dab L1) R, ait = Dy (nt k-H4) 


( Ros ! 
Pian) ( 5 R, + Ry ) Prana (Fy + Ry) + 
‘ { 
* + Wi(n—2) (Ry eS Rx) = > Pity (Ry ae Rx) = U. 
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Thus, we have obtained a system of n equations with n unknown fluxes Yin +1) 
Pin 4 2). -- PL Qn): 

The solution of this sytem of equations by the usual method by means of determinants 
is accompanied by considerable difficulty in calculation. However, the solution is extremely 
simplified by the use of the following approximate method. Usually the resistance of the 
pole piece material R,1/2 is considerably lower than the resistance Rj for flux for air and 
the fluxes Yim +1), Yi(a + 2) ete differ little from one another. Hence, in the first equation 
of system (48) we can with sufficient accuracy disregard the second expression. Flux ~1(n + 1) 
in the first approximation will be 


. rig pysi2iih 
Putt = (49) 


Substituting the resulting value of W1(y + 1) into the second equation of the system and 
disregarding the second expression therein, we find Wy (y + 2), etc. Thus, we can find 
all the fluxes successively. If it is desirable to determine these fluxes with greater 
accuracy, the resulting values for the fluxes in the first approximation must be used in 
determining the fluxes in the second approximation. Thus, flux ~ (n + 1) in the second 
approximation is defined by the expression 


R 
cl 
eee 2 


YP napl = Pea ° (50) 


The other fluxes can be determined in a similar fashion. 

We have discussed a method of determining the internal fluxes ~1 of a magnetic system. 
Completely similar equations can be derived for the leakage fluxes Y. For this purpose 
Kirchoff's second law must be applied to the upper contours with numbers n + 1 and following. 
Kirchoff's first law must be applied to nodes by + 1, by + 2, etc. The resulting system of 
equations is solved in the same manner as described above. 

Having determined fluxes v1(n + k) and ~2(n +k), the fluxes in the remaining sections of 
the magnetic circuit can be determined without particular difficulty. The operating points 
and mmf's of individual magnets of the periodic system can be determined. If the values 
of flux and mmf obtained differ but little from the values determined in Section 2, this 
calculation may be considered final. If there is significant difference in the values, another 
approximation will be required. 


CONCLUSION 


This paper has presented a method of calculation of a magnetic system for periodic 
focusing of the electron beam ina TWT. Its chief difference from the hitherto used Chang 
method lies in the consideration of the reluctance of the pole pieces. The method described 
permits the calculation of the influence of the pole piece material and dimensions on the in- 
duction of the resulting periodic field. It permits calculating the influence of disturbances 
of the periodicity and ends of the device on the value of induction and ensures better coin- 
cidence between the results of calculation and experiment. 


REFERENCES 


1. K.K.N. Chang, Optimum design of periodic magnet structures for electron beam focusing, 
RCARR Gyan oops UG nmlen Oo) 
2. F. Sterzer, W.W. Siekanowicz, The design of periodic permanent magnets for focusing 
of electron beams, RCA Rev., 1957, 18, 1, 39. 
M.A. Babikov, Electrical Devices, PartI, GEI, 1951. 
A. Ya. Buylov, Principles of Construction of Electrical Devices, GEI, 1946. 
Roters, Electromagnetic Mechanisms, GEI, 1949. 


ow 


Submitted to the editors 29 February 1960 


247 


SOME EFFECTS ACCOMPANYING DETECTION 
IN A GASEOUS DISCHARGE 


G.D. Lobov and V.I. Yeremeyev 


This paper presents an experimental comparison of the characteristics of the detection 
of microwave oscillations in a gaseous discharge with the characteristics of enhanced glow 
of the same discharge. Both effects are investigated for helium and neon gas-discharge 
tubes at various powers and discharge currents. 


INTRODUCTION 


The action of an electromagnetic microwave field on a gaseous discharge has already 
been described in the literature. It has been shown, for example, that in applying micro- 
waves to a gas-discharge tube both the value of the discharge current [1] and the glow level 
of the gas [2, 3] vary. In the work described here an increase in the excitation glow of 
atoms and a decrease in the recombination glow in the microwave field were investigated 
in connection with the process of detection in a glow discharge. 

EX PERIMENT 


Experimental investigation of the effects of excitation glow and quenching of recombina- 
tion glow was performed with the setup shown in the block diagram in Figure 1. 


Figure 1. Block diagram of measurement setup. 


1 - Gas-Discharge detector; 2 - microwave gen- 
erator; 3 - synchronizer; 4 - power supply for 
gas-discharge detector; 5 - photomultiplier with 
power supply; 6 - wattmeter; 7 - oscillograph. 


Microwave pulses of duration Tp = 10 usec were applied to a gas-discharge detector 
through a waveguide. For the investigation of the influence of microwave power on the 
various sections of discharge, the waveguide was narrowed to a slot at the point of in- 
sertion of the gas-discharge detector. This permitted application of the microwaves to 
narrow sections of the gaseous discharge. The gas-discharge detector was a glass tube 
with plane metal electrodes between which a glow discharge was excited. The tube was 
filled with an inert gas. 

The glow discharge in the detector could arise from both a d-c voltage and the 
application of special ignition pulses of duration Tig = 75 usec. In such ignition of the 
gas-discharge detector microwave pulses could be applied at different moments, including 
during the afterglow of the gas. 

The brightners of the glow discharge and all changes in this glow under the influence 
of the microwave field were detected by means of a photomultiplier through a slit in 
the narrow wall of the waveguide. The voltages from the photomultiplier and the load 
resistor of the gas-discharge detector were applied to an oscillograph. 
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2. RESULTS 


In all cases during the action of sufficient microwave power on the discharge the 
intensity of the glow at the point of action of the pulse was increased. This increase in 
glow differed for neon and helium. With a decrease in microwave power the glow dropped 
to zero and then there appeared a pulse of reverse polarity (that is, a quenching pulse 
caused by a decrease in recombination glow). 

Figure 2 shows oscillograms of the voltage at the photomultiplier output, obtained during 
the action of microwaves on the discharge 
created by d-c voltage in helium. The oscil- 
lograms are shown in order of increasing 
microwave power. At low powers there is 
seen the quenching pulse of the recombination 
glow at the point of action of the microwave 
pulse (Figure 2a). 


f, 

d 
Figure 2. Oscillograms of pulses of ex- 
tinction and increase of glow for a dis- Figure 3. Oscillograms of extinction 
charge by d-c voltage in helium, p = 9 mm pulses and an increased glow pulse in 
Hg, Tp = 20 usec. Photographed at dif- neon afterglow for different microwave 
ferent microwave powers. powers. p= 21 mm Hg; 
a = P/Py, = 0.01; b - P/Py = 0.05; ¢ - Poe EG: Tig = 75 usec. 
P/Py, = 1; d - gas-discharge detector a-P=0w;b- Bem = = 0.08; 
pulse. Cee, 


With an increase in microwave power the trailing edge of the pulse gradually decays. 
This may be explained by the fact that the electrons, having received energy during the 
pulse, retain it for a while after the pulse, which decreases the probability of recombina- 
tion. The greater the applied microwave power, the greater the decrease in recombination. 
With a further increase in microwave power there appears at the site of the quenching pulse 
a pulse of intensified excitation glow accompanied by a distinct overshoot after the pulse. 

In a discharge by d-c voltage no quenching pulse is observed in a neon-filled tube. 

For more distinct observation of the recombination glow, pulse ignition of the gas- 
discharge detector was employed. Figure 3 shows the oscillograms for a neon tube. 
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Figure 4. Dependence of pulse magnitudes 
at the output of the photomultiplier (AUpn) 
and the gas-discharge detector (AUg) on 
power for a gas-discharge tube with he- 
lium; p = 9 mm Hg, T,, = 10 psec, T., = 
p 1g 
75 psec: 


1 - AUpn(P), microwave pulse occurs 

during afterglow; 2 - AUpn(P) microwave 

pulse occurs during the action of the ig- 
nition pulse; 3 - AUg(P). 
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Figure 5. Dependence of pulse mag- 
nitudes at the photomultiplier output 
(AUpn) and the gas-discharge de- 
tector output (AUg) on microwave 
power applied to a neon gas-discharge 
tube; p = 21 mm Hg, Tp = 10 usec, 
Tig = 75 usec: 
1 - AUpn (P) for afterglow; 2 - UApn(P) 
for ignition; 3 - AUg(P). 


Oscillograms for helium tubes were similar and hence are not presented here. Figures 

4 and 5 show the dependence of the magnitude of the pulses from the photomultiplier and 

the gas-discharge detector on applied power for helium and neon tubes. From these curves 
and the oscillograms it is evident that the quenching effect is observed in the helium tubes 
and in the neon tubes, but is stronger in the helium tubes. Thus, for tubes filled with neon 
the quenching pulse is observed only in the afterglow (Figure 5, Curvesland 2). In the 
helium-filled tube extinction of the recombination glow occurs both upon the action of the 
microwave field during the ignition pulse and after it (Figure 4). 

Several curves were plotted for clarification of the relation between detection and 
enhanced glow. Thus, Figure 6 shows the distribution of the glow and the dependence of 
pulse magnitudes from the gas-discharge detector and photomultiplier on the length of 
discharge in helium. The maximum of the detected pulse coincides with the excitation 
maximum; both of these maxima occur in the negative glow. Since the discharge spread 
beyond the electrodes, the pulses were also observed beyond the discharge gap. Curves 
were plotted for neon-filled tubes and were similar to those shown in Figure 6. 

Figure 7 shows the distribution of the extinction and the excitation of glow along the 
length of the discharge in a tube filled with helium. The maximum of the quenching effect 
and the maximum of the increase in glow are slightly out of coincidence. 

The dependence of the alternating voltage at the gas-discharge detector and the photo- 
multiplier output on discharge current is shown in Figure 8. The detector voltage has a 
distinct maximum as the discharge current is varied. 


3. EVALUATION OF RESULTS 


The experimental relationships derived reveal a similarity in the phenomena of 
detection and variation in glow. A qualitative explanation of these results may be obtained 
on the assumption that the microwave power changes the electron energies. 

If we assume (as in Reference 1) that the current depends on the electron energies, 
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where Ig is the value of the discharge cur- 
rent and Q is the electron energy, then 
an electron energy, increment of AQ will 
cause a current increment 


AT Sal,ee. (1) 


On the other hand, electrons possess- 
ing energy Q experience, according to 
Reference 4, the total number of collisions 
leading to excitation 


Pay 


dt xc 


With a change in energy by AQ the increase 
in the number of exciting collisions will be 


8 (Baxo= gel? +208) 


XC 


Figure 6. Pulse magnitude at output of the 
photomultiplier and the gas-discharge detector 
as functions of discharge length: p = 9 mm Hg, 
|? (2) helium, Tp = 10 usec, Tig = 75 psec: 


1 - AUpn(x); 2 - AUg(x); 3 - distribution of 
If it is considered that the intensity of glow Ug(x). 

glow ® is proportional to the number of 

exciting collisions, then a change in the 

number of collisions will lead to a change in glow by A®, which, in accordance with 

Equation (2), is 


AMex0= Py [+ +2 , 7 ° (3) 


i] 4 8 o 16 2,mm 0 04 08 42 46 2.0 
a K Figure 8. Dependence of pulse values at 
the output of the photomultiplier (1) and 
Figure 7. Distribution of quenching (1) the output of the gas-discharge detector 
and excitation (2) of the glow over the (2) on discharge current; p = 9 mm Hg, 
length of the discharge; p = 9 mm Hg, Helium, Tp 10sec, Tig = 75 psec, 
Helium, Tp = 20 usec, P/Py = 0.05. P/Py = 1. 


Hence it follows that an increase in energy leads to an increase in the excitation glow of 
atoms, wherein in the first approximation the increment of glow is proportional to the 
total excitation glow. This fact has been observed in actual operation (Figure 6). 

The decrease in glow upon the application of low power is most easily explained if 
it is considered that the recombination coefficient decreases with increasing electron 
energy in proportion to Q-l/ 2 [4]. Since the recombination coefficient indicates the 
number of electrons recombined per second, it characterizes the degree to which the 
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glow depends on recombination. Hence it follows that an increase in energy by 4Q will 
lead to a decrease in glow by 


ADyec= a = Drec o c (4) 
Due to the fact that the regions of intense recombination and excitation do not coincide, 
the functions Ad-e, and Adgy, have their maxima at different points of the gas discharge 
(Figure 7). 
As follows from Equations (1), (3) and (4), the magnitudes of the effects in detection 
as well as in enhanced glow are identically dependent on the increment in electron energy and, 
consequently, to this degree must be similar in nature. 


CONC LUSIONS 


The curves and oscillograms show that detection in a gaseous discharge is accompanied 
by a change in glow. The nature of the dependences of the change in glow on various param- 
eters (power, discharge current and site of the radiation) coincides qualitatively with the 
analogous dependences for detection. Toa certain degree this relationship explains the 
physical nature of the operation of the gas-discharge detector. Actually, while detection is 
not always accompanied by a change in recombination glow, the recombination coefficient, 
as this work has shown, does vary. Since a change in this coefficient may lead to a change 
in the number of free electrons, then the discharge current may also change to the same 
extent. : 

A change in the number of excited electrons upon application of the microwave field 
may also serve as a source of change in the discharge current (due to secondary processes). 
On the other hand, similarly, a change in the number of collisions leading to excitation also 
changes the number of collisions accompanied by ionization. This phenomenon will also 
cause a change in discharge current. 
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EMISSION OF HOT ELECTRONS 
FROM p-n JUNCTIONS IN SiC CRYSTALS 


M.1. Elinson, G.V. Stepanov and V.1. Pokalyakin 


Emission of hot electrons from p-n junctions in crystals of SiC is investigated as a 
function of the magnitude of the cut-off voltage at the junction Uj and temperature T. The 
existence of a high current density and a nonuniform distribution of emission from the 
surface is confirmed. The magnitude of the emission current ig is determined by the 
current through the junction ipey. The ratio y has a maximum at the specific voltage Uj. 
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INTRODUCTION 


SiC crystals are suitable objects for the study of emission of hot electrons, since 
€j > x [1] (€j - is the threshold energy of impact ionization, y - is the electron affinity). 
For hexagonal SiC ¢j ~ 4.3 ev, y =4 ev). Moreover SiC crystals are rather inert 
chemically which allows us to expect stability in their surface properties. 

In Reference 2 for the first time the emission of electrons from natural p-n junctions 
in SiC was investigated. A number of interesting features of the emission were found, 
in particular, a high current density je > 1 amp/cm2 and the law ig ~ ike where ig is 
the emission current; ipey is the reverse current through the junction. There exist 
several opinions on the nature of iyey [3, 4, 5], important for the understanding of the 
mechanism of emission. Mainly References 4 and 5 give a concept of incompletely ionized 
donor centers in SiC and their ionization which seems more plausible from a comparison 
with the artificial considerations expounded in Reference 3 concerning p-=n junctions. Un- 
fortunately, in Reference 2 we did not investigate the most important variations of ig with 
electric field in the junction and with temperature. Also the nature of emission centers is 
not clear. 


1. OBJECTS AND METHOD OF INVESTIGATION 


Just as in Reference 2 we investigated natural p-n junctions which form during the 
growth of SiC crystals by the method of sublimation. The existence of junctions was 
established from the well-known method of electroluminescence and by the volt-ampere 
characteristic. The crystals were cut from nodules, polished and selected by the number 
and brightness of luminous points upon switching on the cut-off voltage. The surface 
treatment consisted of the removal of SiOg dielectric films and outgassing. 

The experimental apparatus and measuring circuit are given in Figure 1. Crystal K 
was clamped toa thick molybdenum lead Ct with resilient tungsten points 0. The crystal 
temperature was measured by thermocouple T. The apparatus also contained a collector 
grid C connected to a willemite screen F onto a semitransparent platinum base. This 
base allowed the comparison of spots on the screen with the position of luminous points 
on the crystal. 

On to the crystal one could apply a static as well as a pulsed voltage. The latter 
enabled us to maintain a constant crystal temperature even at high reverse voltages U;. 

The anode voltage was adjusted over wide limits, and the emission current was 
measured by an electrometer with sensitivity aaLireS amp. 


Figure 1. Experimental apparatus and circuit. 


K - SiC crystal; O - tungsten clamping points; 
Ct - thin molybdenum rod for enhancement of 
heat exchange with the crystal; C - collector 
anode grid; F - fluorescent screen on a metallic 
semitransparent base; 1 - pulse oscillator 261; 
2 - pulse amplifier; 3 - device for measuring 
stationary or average current; 4 - device for 
measuring d-c or pulse voltage; 5 - vacuum 
tube electrometer; 6 - d-c current source; 

7 - anode voltage source; 8 - temperature 
measuring instrument. 


2. RESULTS OF INVESTIGATION AND DISCUSSION 


In a number of crystals, emission currents from 10-12 to 10-6 amp were obtained. 
As a rule, emission occurred in due course, which we produced with a continuous sampling 
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of current and with an increase of crystal temperature (~ 400 deg C.). After approxi- 
mately eight hours of warmup, the emission attained a stationary value and became 
exceedingly stable. The formation of emission is connected perhaps with the clean- 
liness of the crystal surface. 

The junction voltage Uj necessary for emission varies within wide limits for dif- 
ferent specimens. Comparison of the emission picture on a lumiscent screen with the 
pattern of luminous points on the crystal shows that the luminous points are emission 
centers. The large number and the great intensity of the luminous points are not 
definitely associated with a large emission current. The junction depth distribution 
near the emitting surface is the determining value. The number of emitting centers 
increase with an increase of U;. The linear dimensions of the emitting centers lie 
within the limits 1 {o 10 uw. The current density calculated from the total area of 
emitting centers is 1-10 amp/ em2. Hence, the most important result of Reference 2, 
concerning the high emission current density is confirmed. 

Figure 2 gives the relations irey = £(Uj) and ig = f(U}) for one emitter. The coarse 
structure of the characteristic for the reverse current (curves 1 and 2) consists of a sec- 
tion of rapid increase i, ,(AB) and a section of slower variation of iyey(BC) (see also 
Figure 3). Diminution in the growth of ipey is somewhat unusual for a cut-off voltage. 
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Figure 2. Comparison of static and pulse characteristics of a 
junction. 


1 - iyey(Uj) for the static condition; 1' - ie(U;) for the static con- 

dition; 2 - irey(Uj) for the pulsed condition, + = 20° C; 2' - ie(Uj) 

for the pulsed condition, T = 20° C; 3 - layer temperature during 
operation in the static condition. 


Curve 1 corresponds to a static regime in which the crystal temperature changed 
greatly with an increase of Uj (the temperature change is described in Curve 3). Curve 
2 corresponds to pulsed operation (Tj = usec), when over the entire range of U; the crystal 
is maintained at room temperature. : 

The corresponding characteristic curves of the emission current ig are given in 
curves 1' and 2'. 

The region AB corresponds to a rapid rise in ig. In region BC, the rise in ig is 
greatly diminished (more than iyey). For static operation, due to a continuous temperature 
increase, a smaller diminution in the growth of ie is observed. Only at T> 400 deg C there 
occurs a drop in ie, possibly connected with an increase of electron scattering by lattice 
vibrations. 

Figure 3 gives the pulse characteristics irey(Uj) and ie(Uj) taken for two crystal 
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temperatures T = 20° C and T = 75° C. 
On the basis of Figures 2 and 3 one may infer that the influence of temperature on 
hot electron emission under these conditions is weaker and of a different nature than that 


assumed in the theory [1], which applied to Ge and Si type crystals, especially in the 
region BC. 


Hence in the system studied, the in- 
crease of electron scattering with a rise 


in temperature T is compensated either ea ey : 
by an increase in the number of electrons, 
or possibly, by a change in the electric 197 
field in the junction upon varying T. 
The behavior of the characteristic ow 
also shows that the voltage at which virtual 
"saturation" of the emission current is 10°§ 
achieved is very far from threshold. 
The following established fact is 05 


especially interesting. If ig is plotted as 
a function of irey (see Figure 4), all ex- 
perimental characteristics, independent of Figure 3. Pulse characteristics iyey(Uj) and 
junction voltage Uj and temperature T, ie(Uj) at different temperatures: 
within the limits of experimental errors, . & 
agree with one another, i.e. , the experi- 5 ees nee i Bh iy he 
mental value of iyey always corresponds irey(Uj) a “tei 
entirely to a specific value of ie. 

This may be interpreted so that the electrical field at the junction E is very large 
and depends weakly on Uj (for example, E VUj;). Therefore the acceleration of electrons 
at the junction is always sufficient and is Pewually independent of Uj and T. The emission 
current is larger, the more the carriers at the junction. 
Figure 4 also illustrates the ratio y = ie/iprey, common for all characteristics. This 
ratio has a maximum in the break region of the characteristic irey(Uj ) (region of points B). 
Consequently, y depends on the behavior of iyey and can be explained by exterior 
causes, for example by junction emission in the condition of space charge. 

Evidently, in the region BC a new strong scattering mechanism comes into play, and 
the number of electrons capable of being emitted increases considerably more slowly than 
the total number of current carriers. 
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Figure 4. Characteristics ig(irev); Y (irev) 
1 - static; 2,3 - pulse, at 20°C; 4 - pulse, at 75°C. 
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The relation between iyey and ig also appears clearly in characteristic curves taken 
at liquid nitrogen temperature (see Figure 5; the curves refer to another crystal). It was 
established that a sudden increase of currents leading to the descending character of the 
relation is connected with heating of the crystal. 
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Figure 5. Static characteristics iyey(Uj) and ie(Uj) at -180°C. 


The rigid relation ie/iyey indicates that the back current passes in considerable part 
through those places of the junction, which are emission centers. It indicates the constancy 
of emitting surface properties. 

The maximum value is y ~ il Qieee i.e., verysmall. This is connected, evidently, with 
electron energy losses in the strongly alloyed substrate n-th layer in which the electric 
field is insignificant. 

The increased efficiency of this type of emission is one of the chief problems for 
further study. 


CONCLUSIONS 


In the study of emission of hot electrons from p-n junctions in SiC the following facts 
were established. 

1. The emission current density, in agreement with Reference 2, is rather large 
je = 1-10 amp/cm2. 

2. The emission is uniformly distributed over the surface. Luminescent points on 
the crystal surface are emission centers. 

3. The temperature variation of emission is less and has a different character in 
comparison with the theoretical relation for homogeneous Ge and Si-type semiconductors, 
which is connected with factors compensating for an increase in electron scattering by 
lattice vibrations when temperature is increased. 

4. The ratio y of emission current to reverse current through the junction is very 
small (~10-4), which possibly is connected with electron scattering near the emitting 
surface and has a maximum at a specific voltage U;. 

The decrease in y after the maximum is related to a new strong scattering mechanism. 

5. The emission current is closely connected with the reverse current through the 
junction, independent of the temperature T and voltage U; at which the experimental value 
of reverse current is obtained. This fact is explained by the existence of a strong electric 
field at the junction, always ensuring sufficient electron acceleration. 


256 


The authors express their gratitude to V.B. Sandomirskii for useful discussion, and 
also to N.V. Sumin and A.M. Fadeeva for assistance in the work. 
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FIELD EMISSION FROM TUNGSTEN SINGLE CRYSTALS 
IN THE PREBREAKDOWN CURRENT REGION 


G.N. Fursei 


Results are given for an investigation of field emission in the 'prebreakdown" period; 
the design of the apparatus and the experimental technique are described. 


INTRODUCTION 


It is well-known that at a current density above a certain critical value, time-stationary 
field emission passes into a "vacuum arc''-avalanche process, by melting the emitter [1]. 
This process leads to the destruction of the cathode, and hence it is interesting to study 
phenomena which lead to a 'vacuum breakdown" and its associated phenomena. Papers 
by Dyke [1, 2, 3, 4, 5], Elinson [6, 7, 8], as well as Shuppe and Gorman [9] in recent 
years have been devoted to this subject. 

These authors used a pulse technique which allowed the removal from the emitter of 
considerably larger current densities than in a stationary regime, and the observation of 
a number of phenomena, which would be difficult to study in the stationary regime. 

It was found that the arc is preceded by a number of phenomena, which are anomolous 
in comparison with the usual flow of the process;namely: 1) spontaneous growth of current during 
the action of a rectangular voltage pulse; 2) deviation from the form of the volt-ampere 
characteristic predicted by the Fowler-Nordheim theory which is well substantiated at 
smaller current densities; 3) the emergence of a bright ring surrounding the usual 
emission picture. 

The goal of our investigation was to check the existence of these phenomena, as well 
as the results published in References 1-8 and in Reference 9 which show a discrepancy. 

This discrepancy is that in References 1-8, the volt-ampere characteristic deviates 
from rectilinearity at small currents, but in Reference 9 - in the region of larger currents. 

Assuming that the technique used by us possesses a somewhat greater resolving power 
and sensitivity, we hope to obtain additional information concerning the development of 
the process. 


1. CONDITIONS AND METHOD OF INVESTIGATION 


In order to carry out the investigation, equipment (Figure 1) consisting of the following 
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main units was assembled. 

1. A generator of high-voltage rectangular pulses of duration 0.5 usec at a repeti- 
tion rate of 50-1000 pulses/sec. The generator could also give separate pulses of length 
0.5-10 psec with smooth amplitude regulation within the limits 0-20 kv and a rise time 
less than 0.1 usec. This generator included a type 261 pulse generator, which served 
as the source unit, and a modulator (M) of two tubes (GI-30 and GMI-83) which operated 
in a keyed mode of operation with partial capacitive discharge. The use of the 261 generator 
considerably simplified the assembly of the equipment because it was necessary merely to 
build a very simple modulator (Figure 2). i 

2. Double-beam oscillograph OK17M with 
a photographic attachments. 

3. Cathode follower (CF) with an attenuator 
at the input, which enabled us to change the in- 
put signal amplitude in steps up to a factor of 100. 

4. Heating circuit for the point (PO). 

A negative high-voltage pulse with an amplitude 
of several kilovolts and length two psec was 
incident on cathode K of the experimental vacuum 
device (Figure 1). The pulse of field emission 
current, which originated was detected by 

Figure 1. Block diagram of the collector C and proceeded through the cathode fol- 

equipment. lower to the first amplifier of oscillograph OK17M. 
The cathode follower was connected directly to 
the collector of the experimental device in order to minimize stray capacitance and in that way 
decrease the input circuit time constant, and also to match the total load resistance with the 
input amplifier of oscillograph OK17M (100 ohm). On the second pair of plates in the oscil- 
lograph cathode ray tube we imposed the voltage pulse taken from a resistance divider. 


+2ky tk IT n 00pF 
@ =a 


© @ 
+450 -800 
Figure 2. Schematic diagram of the modulator. 


Thus on the oscillograph screen one could simultaneously observe the voltage pulse 
and its corresponding current pulse at the collector. 

Triggering of the equipment was carried out by an intensifier pulse from the OK17M 
oscillograph. Synchronized operation of the equipment and oscillograph sweep circuit 
was ensured by an appropriate selection of sweep time in the 261 oscillator and oscil- 
lograph. 

The shutter of the photographic apparatus was carefully opened and in one frame an 
exposure was made at three-four points. 

The amplifier and plates of the oscillograph tube had been previously calibrated by 
means of a GSS-6, and also by means of a square wave pulse generator and a pulse volt- 
meter. Both calibration curves agree within the limits of error of the equipment used 
for the calibration. 

Measurements of the magnitude of current and voltage were made directly on the film 
by means of a comparator. 

Voltage resolution was restricted by the width of the sweep trace on the film and was 
equal to 0.05 mm, which corresponds to a voltage of approximately 50 v. Current reso- 
lution was also restricted by the width of the trace, but since with an insignificant voltage 
variation the current changed substantially the current lines were virtually always "resolved". 
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Errors in the current and voltage measurements were determined chiefly by the 
class of device used to calibrate the oscillograph, and could be estimated as 5-7%. 

Tungsten was chosen as the material whose emission and vacuum properties were 
to be most carefully studied. 

A cathode in the form of a point was prepared by the standard method of automatic 
electrolytic etching which enabled us to manufacture the point to extremely precise values 
of radii. This emitter was placed in a device which was subsequently evacuated. The 
vacuum system was improved by spraying a titanium and molybdenum absorber and by 
evacuation with an Alpert manometer. The measurements were made at a pressure of the 
order of 10°’ mm Hg. Before recording each point, the point was heated in order to 
clean its surface. Surface purity was controlled through emission pictures which were 
observed at the anode slit. 


2. DESIGN OF THE APPARATUS 


The main difference in the technique described above was the use of screening in 
the vacuum device. 

Great difficulty is the measurements is created by the existence of very strong pickup 
from the electrode onto which a high-voltage pulse is fed. To prevent this effect one may 
decrease the load resistance which makes it possible to obtain a time constant at which 
the pulse pickup is differentiated and timewise amounts to only an insignificant part of the 
useful signal (decrease of R), or via internal screening (decrease of C). 

The second method described by us is 
simpler, because, by taking a signal with 
a comparatively large resistance (10 kohm), 
one gould measure rather weak currents 
(10°° amp), by using the wide-band ampli- 
fier of oscillograph OK17M without any 
alteration. In addition, without changing 
the gain, and also without leading to ad- 
ditional calibration, one could change the 
input resistance by steps which allows a 
1000-fold decrease in the signal and a 
simultaneous improvement in the condi- 
tion for the passage of the current pulse. 

In this work we mainly used two types 
of vacuum devices (Figure 3). Their 
principal element was ascreening cylinder. Figure 3. Design of experimental apparatus. 
In design 1 (Figure 3) the point is located 
inside the cylinder close to the aperture 
of the anode slit so that the emitted elec- 
trons could not impinge on clear glass 
which would lead to its charging and to field distortion. The collector section with a 
flexible screen was connected to the input of the cathode follower. Design 3 was different 
from 1 in that the screening role was played by a metallic layer deposited on the glass. 


K - collector; A - anode cylinder; D - slit; 
N - cathode; O - window for observation of 
of emission pictures. 


r--- 


Figure 4. Collector section and cathode 
follower input. 
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Ig ifi,a) The screening cylinder (Figure 4) served 
ie q to shield the collector from pickup inside the 
device. Such a design allowed us to decrease the 
pickup approximately 1000 fold with a 5-15 fold 
reduction in the useful signal. Thus, the screen 
Gs succeeded in almost completely eliminating the 
effect of the pickup, which resulted only in an 
“3 insignificant distortion of the leading edge of 
r the pulsé. The point was situated at a distance 
of 5-10 mm from the aperture, so that it would 
be possible to project onto the collector the 
ereater part of the emission picture. 
af At the maximum sensitivity (R = 10 kilohm) 
one could confidently measure signals corresponding 
to a collector current of the order 10-6 amp. The 
aff time constant of the circuit allowed us to handle 
CT a ae ee 2 ae Dy Yel} pulses of useful signal with very little distortion. 
u At other positions of the input attenuator, corres- 
Figure 5. Volt-ampere characteristic of | ponding to resistances of 3.3, 1 Kilohm, 333, 
the field emission current during pulsed 100, 33, 10 ohms, similar distortions are absent. 
operation. The decrease brought about by the screening 
cylinder was determined experimentally for each 
case. For this purpose we supplied the cathode with a two usec length voltage pulse at a 
repetition rate of 250 pulses/sec and measured galvanometrically the average total anode 
current and collector current. By a similar method we recorded up to 20 points and then 
determined the average value of the attenuation coefficient. 


3. RESULTS 


Control measurements were made on 15 emitters in vacuum devices of a typical type. 
These measurements were allowed us to establish the following. 

1. At large current densities in all 15 cases, the volt-ampere characteristic deviates 
toward the region of smaller current value (Figure 5). 

2. In the "nonlinear" part of the characteristics, periodic current oscillations occur. 
A similar type of undulating pulsation of thermionic emission was observed in the presence 
of strong fields [10]. In Reference 11 these differences are explained by quantum-mechanical 
barrier effects of the interference of electron waves reflected from the potential barrier. 
We are inclined to assume that a similar type of phenomena also occur in our case. By 
confining ourselves for the present to the assumption stated above, we hope to conduct a 
careful analysis of this phenomena in the next paper. 

3. We observed a spontaneous rise in current with time (Figure 6), which was dis- 
covered at a current density of approximately 30-40% of the critical value which agreed 
with the results of References 1-5. 

4, During a calculation of the emitter radii according to the approximate formula of 
Reference 12, we calculated critical current densities which amounted to from 3(10)" to 
6(10)" amp/em2, which also is in good agreement with data obtained by other authors. 

5. Immediately before the emergence of the arc we noticed a spontaneous decrease 
of current at constant voltage, which supposedly can be attributed to the effect of space 
charge (Figure 7). 

6. As can be seen from the oscillogram of the moment of the emergence of an arc 
(Figure 8), the "arc process" is very short-lived and lasts less than 10-7 sec. 

The author gives heartfelt thanks tol. L. Sokol'skaya for recommending this work and 
for reviewing the manuscript. The author thanks Y.U. V. Zubenko and G. P. Shcherbakov 
for practical advice related with high vacuum technology. It is a pleasant obligation to 
express my Sincere gratitude toI.E. Isaev, Yu. V. Markov and the graduate A. A. Rogachev 
for a helpful conversation on the subject of technique. 
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Figure 6. 


Figure 7. 
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Emergence of spontaneous growth of current. 


Moment of spontaneous growth of current, 


decrease of current and the formation of an arc. 


Figure 8. 


aN 


Oscillogram of the voltage pulse and current 
pulse at the moment of breakdown. 
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ON THE THEORY OF THE STATIC CHARACTERISTICS 
OF A PLANAR DIODE WITHOUT AN INCANDESCENT 
CATHODE TAKING INTO ACCOUNT THERMAL 
ELECTRON VELOCITY DISTRIBUTION 


G. Ya. Myakishev 


Expressions are obtained which determine the static characteristic of a planar diode 
without an incandescent cathode for different regimes. In contrast with known papers we 
take into account the velocity distribution of electrons. Conditions are analyzed for a tran- 
sition from one regime to another. 


INTRODUCTION 


At the present time the static characteristic of a planar diode without an incandescent 
cathode is investigated in detail under the condition that the electron velocity distribution 
can be neglected [1]. Such an approximation, generally speaking, is permissible, because 
the mean thermal velocity is usually much less than the velocities electrons gain in an 
accelerating field. 

However with the formation of a virtual cathode, a calculation of the electron velocity 
distribution, especially in the region near the virtual cathode, may in a number of cases 
have great significance, in particular, during a study of diode noise and during the investi- 
gation of the influence of high frequency perturbations on the electron flux and retardation 
by the electric field [2]. The virtual cathode brings about a sorting of electrons according 
to velocity. 

It is interesting therefore to investigate the influence of the thermal distribution of 
electrons with respect to velocity on the static characteristic of a diode without an in- 
candescent cathode. It is of interest also to clarify the limits of existence of different 
types of diode operating regimes. 

We limit ourselves to a consideration of the most important case when the electron 
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beam accelerated by the grid falls into a retarding electric field between the grid and 
the reflector (reflex klystron). The expression obtained may be extended to the case of 
an accelerating field. 


1. STATEMENT OF THE PROBLEM 


Let the grid be located at the coordinate origin, and the reflector at a distance d from it. 
The electron velocity distribution and the grid plane, through which electrons pass in this 
diode, will be taken into consideration by means of the distribution function corresponding to 
a solution of the problem of a planar diode with an incandescent cathode* (for velocities v >0). 
Electron scattering by the grid structure is neglected. We will also assume that the reflector 
completely absorbs electrons incident on it (no secondary emission exists). 

Under these conditions the distribution function has the form 


mov maiees 
f (0,7) =O(@)N, V2 sane “HE for v>0, 
TGs v0 for v<0, 


(1) 


1 is the grid- to -cathode potential difference; No - is the electron concentration; @(v) - 
is a cut-off factor that takes into account the fact that the electron velocity at the grid 
cannot be less than specific values: 

Pt6r yiz Uae 


0 (v) = 
0 for v =< Vor Vi; 


(2) 


Vi =©1+ 2, where (0 is the depth of the potential minimum relative to the cathode in 
the cathode-grid space. In the absence of a virtual cathode in this space Vj =}. 
The boundary conditions for the potential are chosen in the form 


gO)=Gq, g(d)=q, (3) 


where (9 is the reflector potential relative to the cathode. In this case (1 2 (9. 

To solve the problem we use the system of equations adopted by A.A. Vlasov [4], by 
assuming that inter-electron collisions can be neglected. 

As a result of the plane symmetry of the problem and its static nature, the equations 
for the distribution. and potential functions have the form 


Of (a.v) , € dp(e) Af (er) _ 9 
: Ox Mf, GHB Ov df 

29 (&) © 

: oe = Ane \ } (a, v) dv. (4) 


—cao 


The sign of the electron charge is taken into account here 
We introduce further dimensionless magnitudes: 


pes ae 
en 7 ay m o é ©) 


*Such a function better represents the actual physical picture than the Maxwellian drift 
function discussed in Reference 3. 
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Omitting the primes in what follows (for simplicity of notation). The original system 
of equations and boundary conditions may be rewritten in terms of new dimensionless 
quantities in the following manner: 


+co 
of , dp af Ce te 
Gul de var war = \ 14, 
4 ~S +e, 
7(0, v) = O(z) Vea e for v>0, (6) 


(1 for oo Ve. 
~ 10 for v< VV, 
ji BD) =0 for v <0, 


OQ (v) 


PO)=q, 9) =, 


4 / 2 
Ll=d me N is the dimensionless grid-refiector separation. 
2. DIODE OPERATING CONDITION IN WHICH A VIRTUAL 

DIODE EXISTS (REGIME III) 


In this case it is convenient to introduce a new system of coordinates, choosing the origin 
at the potential minimum. The potential will also be reckoned from the minimum. Denoting 
the new coordinates by € and the new potential by 7, one may write 


E=2—2m, n=o—Vm, (7) 


where Xy, is the position of the potential minimum relative to the grid; Vi, is the depth of 
the minimum relative to the cathode. The approximate shape of the potential curve is 
illustrated in Figure 1. 

In the new variables, system (6) takes the form 


@ ze =e a “ =0, nm) =m, (Se) = Ne, 
a1 le f (Eq, v) = 9 (r) : ie for v>0 
dé? — ’ (1)> V on : (8) 
f (Ee, v) =0 for v< 0, 
{for yess) eV, 
O(v)= ae z 


0 for v< V 2V,. 


Here € (1), & (2) and #1, 72 are the coordinates and the grid and reflector potential re- 
spectively in the new system. It is obvious that 


= (S@)) 45 yo) =U 
Ui UE NO x (9) 


Equation (9) together with the system of equations and boundary conditions (8) completely 
define the problem. 

To solve system (8) we use methods developed by L.A. Vainshtein [5] and L.E. Parga- 
manik [6]. 

The general solution of the kinetic equation of system (8) is an arbitrary function of 
the magnitude > — 7 (6): 


f= H\>—n@}. (10) 
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A solution of Equation (10) must satisfy boundary conditions (8) and consequently, f(é,v) 
should be different from zero only over specific velocity intervals. 

In the diode space, in the case investigated, 
it is necessary to distinguish between three re- 
gions: yj <n <7 ,-Vj, first region ; nj - 

- V1; << 0 second region ; and0 <7 < 79 
(Third region) (Figure 1). 

For the second and third region, the solu- 
tions whose accuracy is of the order of the 
factor e?1 agree with the corresponding ex- 
pressions of the theory of a diode with an incan- 
descent cathode [5]. They can be obtained after 


integration of the system of Poisson equations Figure 1. Qualitative picture of the diode 
(8). potential distribution in the presence of a 
For the second region virtual cathode (regime III). 
RS a eS ees 
a1 | Te \z Vy Vy) ay 
iS 
4eti— m1 2 \ A 
+ \erVy Fad), (11) 
Va on J 
Tiesto eA EO 
ae ee a —Y = a 
ti = —e Vy fen Vy+n—Vy)dy 4 
es oa 
ss \ evVytn ndy | dy. (42) 
Vx Sp 
For the third region 
as 
Mm—F1 | 4 eo 4 a 
reve 2 4 Baia Vay Seen VS 

Eo =e \ Vi | eu Vy+n—Vy) ay| dn. (13) 


Tables have been compiled [1] for values of gy , and €5 as a function of 7). 
We consider the solution for the first region. an this region electrons may have velocities 
within the following ranges: 


co DevSV2V,—m+n) for »>0, 


(14) 
—V an Sv >— V20,—n +0) for vy<0. 
Consequently, for velocities within the ranges (14), the distribution function satisfying 
boundary condition (8) is 
£ oe. 4 —s. N—Nit1 15 
JAS: )= Vox ( ) 
For all other velocities f (€,v) = 0. 
The equation for the potential takes the form 
ee 00 eyes = Vain) OL 
alee ee ae \ ees Ra pred @) ae ah 16 
ds" a V 20 \ dv 1 V 2n \ v ( ) 


V3(Vi— aren) —V% 
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iw) 


Making the change of variable y = =. we obtain 
oe) 0 
dn ert te e “dy e eYdy 
SITs a Gas a — 2 ee Ang 
ae VV 2a \\ Voutm shay V2(y+n) 


Integrating Equation (17), we obtain 


Pi M1 “2 ——— = 
2-145 euVytn—Vy)dy+ 
1—T1 2 a — $1 
fea \ evVy+n—Vy)dy+cl®. (18) 
Va —(m—V1) : 


The integration constant C may be found from the continuity condition of the electric field 
at the boundary between the first and second region. Using Equations (11) and (18) we obtain 


0 


4ePr—1 — 
Ce ee | evVyay. (19) 


Tv « 
ie —(ni—V1) 


Integrating Equation (18) and determining the constant from the potential continuity condition 
we finally obtain 


m1 71 oo 
: i A ea wt = 
Gea nee” \ | =\ eV y +n— Vy)dy+ 
™M—Vi1 Va 0 
Z 0 aly N—1 
oa Vy + nay} "dy —e * x 
1 
—(m—V;) 
Duty ifadeals arm Page 4 a 
ce at y)dy + = \evVy Endy} dn. (20) 
0 0 —7 


An equation determining the parameter 71 may be obtained, by substituting values of &€ from 
Equations (20) and (13) into Equation (9). In developed form, this equation is 


Ni mh os 
EE te Dab 


San are) bea A dare 
+7 ( eo Vy+ndy) 


rf oy) 
—(m—V1) 
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0 0 


1 
= 5g) 
o1Vyrn Vi) ay| dy =I. (24) 
These expressions determine the static characteristic of the diode for this case. The solution 


is applicable also for @o < 0. 


3. THE OPERATING CONDITION OF THE DIODE AT WHICH A POTENTIAL 
MINIMUM EXISTS, BUT A VIRTUAL CATHODE DOES NOT ARISE (REGIME II) 


We find a solution for this case by the method used in the preceding paragraph. We 
distinguish between two regions of the diode (Figure 2): 0<7 <j (first region) and 
0<7 < 12 (second region). 
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For the first region 


SA yE oe de ——————— 
es ae z for oDS>u> V20i—n1. +7), 
(22) 
7(§, v) =0 for v<V2(V,—m, +). 
The equation for the potential agrees with Equation (16), if in Equation (16) we omit the 
second term of the right-hand side. After simple transformations analogous to those per- 


formed earlier, we obtain 


ee) 
dy eFi—™ e "dy 


ae Vand V2u+q) at 
After integration of Equation (23) we find € as a function of 7: 
a: 
1-1 1 9 0 ~ 4ae 

=e? | \ e1(Virn—Vi)ay| dn. (24) 


A solution of the second region in no way differs from Equation (24) excluding €9=-€1. 
It is not difficult to see that when 1; = V1 Equation (24) agrees with Equations (13) and 
(20). 


Figure 2. Qualitative picture of the po- Figure 3. Qualitative picture of the po- 
tential distribution in a diode under the tential distribution in the diode in the 
conditions that the potential minimum absence of a potential minimum. The 
exists, but no virtual cathode emerges case when 09 > ~)q 

(regime II). (regime I). 


The equation for determining 71 which agrees with Equation (21) when 71 = Vj is 


wle 
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m—F1 11 © a 
ene eee \ eu(Vy4+n—Vy)dy} dnte” x 
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eT ( ) aa -V de] | aren (25) 
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4, OPERATING CONDITION OF THE DIODE AT WHICH THE POTENTIAL 
DISTRIBUTION DOES NOT POSSESS A MINIMUM (REGIME 1) 


At first we consider the case when 29 2 “Pg. Let the origin of the coordinates 
coincide with the reflector plate (Figure 3). The potential ¢ also will be reckoned from 
the reflector. Then the solution of the kinetic equation satisfying the boundary conditions, 


assumes the following form: 
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f(t, 0) =e * for © >v>V2(Vi;—G+%4 9), 
23 26 
(Go) = for »<V2(V,—GQit+e+ 9). a 


Substituting this function into Poisson's equation and making transformations analogous 
to the previous ones we obtain 


to{e 


il 2 - ——— 7 7 
=e I \ HV TES AY + Cs) (27) 


Vi-(@i—#2) 


where Cj is the integration constant. Integrating with respect to~, we obtain x as a 
function of ~: 


co tan 


Ce oz \ eu Vy + ody + c,| dq.- (28) 


0 Vi=(01—#2) 


1 
> 


The constant C1 is determined from the condition 


es) Fa \ Wy + gay + Chee? (29) 
0 Vi—(¢1—2) 
At small electron densities in the beam l — 0 and consequently, Cj +”. Equations (28) 
and (29) in this case lead to a linear potential drop: 


Pier, 


Oe i (30) 


We note that Equation (29) must agree with Equation (25), if in Equation (25) we set 
N1=€1-. (a potential minimum is found at the reflector; regimes I and II coincide). 
Consequently, in this case 

. o9 a " 
(eS = Fe \ ou V y dy. (d1) 


Vi—(41—#2) 


This value of Cj is minimal for fixed V1, ~1 and 
©9 , because if the opposite were true 1 ceases 
to be real. If, in addition, V1 = 1 -9 (all 
three regimes coincide), Cj = -1. 

We now proceed to the case 9 < Dg. 
For this condition in the diode space, two 
regions should be distinguished (Figure 4): 

1 + 92 > 201 + Og - Vj (first region) 
and 1 + 2 - Vj > > 0 (second region). 

In the first region 


Figure 4, Qualitative picture of the 
potential distribution in the diode in the ” 
absence of the potential minimum. The 1 RSE ee Deas 

case when Qo < - Yq (regime I). eae) Von 
f(x, V) = 0 for all the remaining values 
of Vv. 

f er Vides ae 

After substituting Equation (32) into Onooiery sale Vile stain alsa ane?) 
the right-hand side of the equation for the 
potential and making the same transfor- 
mations as in the preceding cases, we 


obtain and —-V2(Vi— 4; — @2 + 9) >v>—V 2%, 
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ce 0 5 


Po 
Ca ine ean es f ‘ SS 
ia ran eUVy + ody + ie \ eu Vy +ody+ c,| ; (32) 
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C2 is the integration constant. Hence 


Pe q foe) 
i= — on | az. \ ev Vy +ody+ 
Oiter—Vi Va 0 
0 ee 
4 ° 4a 2 
Sara . é Vy +@dy + c,} dp — Xoo, (34) 
—(9:+¢2—Vi) 


where xXq@ is the value of coordinate x at the point  =91 + 2 - V1. 
In the second region 


Get hy 
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e 72 V. 
f(a, ”) = Vax? for 0 >v>—V 29, 


; pa (35) 
f(x, v)=0 fORRe <=) 00), 
The derivative of the potential corresponding to function (35) is 
1 
dg = oxic! ~——. Lops tasty ale 
Sig@iGct ae fo i 


It follows from the continuity of this derivative that C, =Cy. Integrating Equation (36) 
we obtain 
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ore. 
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Using Equations (34) and (37) we write an equation for determining Co: 


2 Atos 5 oo 
of VLR verted + 
Oito2—Vi 0 
0 -+ 2 @te2—Vi 2 co 
ae evVy+ody+ c,| dpe” \ ayer uF odu+ 
Vu a Vu. 
—(4:+¢2—V1) 0 0 
0 ip 
i ose \ ewVytedy+ c,| dp =I. (38) 
Vu 


EO. 


Equation (38) should change to Equation (21) ifm, = 1-9 and@9<0. This will 
happen when C9 = -1. This value of Co is minimal. 


5. CONDITIONS UNDER WHICH DIFFERENT DIODE 
OPERATING CONDITIONS EXIST 


First of all it should be pointed out at what values of the retarding potential (59 different 
diode operating conditions can be realized. 

A solution corresponding to a virtual cathode in the grid-reflector space, can hold at 
both positive and negative potentials Qo. 

Regime II can be realized only when (9 % “Pg, as may be seen directly from Figure 2. 
Regime I may exist both when (9 > -Po and when $2 < Ho. The solutions for these cases 


have a different form. 
We now proceed to the main question: at what initial current densities are the different 
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regimes realized, if the grid and reflector potentials are fixed? We consider the case of 
greatest interest @2 > ~)g, when all three regimes may occur. : 
In the form of solution found by us it is convenient to consider no initial current density, 


and J =d ane is the dimensionless grid- to reflector separation, introduced previously. 


kT 
In this expression Ng is a magnitude directly proportional to the current passing through the 
grid: 
ae 
j= eN, Vee sia (39) 
It is necessary to show at what values of 
the dimensionless distance ] these regimes may 
be realized. The possible values of ] at fixed 
1, Og and Vj are determined from Equations 
(21), (25), (29) and (38) which express the 
dependence of 7 on the parameters 71 is apparent 
from Figures 1 and 2. For regime III 71 2 V,; 
for regime II 01 -2<17 4< Vy. In the case of 
9 ~C=Cimin Gj regime I 
Figure 5. Approximate character of the tie t Sr ea (40) 
dependence of the dimensionless grid-to- nee Van \ ery ydy. 


. Vi—(1—92) 
reflector separation on the parameters er 


nN, and Cj for different diode operating In Figure 5 curve BC illustrates the ap- 
conditions. Curve AD corresponds to proximate variation of / with 7, for regime III. 
regime I; AEMB-regime II; BC-regime According to Equation (21) Lt increases mono- 
Ii. tonically with n;. Consequently, regime III can 
be realized only when l > OK, i.e., under the 
condition 
1 
ova foe L Re, 
i>e* \ ara ev(Vy+n—Vy) a dy + 
g Re, 
Go Vi-(¥:1— 2) BS co a -+ 
+e* \ ral ev (Vy 1—Va) ay] dy. (41) 
0 / 


Curve AB (Figure 5) illustrates the approximate behavior of 1(n;) for regime II (according 
to the results of numerical integration). Of paramount importance is the existence of a 
maximum at point M. The possible values of ] for regime II lie within the interval OL <]<OP. 
In Figure 5 the range of existence of regime I is indicated. According to Equation (29) 1 de- 
creases monotonically as Cj increases. Regimes I and II coincide when Equation (31) and 
the condition nj = © 1 - 2 are fulfilled. Consequently, permissible values of J] can be 
ar by curve AD. The first regime corresponds to values of ] over the range 
O<l<OL. 

Curves BC, AB and AD give the limits of existence of the different regimes and explain 
the hysteresis character of the current variation at the reflector as a function of initial 
current. At small currents regime I is realized. As the current increases continuously 
there occurs a transition to regime II without a current change at the reflector. At current 
values corresponding to ] = OP (point M on curve AB), the diode jumps in steps to regime 
III (point N on curve BC) and the reflector current drops abruptly because a virtual cathode 
is formed. With decreasing initial current, the diode operates in regime II down to J = OB 
when the diode changes by stages to regime II (point E on curve AB) a transition is also 
possible in regime I at several values of 1 and 9. In this case there will be no abrupt 
jump in the reflector current, because at 7; = V, the current of regime II equals the 
current of regime III. Section MB of the curve cannot be realized. 

Such a hysteresis type character of the reflector current should be observed also in 


> 
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the case of a constant current through the grid, and in the case of a changing potential 9. 
We note that the experimentally observable current jumps have a sharply pronounced 
character only during a diode transition from regime II into regime III. During an inverse 
transition the current changes rather quickly, but in spite of everything smoothly. This 
result agrees qualitatively with the theory expounded above. We note again that when 
2 < ~Po regime II is not realized and, consequently there should be no hysteresis relation 
for the current. This conclusion is also found to be in accord with experiment. 
In conclusion I thank A.I. Kostienko, R.L. Stratonovich and M.N. Devyatkov for their 
interest in the work and useful discussions. 
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TRANSITION PROCESSES IN SEMICONDUCTOR 
DIODES USING NARROW 
FORWARD CURRENT PULSES 


é 
Yu. R. Nosov 


We consider pulse parameters describing the transition process, when switching a diode 
from the forward to the reverse direction. Relations are obtained for these parameters as 
a function of the duration of forward current flow. 


INTRODUCTION 


In papers devoted to a consideration of the process of switching a semiconductor diode 
from a conducting to a nonconducting state, the duration of forward current flow is usually 
assumed so large that the hole distribution in the diode base * may be regarded as in the 
steady-state [1, 2, 3]. 

In practice, however the use of diodes in fast-acting circuits requires us to frequently 
operate with extremely narrow forward current pulses. Consequently, it is of interest to 
calculate how one may vary the parameters describing diode transition characteristics in 
the case of a nonsteady state hole distribution in the base. 

On the other hand a development of the technique of measuring pulse parameters of 
semiconducting diodes often requires us to transmit through the diode extremely narrow 


*Hereinafter it is assumed that the diode base is an n-type conductor, so that its 
minority carriers are holes. 
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forward current pulses before switching (in pulse generators an increase of the steepness 
of the rise time is connected, as a rule, with a reduction of the possible pulse length). For 
this purpose it is important to determine the minimum duration of the forward current 
necessary to achieve a steady-state distribution of holes in the diode base. 

The goal of this paper also is to clarify the influence of the duration of a forward current 
pulse on the character of the transient process upon switching the diode from the forward 


to the reverse direction. 
1. FUNDAMENTAL ASSUMPTIONS AND EQUATIONS 


We shall consider the model of a one-dimensional junction diode with an abrupt p-n 
junction and with a semi-infinite base region, the p-region conductivity being assumed much 
larger than the n-region conductivity (see Figure 1). Just as in References 1, 2, 3, the 

injection level is assumed to be small, but the lifetime 
of holes in the base is assumed to be a constant. The 


influence of the electrostatic capacitance of a p-n junction 
P ‘ on the character of transient processes is excluded. ** 
These assumptions are satisfied, to some extent, by 
e & domestic industrial types of pulse diodes which are 
Figure 1. Model of a one- pseudopoint structures: (for example, diodes D9, D11-14, 
dimensional junction diode. D101-103). Microjunction silicon and germanium quick- 


acting diodes conform well with this model. 
The behavior of holes in the diode base is described by the diffusion equation 


oP _p Pp PP Pno (1) 
ot 2 Ge Uy 


where p is the density; D, is diffusion coefficient; Tp is the lifetime; pyo is the equilibrium 
density of holes in the base. 

Neglecting the magnitude of pnd in comparison with p and introducing dimensionless 
time 6 = t/Tp and distance X = x/Lp variables (Lp is the hole diffusion length in the base) 
we obtain a diffusion equation in the form 

SS (lay 


The current density through a p-n junction is expressed by the following equation: 


qD py Op 


| = SS y 
Tiron \ a, (2) 


In the steady state, the hole distribution in the base has the form 


qUe¢ /kT —x/L. —A 
fs e ee Pol - (3) 


Ps Pro? 
Where Uf is the forward bias on the diode; q is the electron charge; k is Boltzmann's 
constant; T is the absolute temperature; pg is the hole concentration near the p-n junction 
in the steady state. Consequently, for a steady-state current flow we have 


: qD Po 
,= en : (4) 
The solution of the problem formulated in the introduction is broken down into two stages: 
l)a calculation of the establishment of a hole distribution in the base, upon connecting 
the diode in the forward direction; 2) a calculation of the junction diode transient charac- 


teristic upon switching it to the reverse direction taking into account the distribution obtained. 


**Nosov [4] has more thoroughly formulated the basic hypotheses underlying the modern 
theory of transition processes in semiconductor diodes. 
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2. DETERMINATION OF THE HOLE DISTRIBUTION DURING 
THE PASSAGE OF A FORWARD CURRENT PULSE 


A step in the forward current density from zero to jr leads to the emergence of a hole 
concentration gradient close to p-n junction given by 
Op lf L, 


aN aie = Po. t 0. (5) 
N= mp 


Up to the starting moment of a forward pulse, the hole concentration in the diode base is 
negligibly small: 


TINEA Oat ne = ==10) (6) 


We solve the diffusion equation (la), assuming the initial distribution (6) and boundary con- 
dition (5). With the substitution 


p= U(X, A) (7) 


Equation (1a) reduces to the form 


au aU 
ze ) 
eau ia ps (8 


Integrating Equation (8) by means of a two-dimensional Laplace-Carson transformation [5] 
and taking into account the substitution (7), we obtain 


a 


P= Poy je ol ave Ge Vo) — Sey One eae : (9) 


Figure 2 gives graphs of the relation £, (X) for different fixed values of 9 
From solution (9) it is easy to obtain 


Pyap = Po at Va. (10} 


Hence the voltage rise at a p-n junction is given by 


a ea Mes (Niee ae ee ferf YO], 
q Pa 


U f () oa 


where Uf-s is the voltage at the p-n junction in the case of a steady-state hole distribution 
in the base. 
Because Uf-g is usually of the order of 0.25-0.5 v and at room temperature kT/q = 
0.025 v, the voltage at the p-n junction reaches 0.8 of its steady-state value, under the 
condition In [erf V6] = - (2-4), which leads to © = 0.02-0.004. So small a value of © explains 
the fact that no steady-state voltage at the p-n junction is experimentally observed, because 
it usually takes place during a time less than the resolving time of the apparatus. 
An examination of the graphs of Figure 2 shows that when © < 1 the hole distribution 
in the base may be described by an exponential relation of the form 


v0) (12) 


JOAv (Gy = CD es 


We will define p(@) and/(6). Using a decomposition into a series of error integrals, we 
obtain from Equation (10) 


P ay 7 V9 p (13) 


2 ; 2 
*Here and hereinafter it is assumed that == 1 —erfz=1— =| pty 
x 6 
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Integrating Equation (1) over the volume of the 
lee) 


base and noting that q f pdx = Qy is the total hole 


charge accumulated in the base per unit area of the 
p-n junction, we obtain the following equation for thi 


BS process of accumulation: 
A, (0) = On (8 — 0c) (1 — e-9), (14) 
=10 —~0-= where Qy (8 ~~) =qpolp. When 6 « 1 we have 
" 
~ abi Oy (8) = Qu (A— 00) @ = qpoLp9. (15) 


Figure 2. Distribution of the hole density 
in the diode base for various durations of Integration of Equation (12) multiplied by q over 
a forward current pulse. the base region gives 


Oy (O) = p(O)L(9)q. (16) 


Equating expressions (15) and (16) and taking into account Equation (13), we obtain 


1(0) = UE by. (17) 

Thus the hole distribution in the diode base at @ « 1 can be described by the expression 

Be eta eS 
Pp (xX, 0) — Po == Ve oe (18) 
Va 
Numerical calculations carried out at 6 = 0.1 show that Equations (18) and (9) differ by 
no more than 3% for arbitrary X values. 

Letting 2/V7 ~ 1 as an approximation we note that the hole distribution in the diode base 
for a small forward pulse duration (t « Tp) is the same as for the steady-state case if the 
hole concentration pg near the p-n junction and the diffusion length Lp are decreased by a 
factor of 1/V0. 

At the end of this paragraph, we will find the establishment of the hole distribution in the 
diode base during the instantaneous transmission of forward voltage* at the p-n junction**, 
i.e., we willintegrate Equation (1a) with condition (6) and the boundary condition 


Ply = Do LOL ate 0s (19) 
The solution is obtained just as for the case of constant current and is 


p(X, 0) = Ble* exfe (75 — V8) + eX erfe es +y 6) (20) 


The current density through the p-n junction during the process of establishment is 


i=ie [Sy et VO | (24) 


and at small values of 9 equals 


; F a 
en ee) (21a) 


*The magnitude of the forward voltage is taken so that, during establishment of equilibrium, 
the same current density jg_, flows through the diode as in the previous discussion. 
**Such a calculation was made in Reference 6 for a diode with a thin base. 
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The accumulation of hole charge in the diode base is described in this case by the 
expression 


t 


? u (0) z \ jatj= ve Ou (O > fered) Vo 4 (22) 


0 


Comparison of Equations (15) and (22) shows that the rate of accumulation of holes in 
the diode base, immediately after switch-on (i.e., © K 1), upon applying a voltage pulse to 
the p-n junction, is 1/VO times larger than for the case of a current pulse. Physically this 
result is explained by the fact that the pulse strength of the current density during the process 
of attaining equilibrium, at all times, is larger than the steady-state current density. 

It should be noted that in real electron 
circuits the diode forward resistance, as 
a rule, is much less than the load resis- 
tance connected in series with it, i.e., the 
conditions for a current generator are 
realized. Therefore, in future considera- 05 
tions we use the solution obtained with 
boundary condition (5), but not (19). 


Graphs of Becton different fixed values 


DWDy 


0 1 


of 8, given in Figures 2 and 3, show that 10 20. X 

with small © the distributions p(X, ©) for Figure 3. Distribution of hole density in the 
the case of a switch-over to the voltage diode base for different forward voltage pulse 
generator regime are substantially dif- durations. 

ferent. 


3. SWITCHING OF A DIODE FROM FORWARD TO REVERSE DIRECTION 


When switching a diode from the forward to the reverse direction, a reverse voltage 
pulse is applied to the p-n junction and the series resistance R (equal to the total resistance 
of the diode base and load). 


~ 
a 


6, 
J a 
Figure 4. Theoretical transition char- Figure 5. Transition characteristic of the 
acteristic of a diode when switching it diode at a finite rate 
from the forward to the reverse direction. of switching. 


The form of the transition characteristic for this process is shown in Figure 4. The 
limiting case for R +0 [3] leads to the disappearance of the ''shelves" and to an infinitely 
large value of reverse current at the instant of switching. Two circumstances, however, 
substantially change the initial stage of the switching process. 

The finiteness of the maximum possible velocity Umax limits the hole current through 
the p-n junction [7]. This hole current is directly proportional to the hole concentration 
near the p-n junction py-g. The other limitation on reverse current is related to the finite 
rate of any real switching. In this case the maximum value of reverse current is also 
determined by the magnitude of p,=g, although the direct proportionality is violated. 

Because of these limitations, a real transition characteristic of the diode has the 
form shown in Figure 5. 

A quantitative calculation of the pulsed properties of a diode is made via the introduction 
of three parameters: 

1. The magnitude of the maximum value of reverse current during switching Imax which 
characterizes the initial stage of the transition process. 


275 


2. Recovery time of reverse diode resistance Treg determined by the time interval 
elapsing from the beginning of switching up to the instant the diode reverse current drops 
to a certain prescribed value. The recovery time of the diode back resistance usually 
equals (3-4) Tp and thus characterizes the "tail" part of the transition process. 

3. The magnitude of electric charge passing through the diode in the reverse direc- 
tion during switching Qp. 

If we assume that Imax ~Px=0 is correct, then 


Tmax (Qo) = Imax (Qo — 20) erf V 80; dmax (Qo) = 0,8 Tmax (Ao > %&) 
LOCKE 10535 
In order to determine the dependence of Treg on the duration of a forward current pulse 


Oo0* it is necessary to find the entire diode transfer characteristic. 
We will integrate diffusion equation (1) taking the initial distribution 


p(X, 0) = p(X, 0,) = Be erte (= —VB,) — e¥ erfe (Fat Vo} (9a) 


and the boundary condition eS 
iO, @)=0, Csr (2. 


Assumption of the boundary condition (23) indicates that a back voltage is instantaneously 
applied to the p-n junction. This assumption is correct, if it is considered removed from 
the initial stage of the transition process (as in our case) and the resistance in the diode 
circuit is so small that the length of the ''shelf'' determined by this resistance is insignificant**. 
In most cases using pulse diodes, the assumption of condition (23) is most reasonable. 

From the general solution of the problem, obtained using the original function for a semi- 
infinite region, we find 


oe ° - 
= ar D(X, Oo) e 4X dX, (24) 


ik, 
ax laa, 2 
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We transform the initial distribution (9a) to the form 


(-) — me 22 aes = — 2s ) zie 
p(X, Oo) = Poe * 5 fe ¥ erfc (V9, | 
ee ‘ a XxX ee -y r 9° 
+ e* erfc (V2 + We) = poe * — 9(X,,)). (9b) 


Solution (22) at p(X, 0) = (KX, Qo) is discussed in detail in Reference 3, where it is shown 
that for large values of 6 we have 


é xX aXe sas 0 wil c 3 Yo : DD 
aches \ Q(X, Ao) € d I7'/2Q°/2 | a /2 r erfe Ao ( 
Q 
Under this same condition 
w PS 
XG Ac = Se = Q% 
e SiO 7) 26 
sang | Patt on aaa aatnetiat Ce 
0 


*Hereinafter the duration of the forward current pulse 99 is a parameter, but 9 is a vari 
able magnitude measured from the beginning of the switching. 

**The subject of the quantitative relation between the length of the "shelf"! and the diode 
circuit resistance is considered in Reference 4, 
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Taking into account Equations (24), (25) and (26) the density of the diode reverse current 
at large O* is expressed by 


: barat 4 op — # cme 205°emie! ai in @'/2 a 
years mara ae ale ee f ( ihe erfc 0 |: (27) 


For a large duration of forward current (69 ~@) this relation has the form [3] 


ae (28) 


Trev = ¢ oe: ” 


The graphs of Fig. 6 illustrate Equation 
(27). Itis easy to see that the recovery time 
of the back resistance reaches 0.8 of its 
steady-state value in the case of a forward 
current pulse length of O99 ~ 0.8. 

Thus the transition process of switching 
both in the initial (specific parameter 
Imax) Stages and in the "'tail'' stages 
(determined by the parameter Tyec) are 
practically independent of the forward 
current duration, if it exceeds 0.8 T,. 

To determine the amount of charge 
passing through the diode in the reverse Figure 6. Variation of jyoy/j¢ with 9. 
direction during switching, one must 
know the entire transition characteristic which for arbitrary values of 69 was not found 
by us. 

We find the transition characteristic of a diode for extremely narrow forward current 
pulses (69 « 1), using the approximate initial hole distribution (18)**. 

The solution found in this manner will differ from that for the case of the actual initial 
distribution, however, by virtue of the correctness of the statement of the boundary 
problem of diffusion equation [8], the difference between these solutions will be no larger 
than the difference between solutions (18) and (19). Carrying out an incomplete inte- 
gration of the diffusion equation, we find by means of an operational method [5], the value of 
the gradient of the hole concentration near the p-n junction 
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Using Equations (2) and (29), we determine the amount of charge passing through the diode 
in the reverse direction after switching, 


co 


; Oy Oo 

a = , — pada Otll HS > = = 4 30) 
On (Mo) =\ Jreyat = (Pole == = Qn (A > ~) i ( 
rd 1+ Yabo 1+ a: Ve 

where Qy (E90 +) is the total charge of the holes accumulated in the diode base. In 

view of Equation (15), we obtain 


fn es Q, (®o) i Va i 
= -) eS ee eo ee / @ 31 
3 (Oo) Q,, (Po) Airs V 1 Ve =! 2 V ®&. (34) 


This calculated expression for € (99) determines the correction to the pulse method of 
measuring the injection coefficient y [9] and should be considered when setting up 
appropriate experiments. 


*Calculations show that when © > 3 the error is less than 20% in comparison with 
S) 2) 
**Here, as previously, boundary condition (23) is assumed. 
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As is well-known, when 09 ~@ & (89) 1/2 and, consequently, the possible values of 
€ (QQ) lie between 1 and 1/2. Since the total charge Qy passing through the diode during 
the passage of the forward current pulse equals Qp = I¢to = Qu (C9 +”) 80, taking into 
consideration Equation (14), we obtain 
als Rens (32) 
nS Ge g Fick 


The coefficient n* characterizes the rectifying ability of the diode when using it in the 
switching regime. As is apparent from Equation (32), when @g = 17 (@9) = 0.3 — 0.5 and 
only when @g = 5 7 (@0) = 0.1. This means that the rectifying ability of the diode during 
operation with narrow pulses of forward current (6g < 1) gets considerably worse. 


CONCLUSIONS 


1. Pulse parameters which characterize switching of the diode from the froward to 
the reverse direction, achieve virtually their steady-state values at a forward current 
pulse duration of tg > 0.8 Tp. 

2. During a forward current pulse, the establishment of the voltage at the p-n junction 
occurs during the time t<0. 02 Tp. t 

3. With narrow pulses of forward current (tf « Tp) the hole distribution in the base has 
an exponential character, wherein the coefficient of fie exponential term and the exponent 
are decreased by a factor of V7 /2V© in comparison with a steady state distribution. 

4. When measuring y it is necessary to introduce the coefficient € = 1 — 


oe Vo, taking into account the hole recombination in the diode base. 


‘ 5. The rectifying ability of the diode which operates in the switching regime, gets 
abruptly worse at a duration of the forward current pulse tg< Tp. 
The author expresses his gratitude to K.S. Rzhevkin and B.G. Mendelev for a discussion 
of the work and a number of valuable comments. 
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*This magnitude is analogous to the diode rectification factor (the ratio of forward 
to reverse current), which characterizes the diode rectifying properties in the static 
regime. 
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MOLECULAR FREQUENCY STANDARD WITH 
SUBTRACTION OF REFERENCE OSCILLATOR ERROR 


V.V. Grigor'yants and M. Ye. Zhabotinskiy 


The report presents a brief discussion of a number of circuits providing standard 
frequencies on the basis of a molecular oscillator. Calculations are performed for a 
frequency standard with subtraction of the error of a reference oscillator which provides 
stable frequencies in the centimeter, decimeter and meter wave ranges. Circuits are 
developed for investigation of frequency-phase stability of the obtained standard signal: 
(1) by means of an auxiliary molecular oscillator; (2) by means of an auxiliary quartz 
oscillator and frequency multiplier. 


The advent of the molecular oscillator [1, 2], opening wide prospects for the creation 
of highly stable frequency and time standards, called for the development of methods 
of transferring its stability into the range of radio frequencies and at the same time in- 
creasing its power. Various authors have recently proposed a large number of circuits 
permitting practical utilization of the high stability of the molecular oscillator. Disre- 
garding purely technical differences, these circuits may be divided into three main groups 
according to their principle of operation. 

I. Circuits for frequency comparison of the molecular oscillator and an auxiliary 
oscillator (usually using a quartz crystal) permitting highly accurate measurement of the 
departure of the latter from the desired frequency [8, 4, 9]. 

Il. Circuits for automatic control of the auxiliary oscillator (by means of the molecular 
oscillator), differing from the circuits of group I in the presence of a control loop which 
acts upon the auxiliary oscillator (automatic frequency or phase control) [5, 6, 7, 8]. 

Ill. A circuit for the creation of a standard frequency by subtraction of the instability 
of a reference oscillator (see Figures 1 and 2), proposed by the authors in 1955 [9]. 

This scheme possesses a number of advantages over those in the above groups: first, in 
the absence of feedback, which eliminates the problem of system stability as a whole; 
second, in the absence of influence on the quartz oscillator, which ensures a purer 
standard frequency spectrum. However, in its initial form this circuit had a number 

of significant disadvantages. 

1. Due to the considerable difficulty in obtaining local-oscillator power directly from 
a frequency multiplier, (see Figure 1) it is necessary to use an auxiliary local oscillator 
(see Figure 2) having sufficiently high frequency stability, determined by the passband of 
the i-f amplifier (100 kc). This leads to complication of the circuit and to additional phase 
modulation of the effective signal. 

2. For excellent filtering the sub- 
tracting sections must possess high 
equivalent Q. 

3. The circuit provides only one 
stable frequency. 


Figure 1 Figure 2 


In conjunction with these disadvantages it was undertaken to develop an improved version 
of the circuit with error subtraction. The problem was approached in three phases: first, 
to eliminate the auxiliary local oscillator; second, to improve filtering and the subtraction 
sections; third, to obtain frequencies in the centimeter, decimeter and meter wave ranges 
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simultaneously. The work resulted in the design of a frequency standard with subtraction 
of the reference oscillator error which meets the above requirements (see Figure 3). 


F,>(n-1) f mult, mult, 


Figure 3 


As is seen from the diagram, the problem of obtaining local-oscillator power is solved 
by automatic phase control of the klystron f; through a quartz oscillator fy and automatic 
phase control of klystron fy through klystron f; with use of the same oscillator fy as the 
reference oscillator. Use of the phase stabilization circuit mentioned permits subtraction 
of the error at three frequencies simultaneously: fo, fy andfy. In addition, construction 
and alignment of the automatic phase control circuit is considerably simplified by virtue 
of the considerable signal levels from the mixers (the mixer at frequency f; provides a 
signal of several tens of millivolts, the mixer at frequency fo provides a signal of several 
hundred millivolts). This automatic phase control circuit is equivalent to a high-power 
frequency multiplier, since phase drifts of the automatic phase control are small [10]. 

Let us examine the relationships which must be fulfilled in this circuit. 

For operation of the automatic phase control of f; through fy with the use of fy as the 
reference signal it is necessary that 


fr= (WE AVigh (1) 
It is shown below that in Equation (1) it is best to choose the minus sign: 
tee Ce Lia (2) 


Further, for automatic phase control of fo through f; let us choose the following frequency 
relationship: 

fp = mf, + mfq (3) 
or, using Equation (2), 


fy = (mn —m-+m) fg = mnfq. (4) 


Frequency fz is chosen in the region of the k-th subharmonic of the molecular oscillator, 
that is, 


hig =fo—kmnfg. (5) 


Let us now discuss subtraction of the error at frequencies fo, fy and f,. 
Subtraction at frequency fg. In this case the intermediate frequency of Equation (5) must 
be divided by the factor k. Then adding f2, we obtain 


JT, —kmnf f , ; 
2 : ies 2 mnfg + mnfg = fo (6) 
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In this case, a stable frequency is obtained in the centimeter wave range, for in practice, 
due to the necessity of obtaining from the fy oscillator a local-oscillator power in the 
range of the molecular oscillator, k=1, 2, 3. 

Subtraction at frequency fj. The division factor of ay _¢ must be equal tokm. Then 


adding f1, we obtain 


, Hence, it is evident that in order to obtain a stable frequency it is necessary to add fg: 


f Nias 
2 fy tig =. 


With m = 3 to 10 we obtain a stable frequency in the decimeter range. 
Subtraction at frequency fy. The division factor of fj_¢ must be equal tokmn. Then 


adding fg, we obtain 


i, —kmnfq ' to 
hkmn "1d ~~ Emn° 


With appropriate choice of n we obtain a stable frequency in the meter range. Itis 
now evident why the minus sign is chosen in Equation (1) and in Equation (3) fo - mfg. It 
is this choice of frequency ratios which has ensured a signal in the automatic phase control 
circuit equal to kmn f,, which in turn permits eliminating from the circuit of the standard, 
one frequency divider with successive division factors of k, m andn. Moreover, an 
additional advantage resulting from Equations (2) and (3) is the difference in the inter- 
mediate frequencies and, correspondingly, in the reference signals in the automatic 
phase control of fj through fy and in the circuit for automatic phase control of fg through 
f{, which decreases the influence of parasitic coupling in the circuit. For tube economy 
and the compensation of phase drifts of the frequency multiplier of the quartz oscillator 
it is desirable that the multiplier have a multiplication factor equal tom. Then from the 
point of the multiplier where the signal of frequency mf, exists a tap may be made to 
the phase detector of the automatic phase control circuit cf fo through fj; otherwise, 
frequency mfg must be specially created. 

We shall now discuss the requirements for the-subtraction sections. For the sake 
of simplicity in evaluation we shall consider sideband filtering at the -3db level to be 
adequate. To the subtracting device (assuming this is an ordinary mixer) there is 
applied the signal from the divider with a frequency of fj-¢/m 1, where mj is the division 
factor, and a reference oscillator signal with frequency kg. The distance between the 
side frequencies is 2f; _¢/my, that is, the Q of the filter for attenuation of the side 
frequencies at the -3db level is (10) | 

mf q 
~ Bice 


Hence it is evident that in order to improve the filtering it is necessary to decrease the 
division factor, preceding subtraction, and to increase fj-_f. 
Equation (10) may be rewritten in another form if we consider that mjfy =f - fj-¢ and 
fj-f « fp, where fg is the frequency of the molecular oscillator, 
Vip se is 
C5 a if Vis 


(11) 


Consequently, in ordinary subtraction the required filter Q does not depend on the waveband 
at which subtraction is performed. Let us define this Q Equation (11), assuming fj-¢ = 50 Me: 
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that is, in the centimeter range where cavity resonators with Q of (3 to 5) - 104 are 
available, we may obtain side-frequency attenuations of 102 to 103 times. The Q of 
cavities in the decimeter range is lower, hence in filtering with a single tuned circuit 

we may obtain side-frequency attenuations of only 10 to 100 times. In the meter range 

it is difficult toensure the required filtering with a single tuned circuit, but this difficulty 
is overcome by the use of the simplest automatic phase control circuit. In addition, 

early in 1958S. P. Dubonosov proposed that subtraction of the reference oscillator error 
be performed in two successive steps, that is, actually to lower the value of my in Equa- 
tion (10). This method permitted satisfactory filtering of side frequencies in any waveband. 
Above, in the examination of error subtraction at frequency fj dual subtraction was per- 
formed automatically, otherwise it would not be possible to obtain stable frequency fo/k. In 
this manner filtering in this range was also facilitated (see Equations (8) and (7)]. 

It should be noted that dual subtraction of error is sometimes unsatisfactory. In such 
cases it is necessary to change to higher order subtraction. In this connection let us 
examine sequential subtraction of reference oscillator error of different orders up to the 
general case of the k-th order subtraction. 

Let us first discuss the relation between the frequency multiplication factor ](7 = kmn) 
in Equation (5) of the quartz reference oscillator (which yields the frequency of the molecular 
oscillator) and the division factors my, for various orders of error subtraction. In this 
case let us consider that fj_¢ = fo - Ify- 


(1) Single subtraction 


ie Uf 7 ho 
my iq = my,’ 
if J = my. 
(2) Double subtraction: 
i aiget vidi to 
are ae Mo alte Iq 7 mom, | 
if 7 =mj(1+ mg). 
(3) Triple subtraction: 
xi , tq Tq ‘ to 
mymzms, | mgm3 = oo fq  mymemg’ 
if 7 =m (1+ mg + magmas). 
(4) k-th subtraction: 
aad: Rate Tq nae Ba fi He 
MymMa...My, | Mp~Mpy...My | ins" Mp alg Iq T mys. arr), 7 
if 7 = mj(1+ m2+ mgmg...mx), 
that is, in the general case 
=O (Oh 77 2S ONO Se 3 5 oe UO e = o ON), (12) 
or with m ~ 10 
l= mym,... Mp. (12a) 


Let us now discuss the requirements imposed on the Q of the filters of the subtraction 
sections according to the subtraction order. As before, for the sake of simplicity let us 
assume that side-frequency division at the -3db level is adequate. 

(1) Single subtraction. 

Let us rewrite Equation (10), for fj_¢ = pfy where p is some coefficient 


| mig Gull my my 
et Diet, opig tneop! eee 
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(2) Double subtraction: 
Mot MMofq my 


(fet ia) 2fq (Pm) 2” 


(3) Triple subtraction: 
gan Tria MMM, f gq ms 
3 9 (ict j Iq yee 2h (Darn) ae 2h 
“\ myn, * me q) 


(4) k-th subtraction: 


myt 

Qn pe SRE I 
ene 
a MyM... M, : Wy ap q 
MM... Mprfq Mp, 


~ 2(p-+m,+...+mym,...m, 4) fq sen 9 (13) 

In the above simplifications it was assumed that p ~1 and all m~ 10. If needed, the 
precise equation (13) may be used (up to the ~ sign), particularly in those cases where 
fj-¢ and fy are widely different in value (that is, p > 1 or p <1). 

In the specific calculations of any scheme for a frequency standard based on the 
principle of subtraction of the reference oscillator error it is necessary to assign two 
initial values: (1) the level of side frequencies in the output signal, which permits 
determining Q;, and my; from Equation (13); (2) the frequency of the output signal, know- 
ing which, we may approximately determine the value of 7. 

Further, knowing ] and my, and assuming in Equation (12a) that my = mg =...= Mg, 
let us determine the required error subtraction order from the formula 

l= mk". (12b) 


After this it is necessary to refine the value of J and all m, from Equation (12) and the 
value of fy from the formula fj_¢ = fg - I fg. 


METHOD FOR MEASURING THE STABILITY OF FREQUENCY 
AND PHASE OF THE STANDARD SIGNAL 


The usual method of studying frequency standards is by comparison of several identical 
circuits (the so-called internal comparison). 


standard 
standard under test 


aux molecular mnderitest 


oscillator 


Figure 4 Figure 5 


Below we examine a circuit for comparison of the output signal with the auxiliary 
molecular oscillator. 
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The circuit develops a standard signal with frequency iipeys n= foot io studying this 
signal it is most simple, multiplying it n times, to compare its n-th harmonic fon = 49) 
with the signal of a second molecular oscillator fj. In this case we obtain the difference 


frequency “ ; 
Q = h = Jo: 


Included in Q*, as also in fo; are the phase errors of the frequency standard circuit as 
well as the phase errors introduced by the auxiliary frequency multiplier. 

For measurement of the frequency stability of the output signal of the investigated 
standard it suffices to compare * with the frequency of an auxiliary low-frequency 
oscillator. In order to evaluate stability with an error on the order of 10-10 the 
auxiliary oscillator may have an error on the order of 10-2 -10-3. 

In order to evaluate phase fluctuations in f§ it is, moreover, necessary to shift 
the signals of both molecular oscillators. In this case we ‘obtain =f, - fo. 

Applying 9* and © to the phase detector, we may measure the phase fluctuations of 
the frequency-standard circuit. 

The block diagram of a system based on the above considerations is shown in Figure 4. 
The stability of the auxiliary local oscillator must be such that it is desirable to employ 
quartz-controlled automatic phase control. Otherwise, extra phase modulation of Q” and 
Q is possible. 

As seen from Figure 3, the standard-frequency circuit includes a local-oscillator 
klystron whose phase is quartz stabilized. If the power of this local oscillator at wave- 
length of 1.25 cm is sufficiently great, it may be used at the same time in the circuit for 
measuring the stability of the frequency standard. The problem is further simplified if 
we do not require balanced reception of the molecular oscillator signal. The resulting loss 
in receiver sensitivity is quite small. 


standard under test 


Figure 6 


Two versions of the circuit using a local-oscillator frequency standard and unbalanced 
mixers are shown in Figures 5 and 6. The circuit in Figure 5 uses two hybrid rings of 
ordinary design. The circuit in Figure 6 uses one hybrid ring with six inputs instead 
of the usual four. 

For measurement of the phase stability of the standard signal, in addition to the circuit 
described above, there was developed a circuit which generally does not require a molecular 
oscillator. A block diagram of this circuit is shown in Figure 7. The frequency of a 
quartz -controlled auxiliary oscillator (equal to the frequency f9/n of the output signal 
of the frequency standard with error subtraction) is multiplied by a factor n and applied 
to its input as a simulator of the molecular oscillator. Accuracy of the output frequency 
of the standard is equal to that of the auxiliary quartz oscillator. By applying these 
two frequencies to the phase detector we may measure the phase stability of the standard 
signal. The requirements for stability of the quartz oscillator (~ 10-7) are determined 
by the i-f bandwidth, taking into account the frequency multiplication factor of the 
quartz crystal. The molecular oscillator simulator is also quite convenient for 
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» 


prealignment of the standard-frequency cir- 
cuit. 

The circuit described above (Figure 7) standard under test 
was used for investigation of an experi- 
mental model of a frequency standard with 
subtraction which provided a stable fre- 
quency in the meter wave range. As was 
to be expected, the converted signal of 
the standard maintained the frequency 
stability of the molecular oscillator. 

Random phase modulation of the signal 
(measured by means of the circuit shown 
in Figure 7) proved to be on the order of 
one angular minute. 

The circuit of the frequency standard 
was tested in a setup with an independent 
molecular oscillator with a carbon ab- Figure 7 
sorption pump. The oscillator operated 
without interruption for more than 40 hours. 


| 
| 
| 
| 
| 
| 
| 
! 
| 
| 


REFERENCES 


1. N.G. Basov, Pribory i tekhnika eksperimenta, 1957, 1, 1-7. 
2. T.P. Gordon, H.Z. Leiger, S.N. Townes, Phys. Rev., 1954, 95, 282. 
3 


A. Ya. Leykin, Report at the Scientific Session of the A.S. Popov Society, May 1957, 


Izmeritel'naya tekhnika, 1959, 7, 41. 


4, N.G. Basov, I.D. Murin, A.P. Petrov, A.M. Prokhorov, I. V. Shtranikh, Izv. vuzov 


MVO SSSR (Radiofizika), 1958, 2, 3, 50. 


5. I.D. Murin, Izv. vuzov MVO SSSR (Radiotekhnika), 1957, 1, 5, 555. 
6 


© CO 


10. 


Institute of Radio Engineering and 


G.A. Vasneva, B.A. Gaygerov, V.V. Grigor'yants, G.A. Yelkin, M. Ye., Zhabotin- 
skiy, Radiotekhnika i elektronika, 1957, 2, 10, 1300. 

E.T. Davis, Phase stabilization to microwave frequency standards, Wescon reports, 
August 1957. 

I.L. Bershteyn, Dokl. na konferentsii po statisticheskoy fizike, Gorkiy, October 1958. 

G.A. Vasneva, V.V. Grigor'yants, M. Ye. Zhabotinskiy, Clock-maser circuits, Dokl. 
na 2-y Vsesoyuznoy konferentsii MVO peat po radioelektronike, Saratov, September 
oo 

I. L. Bershteyn, Report at the Scientific Session of the A.S. Popov Society, May 1957. 
Radiotekhnika i elektronika, 1958, 3, 2, 288; I.L. Bershteyn, V.L. Sibiryakov, 
Radiotekhnika i elektronika, 1958, 3, 2, 290. 


; : 0 
Electronics AN SSSR Submitted to the editors 1 March 196 


285 


BRIEF COMMUNICATIONS 


WIDEBAND AMPLIFICATION 
WITH FREQUENCY CONVERSION 


L.K. Mikhaylovskiy 


The author has experimentally observed the interaction of two microwave signals in a 
nonlinear element, which interaction occurs at any arbitrary relationship between the signal 
frequencies.* This permitted proposing the creation of wideband microwave amplifiers 
based on a new method of amplification [1]. This method is based on control of the ampli- 
fied signal by the parameters of the nonlinear element, ** which parameters, in turn, cause 
a change in value of the complex reflection coefficient from the given element for the auxiliary 
(supply) signal. In addition, as has been shown by theory and experiment, it is possible to 
achieve amplification of the microwave envelope upon conversion of the carrier frequency .*** 

The essentials of the proposed method are as follows. Let two microwave signals be 
applied simultaneously to a given nonlinear element. One of the signals is a control signal 
whose envelope is to be amplified and the other is a continuous (or pulse) signal from an 
auxiliary generator. 

Let the first signal have a carrier frequency fe and mean (pulse) power Pe and let the 
second signal, correspondingly, have frequency fg and power Pg. The signal to be amplified, 
in acting on the nonlinear element (e.g., a crystal detector in a detector head or a magnet- 
ized ferrite within a waveguide), causes a corresponding change in the capacitance and re- 
sistance of the crystal detector and the occurrence of precession of the magnetic moments of 
electron spins in the ferrite. Thus, the control signal is first, as it were, detected by the 
substance of the nonlinear element and then, as a consequence, causes a change in the ex- 
ternal electromagnetic parameters of this nonlinear element for the auxiliary signal. The 
change in internal microparameters occurs in accordance with the change in envelope of the 
control signal but with a definite relaxation time which, in combination with the Q of the 
internal circuit of the nonlinear element, determines the amount of distortion of the envelope 
after its "transfer" to a new carrier. 

The auxiliary signal, of course, also causes corresponding changes in the parameters 
of the nonlinear element and hence, during the action of the first signal, it should possess 
constant power (Pg appl. = const or Pg pulse appl. = const). Hence it follows that the 
operating conditions of the nonlinear element must be so chosen that with a change in value 
of the auxiliary signal power reflected from the nonlinear element the parameters of this 
element must remain practically constant. 

As a result of the above process it develops that, after the arrival of the control signal 
at the nonlinear element, amplitude and phase modulation corresponding to the envelope of 
the amplified signal appear in the reflected part of the auxiliary signal or that part of the 
auxiliary signal which has passed through the element. 


*This effect, detected at microwavelengths in an artificial medium (a solid body), prob- 
ably has much in common with the Luxemburg-Gor'kov effect observed at longer wavelengths 
in a natural medium — the ionosphere. 

**By a nonlinear element is meant any device or substance changing its electromagnetic 
parameters under the influence of electromagnetic waves. 

***Here and henceforth we refer to power amplification of the signal, for voltage ampli- 
fication (or current amplification) has no significance in microwave engineering. 
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Mathematically this process may be described as follows. With Po = 0, let the auxiliary 
signal power reflected from the nonlinear element be py (fgP) Pg and with Pg # 0 let it be 
P2 (fefgPcPg) Pg, where p, and 2 are the corresponding focal oF the complex reflection co- 
efficients for the auxiliary signal, depending both on frequencies fe and fg and on the powers 
=r. 
Pe and Pg [oO = oe ]. Moreover, it is evident that the difference 


Tig aed Pk =f a\ayl 


5 =AP, (1) 


8 
will determine the value of the change in power of the auxiliary signal reflected from the non- 
linear element due to the action of the amplified signal on this element. The gain in this case 
will be 


BPalfofgPoPe) 4 Pg 


Po Po 


(2) 


If the nonlinear element permits the application of sufficient power Pg that even at extremely 
small values of Ap it develops that Ap Pg >Pe, there will then occur amplification of the signal 
envelope of fe with conversion of its carrier to the frequency fg. 

It is evident that upon the action of a signal with power Pg (with Pc = 0) on the nonlinear 
element the signal reflected from this element will contain components with frequencies of 
kfg, where k=1, 2, 3, etc. If, in addition, this element is acted upon by a signal with power 
Pe, then the signal reflected from it will contain combination frequencies upon which there 
will be expended a fraction of power APg (with p2 <p). The remaining part of APg will go 
into the increase in losses in the nonlinear element and the possible creation of oscillations 
of those frequencies for which the element, under the action of signal Pe, will become a nega- 
tive resistance, the modulus of which is equal to or greater than the characteristic impedance 
of the waveguide. In the case where the power of these oscillations exceeds the value of Pe 
there also arises the effect of amplification of the signal envelope with conversion of its car- 
rier frequency. 

With p2 > p; the power increment AP» as well as the combination frequencies and the fre- 
quencies from the negative-resistance effect may obviously arise only due to a decrease in the 
power losses of the auxiliary signal in the nonlinear element and due to a decrease in the power 
of its harmonics kfg (except power Pc). 

The described method of amplification is based on the change in value of the modulus of 
the corresponding reflection coefficient. It is evident that, in principle, equal success may 
be obtained by using the change in phase of the reflection coefficient and the change (under 
the influence of the control signal) of the corresponding transfer constant. In the latter case 
the nonlinear element will be a four-terminal network [2]. 

For an increase in k, nonlinear elements similar to those described above may be con- 
nected in parallel or in series; for an increase in the signal-to-noise ratio the auxiliary 
generator must be pulse triggered and automatically follow the signal to be amplified. 

The frequency ranges of signals interacting by the method described above are deter- 
mined by the properties of the nonlinear element and, in principle, may lie in any portion 
of the electromagnetic spectrum. 

By way of example let us point out that the electrically and magnetically controlled micro- 
wave attenuators and phase inverters which have recently come into wide use are based on 
the effect of a change in the parameters of the substance of the nonlinear element under in- 
fluence of the envelope of a control signal, with the frequency of this signal fe = 0 and the 
frequency of the controlled signal fg lying in the centimeter range.* If we examine the other 
extreme case with fg = 0, when, instead of the microwave supply signal, there is applied to 
the nonlinear element a supply signal in the form of a d-c voltage and, possibly, a permanent 
constant magnetic field, then we likewise obtain amplification if under these conditions and 
with appropriate tuning the negative input resistance of the nonlinear element is ensured. In 
other words, a reflection coefficient greater than unity will be ensured for the amplified 


*The possibility of creating electrically controlled (rather than mechanically controlled) 
oscillators was reported by the author in 1953-1954 at a seminar at Moscow State University 


conducted by S.D. Gvozdover. 
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signal. Moreover, amplification will occur without change in the carrier frequency. Also, 
the closer the absolute value of the negative input resistance to the value of the character- 
istic impedance of the waveguide, the greater the gain. In the limiting case (that is, in the 
case of equality of these two values) oscillation will occur. 

Cross modulation of the two microwave signals by the use of a nonlinear element was 
achieved under experimental conditions wherein the frequencies of both signals lay in the 
centimeter* and millimeter ranges and also with the signal of frequency fe in the milli- 
meter range and the signal of frequency fg in the centimeter range. 

In all of the above cases it was possible to distinguish cross modulation with amplifica- 
tion observed on an oscillograph screen. An experiment was also performed for the pur- 
pose of measuring k with both interacting signals in the centimeter range. The measure- 
ments were made with a thermistor power-measuring device. The value of k in using one 
nonlinear element under the specified conditions was approximately 40. 

I wish to take this opportunity to express my thanks to B.P. Pollak and V.F. Balakov 
for their assistance in the experiments at millimeter wavelengths and to A.V. Zaytsev for 
assistance in the last of the above experiments. 
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SUPPRESSION OF HARMONICS 
IN COAXIAL RESONATORS 


A.I. Yel'kind 


For resonators and filters with square-law tuning in the microwave range coaxial butter- 
fly circuits are used which have a near-linear dependence of resonant wavelength A; on the 
variable length Z of the resonator. A resonator of ideal shape, formed by the intersection 
of a coaxial cavity by two planes, has an absolutely linear characteristic but is usually not 
used due to the prevalence of harmonics. Upon tuning such a resonator to the fundamental 
frequency of the excitation signal, harmonic resonance may occur, which is inadmissible in 
most applications. In practical resonator designs the elimination of multiple frequencies 
and, consequently, of harmonic resonance is associated with one or another deviation from 
the ideal shape. Such deviation results in nonlinearity of the relationship A, = A, (Z) and from 
this point of view must be a minimal deviation. Hence it is of interest to determine the mini- 
mum geometric dimensions of the deviation eliminating harmonic resonance in a resonator 
tuned to the fundamental frequency of the signal. 

Since the natural functions and frequencies for a coaxial resonator of ideal shape are 
known, such a problem is solved most simply by the perturbation method [1]. We present 
below an example of the calculations for a widely used type of resonator (see figure) with a 


“In addition, a completely similar effect of cross modulation was established experi- 
mentally, even as a result of the second harmonic of the reflected auxiliary signal. 
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contacting shorting plunger. The deviation is introduced by the shape of the shorting plunger. 
The dimension h is usually small in comparison with R, which permits considering the re- 
sonator at significant values of / to be slightly deformed relative to the ideal resonator delim- 
ater! By planes aa' and bb’. The spectrum of the natural frequencies of a resonator of 

ength Z 


On = On 9 + do, (1) 


where 


Ong = 201) = nc |Z 


is the spectrum of the undeformed resonator; 
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are the eigenfunctions of the undeformed resonator (Ey) = Hj). The numerator is integrated 
over the deformation volume: for g from 0 to 27; for z from 0 to; for r from R, to R, + hy 
and from R»-h», to Rg — wherein, due to the smallness of h, we use E(r) ~ Ep) and E(r) = 

= E)R;,/R2. Integrating the denominator over the entire volume: 


where W is the characteristic impedance of the, coaxial resonator in ohms. 

In resonator operation the fundamental frequency of the signal coincides with the first 
natural frequency w;. For the frequency of the n-th harmonic of the signal nw; detuning 
occurs relative to the nearest natural frequency wy. From Equations (1)-(3) it follows that | 


na 


. Ata 5 
1 — ©, = n80, — 60, = (Ong / 20) f (R, h) (sin 2a | Z— — sin 2nstl | 2). (4) 


The equation n sin x - sin nx = 0, where x = 21 /Z, has no roots in the interval 0<x<q. 
This indicates that for all harmonics the detuning differs from zero at resonator lengths 

7 > 21, With Z = 21 harmonic resonance is possible, but this case is remote from the 
practical application of resonators. For complete elimination (suppression) of harmonics 
the detuning must extend beyond the half-bandwidth of the resonance curve wy /2Qn. Ina 
resonator tuned to the fundamental frequency of the signal, harmonic resonance can occur 
only at the higher harmonics (n half-wavelengths of z). In this case the Q for the harmonics 
is higher than for the fundamental frequency and increases withn. Qn ~ VnQj [1]. 

It follows from Equation (4) that detuning also increases with the harmonic number. 
Consequently, it suffices to determine the condition for suppression of the lowest (second) 
harmonic. The minimum detuning ensuring practically complete suppression of the second 
harmonic is represented by the value 


201 — @2 = @2/Q). (5) 


Such detuning exceeds by 2-3 times the half-bandwidth of the resonant curve w2/2Q, (Q)~1.4Q,) | 
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The dimensions !, h which yield detuning in accordance with Equation (5) or greater must 
satisfy the condition 


V 


h h, 1 \3 W 
(a+ B)(z) 2 mee? 2 

which is obtained by substituting Equation (4) into Equation (5) and expanding Equation (4) 
into a series in terms of 1/Z up to the third order; :andh satisfying Equation (6) are the 
minimum dimensions ensuring suppression of the second harmonic and, consequently, of 
any of the higher harmonics. With W = 50 ohms, Q; = 103 and h/R = 0.1 condition (6) is ful- 
filled with 1/Z > 0.04. With smaller !/Z or h/R suppression of harmonics is not ensured. 
Since the shape of the shorting plunger is usually not chosen for the purpose of the suppres~ 
sion of harmonics, the other end of the resonator (the bb' end) is given a specially graduated 
shape with the :and h dimensions calculated from (6). 

Taking into account the nonlinearity of A; = A4(Z), the tuning relationships for the reson- 
ator are obtained from Equations (1), (2) and (3), where n=\1, Ay = 27c/w4. 
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CALCULATION OF RIBBON STREAM ELECTRON GUNS 


L.E. Bakhrakh and A.D. Sulimin 


1. SHAPE OF THE ELECTRODES 


We will proceed from that part of the cylindrical diode where we use as a cathode the part 
of the internal side of the external cylinder bounded by a sector of angle 26 (Figure 1). In 
a cylindrical diode with a space-charge limited current, the behavior of the potential U(r) 
can be represented in the form 


— | Ua (1) 


where rg and ry are the anode and cathode radii respectively; r is the distance from the 
origin of the polar coordinates to the point under consideration in the interelectrode space; 
Ug is the anode potential; 6 is a function of the ratio r/r; determined by the expression 

Pes © 11 47 
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(eles 
"y 
In order to preserve the rectilinear motion of electrons in this sector (under the condi- 
tion that inside this sector Poisson's equation holds, but outside it, Laplace's equation 
holds), one must place outside the sector, electrodes of such a form that at the boundary of 
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Figure 1 Figure 2 


of the beam, the following condition is fulfilled: 


U Uy (2), 


j= 


/OU 
erhalree 


(i.e., at y = 0 in the x, y coordinate system (Figure 1), where U)(x) is the potential distri- 
bution corresponding to Equation (1). 
We seek a solution of Laplace's equation in series form 


(2) 


id Sie Cabs 3) 
0 
On the basis of boundary conditions (2) we have 


U (x; 0) = Uoi(@), 


(=) Hie ON 
dy v=0 


To find the other coefficients of the series we substitute Equation (3) into Laplace's equation, 
whence we obtain the recursion formula 


(4) 


Un (2) 
Ce) (x 22) (n+ 4) ’ 


by means of which together with Equation (4) all the coefficients of the series can be expressed 
in terms of the function Ug(x) and its derivatives. 

The coefficients calculated in this manner allow a determination of U(x, y) and con- 
sequently, the sought-for form of the electrodes. The result of similar calculations limited 
to three terms of the expansion is given in Figure 2. The electrode shape found in this 
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manner differs somewhat both from the electrode shape determined by Pierce [1] and from 
the calculationof Radley [2]. Figure 3 gives the behavior measured in an electrolytic tank 
of the potential along the boundary of the flow in the case of the electrode shape calculated 
above. One may see that the potential distribution obtained agrees satisfactorily with the 
theoretical (solid curve) corresponding to Equation (1). 


2. EFFECT OF THE ANODE APERTURE AND SPACE CHARGE ON THE FORM 
OF THE FLUX IN THE ANODE SPACE 


It is well-known that the space-charge limited current in a cylindrical diode of length 
is determined by expression 


14,64 - 10-82 
2p ("a 
Ma Ba e, 


where I, = /,// is the current per unit length along the axis of the cylinder in amperes per 
meter; Yr, is the anode radius in meters; Ug is the anode voltage in volts. Consequently the 
perveance of such a cylindrical gun with an angle of convergence 2 © (© in degrees) will be 


a 


(5) 


T, 14,64 - 10-6 (20) 
Of ry N30) 


A diaphragm with an aperture may be considered as a dispersing electron lens with a 
focal length determined from the relation 


where Eg is the electric field intensity near the anode determined by Equation (1). 
Using the lens equation / ; 
is 5 se ane 
where Sp the distance from the lens plane to the point at which continuations of the electron 
trajectories would converge in the absence of the dispersing force of the space charge, one 
may reach the conclusion [4] that the lens action of the anode aperture is converging only 
for the values rg/r, < 0.53. 

To explain the behavior of the electron stream in the anode space under the action of 
space charge, we will use the equation of motion of the outer electron. For a constricted 
stream half-width of ym and for a constricted stream-to-anode aperture separation of Zn 
this equation gives [3] 


y 
Ym = o(1-—*,-), 7) 
0 
2m — 2K sy’ (8) 


1b : : 
where K = Tan 2.4-104; yo is the half-width of the stream in the anode aperture. At K< 


a 
< y)/4s% the electron stream intersects the central plane, i.e., no constriction is formed 
and, consequently it is impossible to obtain a parallel electron stream at the output of the 
gun. Thus, 


RES, Ae (9 ) 
4s? 
is the condition for the presence of a constriction in the electron stream issuing from the 
gun. It is easy to see that this relation interrelates the quantities ra, 0, Y9/rk and P. 

The solid curves of Figure 4 illustrate the variation of Prag with ra/r_ at different 
values of the parameter 9 according to Equation (5). The dashed curves show the relation 
between these same quantities according to Equation (9), in which we set Yo=Yatge. The 
intersection of these curves at the same values of 6 gives a system of points lying on the 
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Figure 4 


boundary of the region of existence of the constricted stream. The curve joining these points 
is the lower boundary of the region of values rg/rk, Pra and © where it is possible to obtain 
a constricted stream. The upper boundary of this region is defined by the value rg/rk = 0.53. 
On the basis of Equations (5), (6), (7), and (8), Figure 5 shows the values of ym and 
Zp, as a function of 9 at several values of rg and P = 4.104 mye 
The data obtained may be used for an approximate calculation of the basic parameters of 
an electron gun for given characteristics electron stream. Let, for example, the following 
be given: current I anode voltage Uz, minimum stream thickness 2ym, its width ]7 and the 
permissible cathode current density. Then we determine the perveance of the flow accord- 
ing to the relation P =I,/] U3., we choose the anode radius rg and we determine Pr,. From 
the curve of Figure 4 we find for a give Prg the range of possible values © in the focusing 
regions and from Figure 5 we choose the value ©, which corresponds to a prescribed magni- 
tude of yy, and also we determine Zy,. Returning again to the curves of Figure 4 we find the 
sought-for relation rg/r, which corresponds to the resulting value of ©. We find yo =r, tg 9. 


Figure 5 


Electron guns having the electrode shape indicated, with an anode aperture transverse 
eross section of 1 x 8 mm, for 20 = 30°, rg/r_~ = 0.36, rg = 2 mm, have in fact yielded 
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ribbon electron streams in which the current I = 50 ma (Ua = 800-1000 v) when the current 
transmission coefficient was 90-95% and the operating conditions were close to the design 


conditions. 
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SOME SPECIAL FEATURES OF FIELD EMISSION 
CATHODE OPERATION IN A MICROWAVE FIELD 


M.I. Elinson and V.A. Gor'kov 


Field emission from a metal or a semiconductor can be obtained not only by the use of an 
electrostatic field, but also with electrical fields which change with time. * 

In 1958 we suggested a system in which the field emission electrons from the point were 
placed directly in the microwave fields of resonators or of high-power waveguides, and 
analyzed the possibilities of such systems [1] (see, for example, Figures land 2). 

Recently publications have appeared concerning the creation of a number of important 
types of electron devices based on similar concepts [2]. 

The goal of this paper is a qualitative theoretical analysis of certain features of the 
operation of field emission cathodes in microwave electrical fields, which were not consid- 
ered in Reference 2. 

1. It is obvious that the field emission of a cathode operating in a sinusoidal electric 
field, as a result of the strong nonlinearity of the volt-ampere characteristic, will occur 
in the form of transitory electron packets. 

One can show that ~ 42% of the charge emitted during a period is concentrated within 
the phase range 0.097 (16°). During this interval changes in the electric field amount to 
AE = +0.5%. Numerical integration of the one-dimensional (for simplicity) relativistic equa- 


tion of motion mz | V1 — (=) = E,sin ot. allows a calculation of the dimensions of the electron 


packets, their duration and space charge distribution. 

In the phase interval (0.4-0.6) 7 at Ey = 10’ v/em, d = 0.1 cm and w = 27-10! seeu!, 
the spread in the time of electron flight in the oscillator is At = 0.6%, and the velocity spread 
is Av = 0.4%. The packet extension at the resonator output is AJ= 3 mm and its.duration 
Ar = 107! sec [over the phase range (0.45-0.55)  AJ=1.5 mm and Ar = 5.107!2 sec]. 


*Therefore use of the term "electrostatic emission’ may sometimes lead to misunder- 
standing. 
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Figure 1 Figure 2 


Figure 1. Point in a toroidal microwave resonator: 1-resonator; 2- 

power input; 3-point with a heating device; 4-fine-mesh grid in the resonator 
wall; 5-equipotential space in which electron packets are propagated; dis 
the width of the resonator gap. 


Figure 2. Potential distribution in the oscillator gap with the point (solid 
curve) and without the point (dashed curve). At a=10v Ey) ~50E [3], 
where E is the uniform field in the resonator without a point. At E ~108 
v/em Ey) ~5°10' v/cm and is sufficient to produce high field emission cur- 
rent densities. 


p arbitrary units 
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Figure 3. Distribution of space charge density in an electron packet at the 
resonator output. 


0.67 Aj 
Figures 3-5 show the space charge distribution in the electron packets °= 4 2 at var- 
0 


4m 
ious distances from the resonator. The concentration of electron charge in the packet head is 
advantageous for radiating microwave systems. An electron beam in the form of packets can 
also be used in accelerators, klystron-type devices, etc. 
The concentration of electrons in the packet varies in the range from 10'! to 10!4 cm=’ 
depending upon the magnitude of the focusing magnetic field, the degree of ionic compensation 


(if used), and the microwave power. 
2. During operation of a field emission cathode in a microwave system, the character of 


its total bombardment changes substantially. 
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Assuming the angle of flight of electrons in the resonator to be negligibly small, and 
assuming that the ionization cross-section is independent of electron coordinates, in the 
nonrelativistic approximation (for ions), one may solve the problems of the number of 
ions incident per unit time on the cathode and of their energy. By elementary means one 
may obtain 

aoe S,d (1) 
UP DekT io: 
Here n; is the number of ions impinging on the cathode per sec; Igp is the time-average elec- 
tron current; p is the gas pressure; S; is the ionization cross-section; d is the gap width of 
the microwave resonator. It is easy to see that nj is half the number of ions incident on the 
cathode in the static case which is related to the specific character of ion motion (sum of 


the forward and oscillatory motion). 
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Figure 4. Space charge distribution in the electron packet 30 cm from the 


resonator. 
The ion velocity is 
a eee (cos Wtp — cos wt). (2) 
mo 


Here Ey is the microwave field amplitude; w is its angular frequency; wt) is the phase of the 
electron motion; mj is the ion mass. Taking into account that the maximum number of ions 
is formed when wty = 7/2, we obtain 


272 
ad Ops 


Gi = 2m,0 ~ 


(3) 


It is significant that the energy decreases with an increase of ion mass, since the heavy 
ions are particularly effective in the process of cathode sputtering, and that the energy also 
decreases sharply with an increase of field frequency. 

A more accurate calculation of the number of incident ions and their energy leads to the 
following expression: 


PS; Imax 0,44eE 
nw ne [0.294 + cat | (4) 
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Figure 5. Space charge distribution in the electron packet 100 cm from the 
resonator. 


where Imax is the maximum field emission current, 


2 F2 
Cd 


Dae 
2m, @ 


6; @ 1,25 (5) 


For mj = 100 my, Ep = 10’ v/em, w = 7.57 -10!° sec! (A = 8 mm) @,~10 ev; for mj = 
=50 my (argon) 6:~20 ev. It is apparent that the incident ion energy decreases radically 
in comparison with the static case and this example is close to the threshold of cathode 
sputtering. 

In such devices we virtually exclude the problem of cathode sputtering, and the stability 
of the field emission increases substantially. 
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CALCIUM ADSORPTION ON TUNGSTEN 


V.1. Makukha 


In connection with the wide use of various types of composite cathodes in modern vacuum- 
tube devices — particularly such cathodes on a base of alkali-earth metal which combine the 
properties of semiconductor and coated cathodes — great interest has been shown in the 
investigation of the behavior of films of alkali-earth metals on refractory metals, particularly 
on tungsten. However, up to the present only the Ba-W system has been investigated [1-8], 
whereas not a single report has been published on the Ca-W and Sr-W systems, this despite 
the fact that in many composite cathodes the addition of strontium or calcium compounds has 
resulted in considerable increase both in the stability of emission from the cathode and in 
the longevity of the cathode. Hence, investigation of the Ca-W system is of considerable 
interest in clarifying the role of calcium in the mechanism of increasing emission of compo- 
site cathodes. : 

The present paper discusses the results of an investigation of the Ca-W system by means 
of an field emission projector [8] with vacuum conditions on the order of (1-2)-107? mm Hg. 
Calcium beryllate in a tantalum boat was used as the calcium source. The calcium was 
sprayed in small portions on a monocrystal of tungsten from below and by careful heating, 
the individual portions were dispersed over the entire surface of the point. After dispersion 
of each portion the image was "frozen" and the current and voltage were measured and the 
emission pattern was photographed. 

The relative work function was determined from the slope of the characteristic [8] 


ley = £(1/V). 


1, EXPERIMENTAL RESULTS 


Curves for the dependence of the work function @ on the duration of calcium deposition 
ts were plotted for several tungsten points. Figure 1 shows one of these curves. For all 
the emitters the curves of @(tg) passed through a minimum of 2.4 + 0.1 ev. 


30 


———— ; 
) 5 10 15. 20 ts, min 


Figure 1. 9 as a function of duration of 
deposition of Ca on W. 


Figure 2 shows a series of patterns with successive spraying (with the point at room 
temperature) of individual portions of calcium with subsequent dispersion over the entire 
tungsten surface (at a point temperature not above 450° C). 

A thin layer of Ca sprayed on half of the point begins to migrate at 400° C and at 500° C 
spreads over the entire point in 30 seconds. The calcium avoids faces {110} and {100} . First 
faces {112} and {113} are coated (uniting in one bright spot) and the centers of faces {111} 
(Figure 2a) are also coated, In subsequent spraying the Ca spreads over face {116 » becomes 
firmly established on { 111} and evidently forms a narrower thick layer on faces {112 t, for the 
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latter emit not so strongly and then cease glowing altogether (Figure 2b). Figure 2c corre- 
sponds approximately to the optimum coating. The only emitting faces are {116} and {111} 
(the latter in the form of bright triangles). 

The method of ''supplementary spraying’ [9] permitted verification of the fact that at 
the given stage we are dealing with an emission pattern which is characterized by a consider- 
ably reduced work function at the optimum coating. The mean work function in this case was 
2.4 ev. Spraying of large amounts of Ca results in a work function attributable only to Ca 
(2.7-2.8 ev) as shown in Figure 2d, wherein it is evident that crystallite formation has begun 
For some pointsthe field emission pattern of a thick layer of Ca was slightly different at the 
same value of work function. 

The most typical stages of gradual evaporation of Ca from tungsten at 800° C are shown 
in Figure 3. After pre-heating for 15-30 sec Ca leaves faces {116} , appears as bright spots 
on faces {112} and continues to cling firmly to {111} (Figure 3a). With further heating faces 
{112} are freed and the triangles of the octahedron decrease in size but continue to emit 
strongly (Figure 3b). In Figure 3c faces {112} are completely free of calcium and faces 
{111} are also almost clear. Only the peripheries of these faces show strong emission. 
After 15-20 minutes of heating the calcium completely leaves the tungsten and the emission 
pattern is that of pure tungsten (Figure 3d). 


2. EVALUATION OF RESULTS 


The experiments revealed that migration of Ca begins at 400° C and evaporation at 700°C. 
In spraying Ca on tungsten with the latter at room temperature crystallites may form, which 
disappear at a temperature below 400°C. At none of the stages of spraying or evaporation 
were traces of calcium observed on faces {110}. In isolated cases of spraying indistinct annular 
borders are observed at the faces of the cube. 

In evaporation the Ca leaves the faces of the monocrystal in the following order: {110}, 
{100}, {112', {111} (our discussion is limited to the principal crystallographic axes). Con- 
sequently, the sequence of faces in order of decreasing adatom bonds of Ca will be 11) 12s 
{100} i {110} , that is, the strongest adsorption bonds of calcium occur at faces {111}. 

The use of the various representations [1] leads to the following results. 

1. The representation of "close contact'' of atoms explains the absence of adsorption of 
Ca on faces {110} and {100} but does not explain the presence of adsorption on faces {111 
and {112}. Actually, in representing the lattice of Ca atoms on these faces we do not obtain 
a pattern of close contact. 

2. Crystallogeometric consideration (further complicated by the requirement of preserving 
in the adsorbed film the lattice constant characteristic of the adsorbed substance in the crystal 
state) cannot serve to explain the adsorption of Ca on faces {111} and {112}. This is easily 
seen from the data of Table 3 in Reference 1. 

3. Proceeding from the approach developed by Drechsler [7] and using his plot for the 
energy of association at the tungsten crystal facets as a function of the values of the metal 
radius of the adatom p and of the sublayer atom r, we may obtain for Ca the following 
sequence of faces in order of decreasing association of a single atom: {123}, {112}, {112}, 
fit 1) , {100$, {110} — which also fails to explain the preferential adsorption of Ca at faces 
{111}. The adsorption at faces {112} may be explained by the fact that the association of a 
single atom of Ca in the "wells" of this face is greater than at the ''valleys'’ and in the "wells" 
of other faces of the tungsten monocrystal. Evidently, in the case of the adsorption of Ca 
Drechsler's single-atom theory may be applied to faces {112} with greater sucess than in 
the adsorption, for example, of Ba, which does not result in such a porous" film as does 
Ca. 

Comparison of the results of experiment with the crystallogeometric representations 
for the cases of adsorption of calcium and barium [1] is given in the table. 

The results permit explanation of the increase in emission of a barium film on tungsten 
(described in Reference 10) upon introducing calcium into it by filling with calcium those 
faces of the tungsten monocrystal on which the barium is not adsorbed. This phenomenon 
may also partially explain the role of calcium in composite cathodes on a barium base. 
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Figure 2. Emission patterns of calcium sprayed on tungsten. 


a - Spraying for 1.5 min, ¢= 3.4 ev; b - spraying for 4 min, 
¢ = 3.1 ev; c - spraying for 10 min, ¢= 2.4 ev; d - spraying 
for 25 min, = 2.7 ev. 


Figure 3. Evaporation of calcium from tungsten (at 800° C). 


a -After 15 sec, = 3.1 ev; b - after 40 sec, 9 = 3.4 ev; 
c - after 5 min, = 3.9 ev; d- after 15 min, g = 4.5 ev. 
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Comparison of Experimental Data with Crystallogeometric Theories of Adsorption 


Adsorption in the form of 
"coherent" films: faces in order 

(1) in compound theory, of decreasing 

(2) in" close contact" theory |single-atom bond 


Sequence of 


Lattice 
constant, A 


{123}, (112}, {141}, 


(1) not possible at faces 
{100}, {110} 


{110} and {100}, {112}, {111} 


(2) also not possible at 
these faces 


(1) not possible at faces 


: (123}, {112} {100}, | at faces | 
| {110}, {100}, {111} (1143, {140} fla and 
{100 


‘2) possible at faces {100}, isolated 


{111} 
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ADSORPTION OF STRONTIUM ON TUNGSTEN 


V.1. Makukha 


The method described in Reference 1 was used for an investigation of strontium on a 
tungsten monocrystal. 

For several tungsten points curves were plotted for the dependence of the work function 
@ on the duration of strontium spraying. One of these curves is shown in Figure 1. The 
curve 9(tg) passes through a minimum at ~ 2.2 ev. 
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a 5 10 1§ ts, min 

Figure 1. Dependence of ¢ on the 

duration of strontium spraying on 
tungsten. 


Figure 2 shows a series of patterns with sequential spraying of individual portion os 
strontium with subsequent spreading over the entire surface of the point. 

A thin layer of Sr sprayed on half the point begins to migrate at 400-500°C, First 
faces {233} are coated and then the periphery of the cube. In subsequent spraying the 
strontium is firmly established on {111} (Figure 2a). New portions of Sr lead to the pat- 
tern shown in Figure 2b, wherein emitting faces {111} and {2335 create an emission 
pattern in the form of triple lobes. With further spraying these lobes become bright right 
triangles at {111} . The photograph in Figure 2c corresponds to the approximately optimum 
coating. Faces {111} and {233} emit (the lobes are rotated, as it were, through 180° about 
an axis passing through the center of face {111} and the periphery of faces aint emit in the 
form of bright rings. The method of supplementary spraying [2] confirmed the presence of 
an optimum layer in the given case. The mean work function in this case was ~ 2.2 ev. 


After further spraying, the slope of the characteristic of lear = £(1/V) for a thick 


layer of strontium (see Figure 2d) yielded a work function of 9= 2.5-2.6 ev. 

The most typical stages of evaporation of Sr from thick coatings at 800° C are shown in 
Figure 3. After 10-15 sec of heating the Sr is firmly established on faces {233}, bright spots 
appear on {113}, the centers of the octahedron faces become dark (Figure 3a). Further 
heating results in the formation of bright triangles on {111} (Figure 3b), which are extremely 
stable in subsequent stages of evaporation. Subsequent evaporation frees the adatoms from 
faces {111} (Figure 3c). The image is characteristic for the stage preceding total evapora- 
tion of the strontium. After 20 min of heating at 900° C the Sr completely leaves the tungsten 
monocrystal and the point image acquires the appearance characteristic of pure tungsten 
(Figure 3c). 


EVALUATION OF RESULTS 


Migration of strontium begins at 450-500° C, evaporation at 750-800°C. During the 
spraying of strontium on a tungsten point at room temperature crystallites are formed which 
disappear at ~ 400°C. Traces of Sr on faces |110\ were not observed at any stage of 
spraying or evaporation. Around the faces of the cube the strontium is arranged in the form 
of a bright "cross.'"' There were observed cases of extremely short-lived coatings of faces 

:112; with strontium. The emission patterns at different stages of strontium evaporation 
are similar to the images at certain stages of spraying. In order of decreasing Sr adatom 
bonds the tungsten faces may be listed in the following sequence: /111;, {113}, {112} , 
{100} , {110}, that is, the firmest adsorption bonds of strontium occur at faces ld 

The use of certain representations [3] leads to the following results. 

1, Proceeding from Drechsler's point of view [4] and using the plots prepared by him, 
we may obtain for Sr the following sequence of faces in order of decreasing single-atom 
bonds: (123), {112}, {100}, {111) , {110}, which does not explain the preferential 
adsorption of Sr on faces {111} . 

2, Crystallogeometric considerations (further complicated by the requirement of pre- 
serving in the adsorbed film the lattice constant characteristic of the absorbed substance in 
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Figure 2. Emission patterns of strontium sprayed on tungsten 


a -After 1.5 min, *= 3.5 ev; b - after 2 min, 
c - after 3.5 min, ¢= 2.2 ev; d - after 7.5 min, 


a -After 10 sec, ©= 2.5 ev; b - after 40 sec, 
c - after 5 min, > = 4.2 ev; d - after 25 min, 
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o=3 eV; 
@= 2.5 ev. 


Comparison of Test Data with Crystallogeometric Theories of Adsorption 


° 


Adsorption in the form of Sequence of 

"coherent" films: faces in order 
(1) in compound theory, of decreasing 
(2) in" close contact'' theory |single-atom bond 


° 


eter, A 


Lattice 
constant, A 


Atom diam- 


(1) not possible at faces {123\ | 4912) eat face 

{110} , {112}5" {100} 5 {111}, {100}, 1141) mr {111} 
110} 

(2) possible almost solidly at less firml 

faces {100} and {111} at faces 


223) 


} 


the crystal state) cannot serve to explain the adsorption of Sr on faces {111} even in explain- 
ing the absence of the latter on faces {100}, {112}, {110}. This is easily seen from the 
data in table 3 in Reference 3. 
3. The representation of ''close contact" of atoms explains the presence of adsorption on 
faces {111} and the absence of the same on such fundamental crystallographic directions as 
{112} and {110} , but does not explain the absence of adsorption at the faces of the cube. 
Actually, in depicting the atomic lattice of strontium at faces {111\ we obtain their almost 
"'soild' contant (slightly less so than for Ba). Such close contact is obtained also at faces 
110 although Sr is not adsorbed at the latter. 
The results of the experiment and the application of crystallogeometric representations 
for the adsorption of Sr are given in the table (similar data for Ba and Ca are given in 
Reference 1). 
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DETERMINING THE FREE BARIUM CONTENT 
IN AN OXIDE CATHODE BY A RADIOACTIVE 
INDICATOR METHOD 


V.P. Vasil'yev and V.P. Pomazkov 


Of great interest in the physics of oxide-coated cathodes is the determination of the free 
barium content in the oxide layer. The methods developed to date have been based on the re- 
action of water vapor with free barium and the measurement of the pressure of the liberated 
hydrogen (the Berdennikov method) or the same reaction with nitrogen and other materials 
[1]. One of the recent experiments performed according to the Berdennikov method [2] had 
a sensitivity of 1-10-8g. 

The method proposed by us is based on the reaction of free barium with one of the halo- 
gens — iodine. The result of the reaction is a somewhat stable compound of Bal,. The 
radioactive isotope of iodine I'*! was used in the form of an aqueous solution without a carrier 
which permitted a considerable increase in the sensitivity of the method. 

The following procedure was used in purifying the iodine. First, utilizing the excellent 
solubility of iodine in benzene, several drops of aqueous solution of iodine were mixed with 
benzene and the latter was precipitated from the water. Here we determined the amount of 
stable iodine (the carrier) required to create the necessary vapor pressure in test vacuum 
tube. The solution of iodine in benzene was placed in a glass tube with 4-5 outlets and the 
iodine was subjected to vacuum distillation. The end of the tube in which the iodine finally 
condensed had a sealed capillary with a striker. Before the iodine was completely distilled 
the entire tube was carefully degassed. Then the tube with the radioactive iodine was welded 
onto the test vacuum tube. 

The vacuum tube was a diode with an end cathgde of diameter 5 mm and a plane anode. 
The material for the core was purest nickel (99.99%) annealed in vacuum. The face of the 
cathode was coated with a double carbonate: 55% BaCO3, 45% SrCOg on an organic binder. 
The thickness of the coating was approximately 100-120 microns. The cathode was activated 
under a vacuum of not more than 5:10-° mm Hg (that is, under ordinary plant conditions) in 
two stages. In the first stage the temperature reached 850° C and in the second stage 1050° C 
with testing of the emission current. In static operation the current reached 300-500 ma/cm?, 
The cathode was then cooled and the tube was sealed off at a vacuum of 10-’ mm Hg. 

Immediately after sealing off the tube the capillary was broken and the iodine vapor 
penetrated into the tube. The tube was placed in an oven at a temperature of 100-120° C for 
several hours, after which the oven temperature was raised to 200° C and the end of the 
ampule was cooled. As a result all the uncombined iodine collected at one point in the device. 
After this the tube was opened and the radioactivity of the cathode and anode was determined. 
If the distribution of radioactivity in the depth of the oxide layer was to be determined, the 
cathode was coated with paraffin and slices of 10-micron thickness were cut with a special 
microtome. 

The high specific activity of I! permits reliable measurement of its quantity up to 107'%g. 
From the amount of iodine in the oxide coating it is easytodetermine the amount of free barium. 
The average measured amount of free barium on the cathodes was 5-107'g, which corresponds 
to a concentration of approximately 0.05% relative to the weight of the entire layer. More- 
over, free barium was detected on the anode surface facing the cathode and was measured 
at 3-10°8¢, 

In order to determine iodine absorption in the oxide layer control tests were performed. 

For this purpose a tube was constructed with two oxide cathodes and two anodes, One of the 
cathodes was activated to the point of stable emission (as described above), after which 
oxygen was injected into the tube at a pressure of several atmospheres. After several hours the 
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tube was again evacuated and the second cathode was activated, after which the usual meas- 
surements with radioactive iodine were performed. The amount of free barium on the un- 
contaminated cathode was close to that given in the data above, while on the contaminated 
cathode it did not exceed 1-107? g and at the contaminated anode 5-107! g, 

The results show that absorption of iodine in the oxide layer is two orders of magnitude 
less than the amount of iodine associated with free barium in the layer. The data on the 
distribution of free barium in the oxide layer indicate an increase in its concentration from 
surface to core. The amount of free barium at the anodes is apparently smaller due to the 
fact that the vacuum conditions permit part of the free barium which evaporates from the 
cathode to combine with residual gases. As a result of the use of tagged atoms in deter- 
mining the concentration of free barium in an oxide-coated cathode it may be stated that this 
method is at least twice as sensitive as other existing methods. 
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